Skip to main content

Full text of "Anomalous superfluid density in quantum critical superconductors"

See other formats

Anomalous superfluid density in quantum 
critical superconductors 

K. Hashimoto * \ Y. Mizukami * , R. Katsumata * , H. Shishido * ^ M. Yamashita * ^ H. Ikeda * , Y. Matsuda * , 
J. A. Schlueter ^, J. D. Fletcher A. Carrington ' , D. Gnida *% D. Kaczorowski " , and T. Shibauchi * 

* Department of Physics, Kyoto University, Kyoto 606-8502, Japan, ^Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439, U.S.A., I'H. 
H. Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol, BS8 1TL, United Kingdom,* 'Institute of Low Temperature and Structure Research, 
Polish Academy of Sciences, P. O. Box 1410, 50-950 Wroclaw, Poland, * Present address: Institute for Materials Research, Tohoku University, Sendai 980-8577, 
Japan,* Present address: Department of Physics and Electronics, Osaka Prefecture University, Sakai, Osaka 599-8531, Japan, § Present address: RIKEN, Wako, 
Saitama 351-0198, Japan, and ttjo whom correspondence should be addressed. Email: 

Submitted to Proceedings of the National Academy of Sciences of the United States of America 

When a second-order magnetic ptiase transition is tuned to zero 
temperature by a non-thermal parameter, quantum fluctuations 
are critically enhanced, often leading to the emergence of un- 
conventional superconductivity. In these 'quantum critical' su- 
perconductors it has been widely reported that the normal-state 
properties above the superconducting transition temperature Tc 
often exhibit anomalous non-Fermi liquid behaviors and enhanced 
electron correlations. However, the effect of these strong critical 
fluctuations on the superconducting condensate below Tc is less 
well established. Here we report measurements of the magnetic 
penetration depth in heavy-fermion, iron-pnictide, and organic su- 
perconductors located close to antiferromagnetic quantum critical 
points showing that the superfluid density in these nodal super- 
conductors universally exhibit, unlike the expected T-linear de- 
pendence, an anomalous 3/2 power-law temperature dependence 
over a wide temperature range. We propose that this non-integer 
power-law can be explained if a strong renormalization of effec- 
tive Fermi velocity due to quantum fluctuations occurs only for 
momenta k close to the nodes in the superconducting energy gap 
A{k). We suggest that such 'nodal criticality' may have an impact 
on low-energy properties of quantum critical superconductors. 

superfluid stiffness | d-wave superconductivity | spin fluctuations | mass 
enhancement | quasiparticle scattering 

The physics of materials located close to a quantum critical point 
(QCP) is an important issue because the critical fluctuations as- 
sociated with this point may produce unconventional high tempera- 
ture superconductivity |1. 2|. Quantum oscillations |j3]|4] and spe- 
cific heat measurements |5 | have shown that, in some systems, as the 
material is tuned towards the QCP by controlling an external param- 
eter such as doping, pressure, or magnetic field, the effective mass 
strongly increases due to enhanced correlation effects. Along with 
this the temperature dependence of the resistivity shows a strong de- 
viation from the standard AT^ dependence in the Fermi liquid (FL) 
theory of metals, and often shows an anomalous T-linear behavior 
which corresponds to the A coefficient diverging as zero temperature 
is approached. 

Although there are many studies of non-FL behavior in the nor- 
mal metallic state 1 1 , 2 1, relatively little is known about how the QCP 
affects the superconducting properties below the critical temperature 
Tc. The superconducting dome often develops around the putative 
QCP so that when the temperature is lowered below T^, the super- 
conducting order parameter starts to develop and the Fermi surface 
becomes gapped. It is therefore natural to consider that the low- 
energy quantum critical fluctuations are quenched by the formation of 
the superconducting gap A, which means that the system avoids the 
anomalous singularities associated with the QCP. Perhaps because of 
this reasoning the superconducting properties are usually analyzed by 
the conventional theory without including temperature/field depen- 
dent renormalization effects resulting from the proximity to the QCP. 
For example, in Refs. |6 7| the NMR relaxation rate l/Ti in the su- 
perconducting state is fitted to the temperature dependence expected 

for particular gap functions with the assumption that the normal-state 
1/TiT is virtually temperature independent below Tc even when it 
has strong temperature dependence above Tc due to the magnetic 

In superconductors near the QCP, the electron pairing is of- 
ten unconventional with a superconducting energy gap A(fc) which 
changes sign on different parts of the Fermi surface ||8] |9]- This 
sign change stems from a repulsive pairing interaction for example 
resulting from antiferromagnetic spin fluctuations. In many cases, 
this leads to the presence of nodes in the gap A(fe) where the gap 
changes sign. The low-energy excitations from the ground state in 
these superconductors are governed by these nodal regions. In such 
nodal superconductors, the effect of quantum critical fluctuations on 
the excited quasiparticles should be k dependent. Then the question 
arises as to how this effect modifies the low-energy properties in the 
superconducting state. 

The penetration depth A(r) is a fundamental property of the su- 
perconducting state which parameterizes the ability of a supercon- 
ductor to screen an applied field by the diamagnetic response of the 
superconducting condensate. As fermionic quasiparticles are ther- 
mally excited from the condensate a paramagnetic cuiTent is created 
which reduces the screening and increases A. So measurements of 
A(r) give direct information about density and Fermi velocity of 
these quasiparticles |I0|. When the effective mass is enhanced by 
the quantum critical fluctuations, the effective Fermi velocity is ex- 
pected to be suppressed accordingly. In a one-component Galilean 
invariant superfluid, electron correlation effects may not cause the 
renormalization in the low-temperature penetration depth |11|. In 
superconducting materials, however, strong electron correlations do 
affect the renormalization resulting in an enhanced penetration depth, 
which has been reported both theoretically |12 13] and experimen- 
tally 1 141 ITsi . To discuss the energy-dependent effect of quantum 
criticality on superconducting quasiparticles, the temperature depen- 
dence of penetration depth at low temperatures is thus of particular 

Reserved for Publication Footnotes 0.1 073/pnas.1 221 9761 1 

PNAS I Issue Date | Volume | Issue Number | 1-9 


Penetration Depth. Here we begin by presenting results for the 
iieavy-fermion system Ce„nna,i+2 which is located close to a QCP 
(T is the transition metal element, n is the number of Celns layers 
alternating with the Tln2 blocks along the c axis). The most studied 
member of this series is the n — 1 member CeCoIns (T^ = 2.3 K) 
IT6I . in which clear evidence for non-FL behaviors in the normal 
state I17l|18 !| and nodal superconductivity has been found ||6lll9ll20l . 
The recent discovery of superconductivity at ambient pressure in 
the n = 2 member Ce2PdIn8 (Tc = 0.68 K) (21], which exhibits 
very similar non-FL properties [22, 23] |M] |25] [26, 27] to those in 
CeCoIns, allows a detailed comparison of the superconducting prop- 
erties to extract common features in these superconductors near the 
antiferromagnetic QCP. The magnetic penetration depth measure- 
ments have been made for high-quality single crystals of Ce2PdIn8 
and CeCoIns down to ~ 60 mK by using a tunnel diode oscillator 
with the resonant frequency of ~ 13 MHz (see 5/ Text: SI 1). The 
weak ac field is applied along the c axis, which excites supercurrents 
in the ab plane. 

The temperature dependent change in the in-plane penetration 
depth AA(r) = X{T) - A(T = 0) in both CeaPdIng and CeCoIng 
(Fig. lA) exhibits strong temperature variation at low-temperatures, 
much steeper than the flat exponential dependence expected for a 
fully gapped superconductor. The AA(T) data for CeaPdIng are re- 
producible in different crystals and the data for CeCoIns are fully 
consistent with the previous studies |28. 29, 30 1. The strong temper- 
ature dependence indicates substantial excitations of quasiparticles 
at low energies, evidencing the presence of line nodes in the energy 
gap. This is consistent with the residual density of states (DOS) in 
the low temperature limit observed by thermal conductivity 1 19 23), 
specific heat 1201 1241 . and NMR measurements |6 27 |. In particular, 
a d^2_y2 order parameter with nodes along the (110) directions has 
been established in CeCoIns by angle-resolved thermal conductiv- 
ity 1 19 1 and specific heat measurements [20] . The strong similarity 
between CeaPdIng and CeCoIns found in the low-temperature A(r) 
points to common nodal structure in these superconductors. 

In a pure d-wave superconductor with line nodes such as high- 
Tc cuprates, it is well-established that AA(T) shows a linear T- 
dependence at low temperatures, which stems from the linear energy 
dependence of low-energy DOS of quasiparticles. However, a clear 
deviation from the T-linear dependence is observed in the present 
heavy-fermion superconductors (Fig. lA). The data also strongly de- 
viate from the dependence expected for the dirty limit case 

I3II (Fig. 16). We rather find that the power-law dependence T°' 
with an unusual exponent a — 3/2 can describe the observed 
low-temperature variation in a wide range in both superconductors 
(Fig. IQ. 

A few explanations for the super-linear temperature dependence 
of AA(r) in d-wave superconductors have been put forth, including 
the effect of impurity scattering 131], non-local effect near the nodes 
|32l and phase fluctuations |33 |. Often an interpolation formula de- 
scribing a crossover from T to dependence AA(T) oc T'^/{T + 
T* ) is used to describe the experimental data. In particular, the im- 
purity effect leads to the crossover temperature T* ~ 0.83\/rAo, 
where F is the impurity scattering rate and Ao is the maximum gap. 
This successfully accounts for the systematic change of AA(r) with 
impurity scattering observed in Zn-doped (Ref. 1 34 1). 
In the present heavy-fermion case, however, fitting to this crossover 
formula in these relatively clean superconductors yields T* /Tc val- 
ues (~ 0.5 for CeaPdIng and ~ 0.2 for CeCoIns) substantially larger 
than the estimates from these theories (see 5/ Text: SI2 and Ref. 1 30|). 
Therefore, the T'^/^ dependence in a wide T/Tc range commonly ob- 
served in these superconductors with quite different Tc, which is dis- 
tinctly different from the T-linear dependence in e.g. 
(Ref. 1 34 1), rather suggests some inherent mechanism related to their 
closeness to the antiferromagnetic QCP. 

Superfluid Density. To discuss the precise temperature evolution of 
quasiparticle excitations, we analyze our data in terms of the super- 
fluid density 1/A^(T). We use the reported values of A(0) (280 nm 
for CeCoIns |30| and 1010 nm for CeaPdlng |35|), from which the 
normalized superfluid density paiT) = A^(0)/A^(T) has been ob- 
tained (Fig. 2A). The factor of ~ 3.5 difference in the slope of 
d\/d{T/Tc) (Fig. lA) is consistent with the difference in A(0) in 
these two superconductors, which is also consistent with the larger 7 
value in CeaPdIng |24|. This results in an almost collapse of the full 
temperature dependence psiT) into a single curve (Fig. 2A, inset), 
and the low-temperature variation shows the (T/Tc)^''^ dependence 
with nearly identical slopes. 

This (T/Tc)^/^ dependence can be also seen in other classes 
of materials which are close to antiferromagnetic order. In 
Fig. 2B we show data for the organic superconductor k-(BEDT- 
TTF)2Cu[N(CN)2]Br which is consistent with but measured to lower 
temperature than that in a previous report [|36 :: and the iron-pnictide 
superconductor BaFe2(Aso,7Po,3)2 1151 . Both these materials are 
known to have line nodes in their superconducting gap 1371 1381 . 

Fig. 1. Temperature dependence of the magnetic penetration deptti in heavy-fermion superconductors near the antiferromagnetic QCP. (A): Low-temperature changes 
in the magnetic penetration depth AA(T) = A(T) — A(T = 0) of single crystals of Ce2Pdlng and CeCoIns. The curves are shifted vertically for clarity and the 
data for CeCoIns are multiplied by 3.5. The inset shows the ac susceptibility over the whole temperature range measured by the frequency shift of the tunnel diode 
oscillator, showing sharp superconducting transitions. The dashed line is an exponential temperature dependence expected for a fully gapped s-wave superconductor. 
(B): The same data plotted against (T/Tc)^. The solid line represents a T'^ dependence. (C): The same data plotted against (T/Tc)'^^^. The solid line represents 
a T^/^ dependence. 

2 I 

Footline Author 










0.025 0.050 0.075 0.100 

Fig. 2. Universal T^^^ dependence of superfluid density in unconventional superconductors in the vicinity of the antiferromagnetic order. (A): The normalized 
superfluld density ps as a function of (T/Tc)^^^ at low temperatures for Ce2Pdln8 and CeColns. The lines represent T^/^ dependence. The Inset shows 
the overall temperature dependence up to (T/Tc) = 1. (B): A similar plot for Iron-pnictide superconductor BaFe2(Aso.7Po.3)2 and organic superconductor 
K-(ET)2Cu[N(CN)2]Br. The solid lines are the fits to T^/^ dependence. The low-temperature data for BaFe2(Aso.7Po.3)2 are vertically shifted for clarity 

In the BaFe2(Asi_i:Pa;)2 series, there is clear evidence for the an- 
tiferromagnetic QCP being located at a:: ~ 0.30 |15|. In n- 
(ET)2Cu[N(CN)2]Br, although the proposed phase diagram suggests 
that the boundary between the superconducting and antifeiTomag- 
netic states is a first-order phase transition |39 1, the anomalous criti- 
cal exponent near the Mott critical end point |40 | suggests the pres- 
ence of strong antiferromagnetic quantum fluctuations. The normal- 
ized superfluid density ps(T) in these superconductors shows very 
similar [T /Tc)'^^'^ dependence at low temperatures with slight devi- 
ations at the lowest temperatures below T/Tc ~ 0.05 ((T/Tc)'*^^ < 

These results imply that the (T/Tc)^^'^ dependence of superfluid 
density in a wide T/Tc range is a robust property in unconventional 
superconductors, in which strong antiferromagnetic fluctuations are 
present (for comparisons between iron-pnictides and cuprates, see 5/ 
Text: S13). 


If the quantum fluctuations survive the Fermi-surface gapping, the 
effective mass is expected to diverge in the zero temperature limit. 
In such a case, the strong enhancement of mass leads to a reduction 
of superfluid density when approaching the zero-temperature limit. 
It has been recently found in iron-pnictides that the zero-temperature 
superfluid density 1/ (0) shows a strong reduction at the QCP | IS*], 
indicating the strong quantum critical fluctuations directly affect the 
superconducting condensate. However, the temperature dependence 
of the normalized superfluid density ps (T) actually continues to rise 
with decreasing temperature (Figs. 2A and B, insets). This suggests 
that the effect of temperature dependent quantum fluctuations enters 
as corrections in the temperature dependence of superfluid density at 
low temperatures. The anomalous non-integer power-law tempera- 
ture dependence of Ps{T) universally observed in quantum critical 
superconductors thus calls for further theoretical understanding. 

Nodal Quantum Critlcallty. Below we will propose a possible sce- 
nario that this universal behavior can be naturally explained by in- 
voking a strong momentum dependence of renormalization due to 
the nodal gap structure. Suppose that the gap formation below Tc 
quenches the low-energy quantum critical fluctuations, which pre- 
vents the effective mass enhancement in the superconducting state. 
Then in the nodal superconductors, because of the strong momentum 
dependence of the gap magnitude we may consider that the degree of 

quenching of the quantum fluctuations in the superconducting state 
has a strong momentum dependence as well; the low-energy fluctu- 
ations are expected to be strongest near the nodes where the FeiTni 
surface is not gapped (Fig. 3A). To model this effect, we consider the 
angle dependence of the effective Fermi velocity v*p ((j)) along the un- 
derlying FeiTni surface. We consider a simple cylindrical Fermi sur- 
face and assume that the renormalization in Vp{k)/vp is inversely 
related to the enhancement in the effective mass m* (k) (Fig. 36). 
(Here 1/vp is given by the dynamic effective mass which is different 
from, but closely related to the thermodynamic mass II12I ). The ef- 
fective mass enhancement on approaching the QCP can be described 
by m*^ oc (p — Pqcp)~^, where p is a nonthermal parameter con- 
trolling the distance from the QCP at pqcp- The critical exponent 
P has been estimated experimentally by using magnetic fields as the 
parameter p, and in several materials a value of /3 close to unity has 
been reported |41) |42j. In the present case, we take the magnitude 
of the superconducting gap | A| as the control parameter, because the 
gap magnitude determines the degree of quenching of low-energy 
fluctuations. We thus assume vp{k) oc j A(fc)j'^^^, which allows us 
to calculate the temperature dependence of superfluid density by the 
integral over the Fermi surface S as 1 10, 13] 

where Vp (vp) is the effective velocity in the superconducting 
(normal) state, the subscripts i,j denote the directions of the cur- 
rent and vector potential (we take both along a), Y{k,T) — 
-2 — = dE is Yosida function, and f(E) is the 

Fermi-Dirac function for the quasiparticle energy E. By using the 
d^2_y2 formula of A(fc) — Ao cos(2(^) and l3 — 1, the normal- 
ized superfluid density ps is calculated with Vp{(l>) oc | cos(20)|^''^ 
(Fig. 3C), which results in the y^/^ dependence at low temperatures 
(Fig. 3D). 

As real materials will never be located exactly at the QCP there 
will be a cutoff for the diverging m* near the nodes. This can be mod- 
elled by introducing a minimum value for Vp. This leads to an up- 
ward deviation from the t'-^^^ dependence of ps at very low tempera- 
tures approaching T-linear behavior at sufficiently low temperatures 
(Fig. 3D). This can explain the essential features of the experimental 
observations. We note that the presence of the cutoff can be also ex- 
pected even at the exact QCP, because the dynamical susceptibility in 
the zero temperature limit should diverge only on certain regions of 

Footllne Author 

PNAS I Issue Date | Volume | Issue Number | 3 

Fig. 3. Nodal quantum criticality in unconventional superconductors. (A): The momentum dependent gap A(fc) (whose magnitude is illustrated by thin lines with 
grey shades) opens on the Fermi surface (thick line), and has nodes (red circles) at certain directions. In d^2_j,2-wave superconductors, for example, A(fc) has 
strong in-plane anisotropy Ao cos(20) as a function of azimuthal angle (f. In quasi-two dimensional systems, the Fermi surface is approximated by a cylinder, and 
thus the gap has nodal lines at perpendicular to the planes. At the nodes, the gap is zero and thus the quantum critical fluctuations may be present (red shades) on 
the ungapped Fermi surface. (B): The nodal quantum fluctuations lead to the momentum dependence of the renormalization in l/vp{k) near the nodes (blue lines). 
(C): The angle dependence of the renormalized Fermi velocity Vp{(f>) relative to the unrenormalized one vp along the Fermi surface, assumed for calculations of 
the superfluid density in (D). Near the nodes, we illustrate different cutoff levels, which model finite distances from the QCP or disorder. {D): Calculated normalized 
superfluid density as a function of (T/TcY/'^ with different cutoff levels, which explains the deviation from the T'^/'^ dependence at very low temperatures. The inset 
is the full temperature dependence up to (T/Tc) = 1. 

the Fermi surface (determined by the momentum dependence of the 
spin fluctuations), which in general may be different from the nodal 
points |9|. Of course, disorder would also produce some additional 
changes to the temperature dependence but we have not included this 
in the present model. 

An additional factor could also come from the temperature de- 
pendence of the renormalization of m* . As only the theiTnally ex- 
cited quasiparticles will be renormalized, the angular range near the 
nodes where this occurs (i.e., where A(fc) < ^bT) is quite limited 
at low temperature and so this has a rather minor effect. In our tem- 
perature range of interest, T < 0.2Tc, this region is actually limited 
to a narrow angle range of ~ ±3° near the nodes. If we add a cut- 
off constant of 0.3 covering this angle range (Fig. 3C), we found that 
the T-dependence is affected only in the lowest temperature range 
of [T/Tcf-^ < 0.01, above which the T^^^ dependence of ps still 
holds (Fig. 3Z)). This exercise implies that the inclusion of the tem- 
perature dependence of renormalization will not change ps (T) sig- 
nificantly in the temperature range of interest (for more discussion, 
see 5/ Text: SI4). 

We also note that in iron-pnictides the gap symmetry is most 
likely s-wave f43|, and the model based on the d-wave gap may 
not be applicable. However, the fundamental physics that the low- 
energy excitations are governed by the nodal regions should be es- 

sentially the same. Although the detailed structure of the momentum- 
dependent Fermi velocity (such as the precise value of the critical 
exponent /?) will affect the detailed Ps{T) (see 5/ Text: SIS), it is 
striking that such a simple model can capture the salient feature of 
the unusual behavior of quasiparticle excitation in the superconduct- 
ing state of these materials. In a FL theory, the renormalization of 
the effective FeiTni velocity vp in Eq. [1] can be described by the in- 
terplay between fc-dependent quasiparticle interaction and nodal gap 
structure 1 13 1. More detailed theoretical calculations will be needed 
to account for: strong energy dependence of critical fluctuations, the 
effects of disorder (such as quasiparticle scattering interference |44 |) 
and the inelastic quasiparticle scattering rate. 

We suggest that the nodal quantum criticality proposed here is 
an important aspect of unconventional superconductivity close to the 
magnetic QCP. Besides the penetration depth analyzed here there will 
be implications for most other superconducting properties such as 
thermal conductivity and the NMR relaxation rate l/Ti{T) which 
has been long known to have strong deviations from the T"^ law, 
which in the usual analysis gives overestimates of the residual quasi- 
particle DOS in quantum critical superconductors ||6]2]. It should be 
straightforward to extend our analysis to these other properties. 

4 I 

Footline Author 

ACKNOWLEDGMENTS. We thank A. Chubukov, H. Fukazawa, R. W. Giannetta, 
K. Kanoda, S. Kasahara, H. Kontani, S. K. Goh, P.J. Hirschfeld, A. H. Nevidom- 
skyy, T. Nomoto, R. Prozorov, I. Vekhter, Y. Yanase, and Y. F. Yang for discussions. 
Tliis work is supported by KAKENHI from JSPS, Grant-in-Aid for GCOE program 
"The Next Generation of Pfiysics, Spun from Universality and Emergence" from 

IVIEXT, Japan, EPSRC in the UK, the National Science Centre (Poland) under 
the research Grant No. 201 1/01 /B/ST3/04482 in Poland and by Argonne, a U.S. 
Department of Energy Office of Science laboratory, operated under Contract No. 
DE-AC02-06GH11357 in the USA. 

1 . See, for example, Varma CM, Nussinov Z, & van Saarloos W (2002) Singular or 
non-Fermi liquids. P/iys Rep 361:267-41 7. 

2. Gegenwart P, SI Q, & Steglich F (2008) Quantum criticality in iieavy-fermion metais. 
WafP/iys 4:186-197. 

3. Shishido H, Settai R, Harlma H, & OnukI Y (2005) A drastic change of the Fermi 
surface at a critical pressure in CeRhlns : dHvA study under pressure. J Phys Soc 
Jpn 74:1103-1106. 

4. Shishido H, etal. (2010) Evolution of the Fermi surface of BaFe2(Asi_a.Px)2 on 
entering the superconducting dome. P/iysRevl.ett104:057008. 

5. Stewart GR (2001) Non-Fermi-liquid Isehavior in d- and f-eiectron metals. Rev Mod 
Phys 73:797-855. 

6. Kohori Y, efa/. (2001) NIUIR and NQR studies of the heavy fermion superconductors 
CeTlns (T = Co and Ir). Phys /?ev 664:1 34526. 

7. Nakal Y, etal. (2010)-^^ Rand ^-^^As NMR evidence for a residual density of states at 
zero energy In superconducting BaFe2(Aso.67Po.33)2- Phys BevB 81 :020503(R). 

8. SIgrlst M, & Ueda, K (1991) Phenomenologlcal theory of unconventional super- 
conductivity. Rev Mod Phys 63:239-31 1 . 

9. Hirschfeld PJ, Korshunov MM, & Mazin II (2011) Gap symmetry and structure of 
Fe-based superconductors. Rep Pmg P/iys 74:124508. 

10. Chandrasekhar BS, EInzel D (1993) The superconducting penetration depth from 
the semlclassical model. Ann Phys 2:535-546. 

11. Leggett AJ (1965) Theory of a superfluid Fermi liquid. I. General formalism and 
static properties. Phys Rev140:A1869-A1888. 

12. Varma CM, Miyake K, & Schmitt-RInk S (1986) London penetration depth of heavy- 
fermion superconductors. Phys Rev LeH 57:626-628. 

13. Jujo T (2002) Effect of Umklapp scattering on magnetic field penetration depth In 
high-Te cupiates. J Phys Soc Jpn 71 :888-904. 

14. Gross F, el al. (1986) Anomalous temperature dependence of the magnetic field 
penetration depth in superconducting UBCis. ZPhys 664:175-188. 

15. Hashimoto K, etal. (2012) A sharp peak of the zero-temperature penetration depth 
at optimal composition In the iron-based superconductor BaFe2(ASi_^P^)2. Sc/- 
ence 336:1554-1557. 

16. Petrovic C, etal. (2001) Heavy-fermlon superconductivity in CeColn^ at 2.3 K. J 
Phys Condens /Wa(fer13:L337-L342. 

17. Nakajima Y, et al. (2007) Non-Fermi liquid behavior In the magnetotransport of 
CeMln5 (M : Co and Rh): Striking similarity between quasi two-dimensional heavy 
fermion and high-T^ cupiates. J Phys Soc Jpn 76:024703. 

18. Paglione J, ef a/. (2003) Field-Induced quantum critical point In CeColn -,. Phys Rev 
LeH 91:246405. 

19. Izawa K, et al. (2001) Angular position of nodes in the superconducting gap of 
quasi-2D heavy-fermion superconductor CeColn5. Phys Rev Lett 87:057002. 

20. An K. et al. (2010) Sign reversal of field-angle resolved heat capacity oscillations 
In a heavy fermion superconductor CeColn5 and <i^2 _ 2 pairing symmetry. Phys 
Rev Lett -104:037002. 

21 . Kaczorowski D, Pikul AP, Gnida D, & Tran VH (201 0) Emergence of a superconduc^ 
ing state from an antiferromagnetic phase In single crystals of the heavy fermion 
compound Ce2Pdlns. Phys Rev Lett103:027003(2009);/i»'(/. 104:059702. 

22. Kaczorowski D, Gnida D, Pikul AP, & Tran VH (2010) Heavy-fermion superconduc- 
tivity in Ce2Pdlng. Solid State Commun 150:411-414. 

23. Dong JK, etaL (2011) Field-induced quantum critical point and nodal supercon- 
ductivity in the heavy-fermion superconductor Ce2Pdlng. Phys RevXI :011010. 

24. Tokiwa Y, Gegenwart P, Gnida D, & Kaczorowski D (201 1 ) Quantum criticality near 
the upper critical field of Ce2Pdln8. Phys Rev B84:140507(R). 

25. I\1atusiak M, Gnida D, & Kaczorowski D (201 1 ) Quantum criticality in Ce2 Pding : A 
thermoelectric study. Phys RevS84:115110. 

26. Gnida D, IVIatuslak M, & Kaczorowski D (201 2) Anomalous magnetotransport in the 
heavy-fermlon superconductor Ce2 Pding . Phys Rev B 84:060508(R). 

27. Fukazawa H, Nagashima R, Shimatani S, Kohori Y, & Kaczorowski D (2012) Nodal 
superconductivity and non-Fermi-liquid behavior in Ce2Pdlng studied by ^^^in 
nuclear quadrupole resonance. Phys Rev B 86:094508. 

28. Ormeno RJ, Sibley A, Gough CE, Sebastian S, & Fisher IR (2002) Microwave con- 
ductivity and penetration depth in the heavy fermion superconductor CeColng. 
Phys Rev Lett 88:047005. 

29. Chia EEIVI et al. (2003) Nonlocallty and strong coupling in the heavy fermion su- 
perconductor CeColns: A penetration depth study. Phys Rev 6 67:014527. 

30. dzcan S, e( al. (2003) London penetration depth measurements of the heavy- 
fermion superconductor CeColn5 near a magnetic quantum critical point. Eum- 
phys Leff 62:412-418. 

31. Hirschfeld PJ & Goldenfeld N (1993) Effect of strong scattering on the low- 
temperature penetration depth of a cf-wave superconductor. Phys Rev B 48:4219- 

32. Kosztin I & Leggett AJ (1997) Nonlocal effects on the magnetic penetration depth 
In d-wave superconductors. Phys Rev Le(f 79:135-138. 

33. Chen 0, Kosztin I, & Levin K (2000) Unusual thermodynamlcal and transport sig- 
natures of the BCS to Bose-Elnstein crossover scenario below t^. Phys Rev Lett 

34. Bonn DA, et aL (1994) Comparison of the influence of Ni and Zn impurities on the 
electromagnetic properties of YBa2CU306.g5. Phys Rev 650:4051 -4063. 

35. Tran VH, Hillier AD, Adroja DT, & Kaczorowski D (201 2) Antiferromagnetic spin fluc- 
tuations In the heavy-fermion superconductor Ce2 Pding. Phys Rev 686:094525. 

36. Carrlngton A, ef al. (1999) Low-temperature penetration depth of n- 
(ET)2Cu[N(CN)2]Brand k-(ET)2Cu(NCS)2. Phys Rev LeH 83:4172-4175. 

37. Hashimoto K, ef al. (2010) Line nodes In the energy gap of superconducting 
BaFe2(Asi_a;Px)2 single crystals as seen via penetration depth and thermal 
conductivity. Phys Rev 6 81 :220501 (R). 

38. Taylor OJ, Carrlngton A, & Schlueter JA (2007) Specific-heat measurements of 
the gap structure of the organic superconductors K-(ET)2Cu[N(CN)2]Br and n- 
(ET)2Cu(NCS)2. Phys Rev Left 99:057001 . 

39. Kanoda K (2006) Metal-insulator transition in k-(ET)2X and (DCNQI)2M: Two con- 
trasting manifestation of electron correlation. J Phys Soc Jpn 75:051007. 

40. Kagawa F, Miyagawa K, & Kanoda K (2009) IVIagnetic IVIott criticality in a ^-type 
organic salt probed by NIVIR. A/at Phys 5:880-884. 

41. Gegenwart P, ef al. (2002) IVIagnetic-Fleld Induced Quantum Critical Point in 
YbRhjSis. Phys Rev Le/f 89:056402. 

42. Matsumoto Y, et a/. (201 1 ) Quantum criticality without tuning in the mixed valence 
compound /J-YbAIBj. Science 331:31 6-31 9. 

43. Yamashita lU, et a/. (2011) Nodal gap structure of superconducting 
BaFe2(Asi_xPx)2 from angle-resolved thermal conductivity in a magnetic field. 
Phys Rev 6 84:060507(R). 

44. Wang Q-H, & Lee D-H (2003) Quasiparticle scattering interference in high- 
temperature superconductors. Phys Rev B 67:02051 1 (R). 

Footline Author 

PNAS I Issue Date | Volume | Issue Number | 5 

Supporting Information 

Hashimoto et al. 10.1073/pnas.l221976110 
SI Text 

SI1-Sample information and experimental techniques 

High-quality single crystals of Ce2PdIn8 and CeCoIns were grown 
by the self-flux method IS1IIS2I . The typical lateral size of the crys- 
tals is 150 X 150 ^ra^ . For Ce2PdIn8 very thin crystals with a thick- 
ness of less than 10 were prepared to avoid contamination of 
Celns. We have checked that our crystals used in this study show 
no antiferromagnetic signals. The magnetic penetration depth mea- 
surements down to very low temperatures (~ 60 mK) have been per- 
formed by using a tunnel diode oscillator with the resonant frequency 
of 13 MHz [S3jlS4], which is mounted on a '^He-'^He dilution refrig- 
erator. The sample is placed on a sapphire rod, the other end of which 
is glued to a copper block on which a Ru02 thermometer is mounted. 
The sample is placed at the center of a solenoid which forms part of 
the resonant tank circuit. A weak ac field is applied along the c-axis 
so that the supercurrent flows within the afe-plane. Changes in the 
resonant frequency are directly proportional to changes in the mag- 
netic penetration depth, AA(T) = GAf{T). The calibration factor 
G is determined from the geometry of the sample |S5 1. 

Samples of K-(BEDT-TTF)2Cu[N(CN)2]Br were grown via the 
standard electrocrystalisation method |S6I. Measurements were also 
made using a tunnel diode oscillator. Special care was taken to cool 
the samples slowly across the temperature region 70 - 90 K where a 
structural ordering transition occurs. The cooling rate was restricted 
to 0.02 K/min, which is sufficiently slow enough to avoid suppression 
in Tc and reduction in electronic contribution to the heat capacity 


The data in Fig. 26 was taken for a sample with dimensions 
0.8 X 0.7 X 0.4 mm^. In these organic samples even weak radio 
frequency (RF) fields can produce self-heating effects at the lowest 
temperatures (T < 300 mK). To eliminate this effect, measurements 
were made with the sample at different positions on the axis of the 
RF coil, effectively reducing the applied field. Typically, a reduction 
in field of a factor three was sufficient to eliminate all signs of self- 
heating. This procedure was tested in three samples and consistent 
results were obtained. 

SI2-Comparisons with existing theories 

The origins of the deviation from the T-linear temperature depen- 
dence of A(r) in a d-wave superconductor are discussed in the frame- 
work of two different theories. Hirschfeld and Goldenfeld |S8 | have 
shown that a small amount of impurity causes a residual density of 
states near Ep, which can change the linear-T dependence (pure 
regime) to a quadratic T-dependence (impurity-dominated regime) 
at a crossover temperature T*^^. This characteristic temperature 
is given by T-*-^p ~ 0.83\/rAo, where F is the impurity scatter- 
ing rate and Ao is the maximum gap. On the other hand, Kosztin 
and Leggett | S9 1 have pointed out that since the coherence length ^ 
diverges near the nodes in a d-wave superconductor, Cooper pairs 
near the nodal directions, where ^ 3> A, cannot participate in the 
field screening. Therefore, the effective penetration depth should be 
larger than the local limit value below a characteristic temperature 
T^i^ ~ Ao^(0)/A(0), which leads to a quadratic T-dependence. 
Thus, for both theories the interpolation formula between two regions 
is given by AA(T) oc T"^ /{T + T*). Here we note an important dis- 
tinction between the temperature dependence of AA(T), which is 
directly measured, and the normalised superfluid density, Ps(T) = 
A^(0)/A^(r). The d-wave form for the superfluid density, ps{T) = 
1 - ^lT/T,, gives AA(T)/A(0) = \KT/T,) + Ip^'iT/T^f + ■ ■ ■, 


— «===J 1 1 






• CeCoIng 

• CesPdIng 



1 1 1 


0.05 0.10 0.15 0.20 0.25 


Fig. S1. Expanded view of normalized superfluid density plotted against T/Tc 
for Ce2Pdlng and CeColns. The solid lines are fit to the crossover behavior in a 
d-wave superconductor described in the text. 

which lead to a small quadratic component to A(r) that depends on 
^, where /j, is proportional to the inverse of the slope of the gap A((^) 
at the nodes, fj, = (4 In 2)fcsrc/(dA(0)/d0|nodc) 1S10|. Then, 
even in the pure d-wave model we have slightly concave tempera- 
ture dependence for A(T). Therefore, we take the following formula 
to describe the crossover behavior expected in these theories; 

p.(T) = i-M(rVTe)/(r + T*). [s-i] 

Figure SI shows the low-temperature superfluid density Ps(T) 
both for Ce2PdIn8 and CeCoIns. We carefully examine the possibil- 
ity of the two theories mentioned above by evaluating T*^^ from the 
present data. The obtained values of T^xp/Tc from the fits are 0.48 
and 0.22 for Ce2PdIn8 and CeCoIns, respectively. If we assume that 
the origin of nonlinearity of A(r) comes from the residual density of 
states induced by impurities, we can estimate the impurity scattering 
rate F via T* ^ 0.83\/rA^. The obtained value of F for Ce2PdIn8 
is O.I 1 K. By using this and Abrikosov-Gor'kov theory|Sl 1 1, we es- 
timated that the Tc reduction is 11% for Ce2PdIn8. However, recent 
specific heat measurements have suggested that the impurity level of 
Ce2PdIn8 is comparable to 0.22% Cd-doped CeCoIns 1S12|. In this 
case, Tc reduction is less than 5 %. Similarly, the obtained value 
of T* for CeCoIns using Eq. [S-1] is much larger than the value 
estimated from thermal transport measurements (< 30 mK) |S13|. 
Therefore, impurity effect are unlikely to explain the nonlinearity of 
A(r) in these systems. As for the nonlocal effects, the reported val- 
ues of ^0 ~ 8nm |S2 S14 S15| for both compounds give T^j, ~ 
14 mK and 140 mK for Ce2PdIn8 and CeCoIns, respectively, which 
indicates that the non-locality effect should be important at much 
lower temperatures than our T^^p- Thus, we can exclude these ef- 
fects as primary origins of the deviation from the T-linear behavior 

Reserved for Publication Footnotes 

6 I 

Footline Author 

We also note that the T^/"^ dependence of low-temperature super- 
fluid density has been derived from the Bardeen-Cooper-Schrieffer 
(BCS) to Bose-Einstein condensation (BEC) crossover below Tc 
|S16|. However, this assumes an s-wave state, which is not appli- 
cable to the nodal superconductors focused in the present study. 

SI3-Comparisons between iron-pnicteides and cuprates 

We comment on the difference in the high- Tc cuprate superconduc- 
tors. In clean single crystals of optimally doped YBa2Cu307_4, the 
superfluid density exhibits a rather wide range of T-linear behav- 
ior and the crossover temperature is as low as T*^^/Tc ~ 0.01. In 
hole-doped cuprates, the superconducting dome is separated from 
the antiferromagnetic order, and a possible QCP of the pseudogap 
phase has been discussed near the optimum doping of the supercon- 
ducting dome. The absolute value of the superfluid density shows a 
broad maximum near the optimal doping IIS17I . which is opposite to 
the iron-pnictide BaFe2(Asi_a;Pi:)2 case, where the sharp minimum 
is observed at the antiferromagnetic QCP IS181 . The latter case is 
consistent with the enhanced mass (which enters inversely to the su- 
perfluid density) associated with the QCP. This difference in doping 
dependence seems consistent with the difference in the temperature 
dependence of superfluid density; the y^/^ dependence is a common 
property associated with the mass enhancement due to the antiferro- 
magnetic QCP. 

To see the doping dependence of superfluid density in the iron- 
pnictide BaFe2(Asi_2,Pa;)2, we plot the low-temperature ps as a 
function of {T/TaY'^ in a wide doping range (Fig. S2). The anti- 
ferromagnetic QCP composition in this system has been located at 
X = 0.30 ISI8I . In high composition samples (x > 0.38) which are 
relatively far from the QCP, the low-temperature deviations from the 
T^^^ dependence is more pronounced than those in samples closer 
to the QCP (0.27 < x < 0.33). These results further support that 
the non-integer power law is related to the closeness to the magnetic 
QCP. In hole-doped cuprate superconductors, the antiferromagnetic 
QCP is actually located outside the low-doping end point of the su- 
perconducting dome. It deserves further studies to understand how 
the phase diagrams relate to these differences of low-energy quasi- 
particle excitations in cuprates and iron-based high- Tc superconduc- 

0.05 0.10 


Fig. S2. The normalized superfluid density ps as a function of (T/Tc)^-^ at 
low temperatures for BaFe2(Asi_£Pi:)2 single crystals. The original penetration 
depth data are taken from Ret. IS181 . 

SI4-Temperature dependence of the renormalization 

Near the QCP, the temperature dependence of the renormalization 
is an important factor. In the normal state above the upper critical 
field, the specific heat is largely enhanced at low temperatures |S19 1, 
indicating the strong temperature dependence of m* which may be 
approximated by 1/Vt at low temperatures (Fig. S3A). (Note that 
due to the field-induced QCP, this dependence is the strongest near 
the upper critical field and the actual temperature dependence at zero 
field may be weaker.) In the superconducting state, owing to the mo- 
mentum dependence of the renormalization we consider, l/vp{T) 
will level off at momentum dependent temperature. Near the nodes, 
the temperature dependence of renormalization will continue to low 
temperatures, whereas in the antinodal directions it will become in- 
dependent of temperature near Tc. 

If we assume v*p((j),T)/vF = m&x{C^T/Tc, \ cos(2<j!>)j^''^), 
we can calculate the temperature dependence of averaged mass en- 
hancement in the superconducting state (Fig. S3A). Here we separate 

A 1.15 



c^"^ 1.05 


• m*(7)/m*(0.2TJ 











- ' 1 ' 1 ' W 

— V((>)/''f = 


- - - 1/a/ 


- - 

= max(C(7■/7J"^|cos(2^)|'") / 



r'' , 1 , 1 

, 1 , \ 


Fig. S3. (A): Temperature dependence of the effective mass normalized by its 0.2Tc value In CeColns (black circles) determined from the specific heat 
measured at 5 T in the normal state (above the upper critical field) S19 can be fitted to a l/VT dependence (black line). The inverse of the effective 
Fermi velocity averaged over the Fermi surface In the superconducting state (corresponding to the first term in Eq. [1] In the main text) calculated assuming 
Vp{4>, T)/vp = max(Ci/T/Tc, | cos(2</))|-'/^) with a proper constant C = Vb\ cos 2(45° it 3°)|-'-''^, has a much weaker temperature dependence (red 
line). (B): Temperature dependence of the normalized superfluid density with (solid lines) and without (dashed lines) assuming the temperature dependence of 
renormalization. The total superfluid density (black) can be separated to the diamagnetic (red) and paramagnetic (blue) components. 

Footline Author 

PNAS I issue Date | Volume | Issue Number | 7 

the two terms in Eq. [1] in the main text as A = ^ — Ap ^, 
where the first term A^^ = |^|^ / "i^^f '^^ corresponds to the dia- 
magnetic current and the second term comes from the paramagnetic 
current carried by the thermally excited quasiparticles. In Fig. S3S 
the temperature dependence of each contribution is presented with 
and without the temperature dependence of renormalization. These 
results indicate that in the superconducting state the temperature de- 
pendence of the averaged renormalization is weak, which seems to be 
consistent with the recent high-frequency conductivity results IS20I . 
Moreover, the temperature dependence of the total superfluid density 
A^^ (T) remains almost unchanged even we consider the temperature 
dependence of the renormalization in the superconducting state. 

SI5-Effect of quantum critical exponent 

To see how the power-law of superfluid density varies with the quan- 
tum critical exponent /3, we have computed pa (T) with several values 
of /3 up to ~ 2. The low-temperature 1 - Ps{T) is fitted to the T" 
dependence up to T = Q.2Tc and the obtained power a is plotted 
against the exponent /3 in Fig. S4. We find that the power a varies 
linearly with /3, which can be approximated by a ~ 1 + /3/2. 

The fact that we universally obtained a close to 3/2 in the present 
quantum critical superconductors suggests that the exponent /? in 
these materials is not far from unity. Indeed, in /3-YbAlB4 in which 
zero-field quantum criticality has been suggested IIS2II . thermody- 
namic measurements shows /? ~ 1, which leads to a « 1.54 
(Fig. S4). In the study of YbRh2Si2, where the quantum critical field 
is 0.66 T, /? ~ 1 has also been reported |S22|. We note that in a sim- 
ilar system YbRh2(Sio.95Geo.o5)2 where the critical field is reduced 
to 0.027 T, smaller value of /3 ~ 2/3 has been suggested |S23|, but 
this value gives a ~ 1.36, which is in practise difficult to distinguish 
from the 3/2 power-law dependence. We also note that in CeCoIns 

and Ce2PdIn8 superconductors the estimated values for (3 from trans- 
port coefficient A near the upper critical field Hc2 are rather scattered 
(~ 1.37 for CeCoIns tS24i| and ~ 0.57 for Ce2PdIn8 | S14|), but the 
first-order transition of the low-temperature Hc2 due to the strong 
Pauli paramagnetism P S25I IS 121 may provide some complexity in 
evaluating the precise value of critical exponent. 

0.5 1.0 1.5 2.0 

Fig. S4. Low-temperature power a of the superfluid density obtained from ttie 
fitting up to T/Tc = 0.2 for the calculated ps{T/Tc) with different values of 
the quantum critical exponent /3. The line is a fit to a linear dependence. 

51. Shishido H, etal. (2002) Fermi surface, magnetic and superconducting properties 
of LaRhlns and CeTlnj (T: Co, Rh and Ir). J P/iys Soe Jpn 71:162-173. 

52. Kaczorowski D, Pikul AP, Gnida D, & Tran VH (2009) Emergence of a superconduct- 
ing state from an antiferromagnetic phase in single crystals of tlie tieavy fermion 
compound C e2Pdln8. Phys Rev LefM 03:027003; ibid. Phys Rev Lett 104:059702 

53. Carrington A, Giannetta RW, Kim JT, & Giapintzakis J (1999) Absence of nonlinear 
Meissner effect in YBa2Cu306 95. Ptiys /?evS59:R14173-R14176. 

54. Hashimoto K, et al. (2010) Line nodes in the energy gap of superconducting 
BaFe2(ASi_^P:c)2 single crystals as seen via penetration depth and thermal con- 
ductivity Pljys Rev B 81 :220501 (R). 

55. Prozorov R, Giannetta RW, Carrington A, & Araujo-Moreira FM (2000) Meissner- 
London state in superconductors of rectangular cross section in a perpendicular 
magnetic field. Phys Rev 6 62:1 1 5-1 1 8. 

56. Kini M, et al. (1990) A new ambient-pressure organic superconductor, k- 
(ET)2Cu[N(CN)2]Br, with the highest transition temperature yet observed (inductive 
onset T„ = 11.6K, resistive onset = 12.5 K). Inorg Chem 29:2555. 

57. Taylor OJ, Carrington A, & Schlueter JA (2008) Superconductor-insulator phase 
separation induced by rapid cooling of K-(BEDT-TrF)2Cu[N{CN)2]Br. Pliys Rev B 

58. Hirschfeld PJ & Goldenfeld N (1993) Effect of strong scattering on the low- 
temperature penetration depth of a cf-wave superconductor. Pliys Rev 6 48:421 9- 

59. Kosztin I & Leggett AJ (1997) Nonlocal effects on the magnetic penetration depth 
in d-wave superconductors. Phys Rev Le» 79:135. 

510. Xu D, Yip SK, & Sauls JA (1995) Nonlinear Meissner effect In unconventional 
superconductors. Phys Be/S 51:1 6233-1 6253. 

511. Abrikosov AA & Gor~fkov LP (1960) Contribution to the theory of superconduct- 
ing alloys with paramagnetic impurities Zh E/csp Teor F/z 39:1 781. [Sov. Phys. JETP 
/.ett.12, 1243(1961)]. 

SI 2. Toklwa Y, Gegenwart P, Gnida D, & Kaczorowski D (2011) Quantum criticality near 
the upper critical field of Ce2Pdln8. Pliys RevB84:140507(R). 

513. Movshovich R, et ai. (2001) Unconventional superconductivity in Celrlns and 
CeColhs : Specific heat and thermal conductivity studies. Piiys Rev Lett 86:51 52- 

514. Dong JK, et ai. (2011) Field-Induced quantum critical point and nodal supercon- 
ductivity in the heavy-fermion superconductor Ce2Pdln8. Phys Rev X 1 :01 1 01 0. 

SI 5. Settal R, et al. (2001) Quasi-two-dimensional Fermi surfaces and the de Haas- 
van Alphen oscillation in both the normal and superconducting mixed states of 
CeColng. J Phys Condens /Wa(fer13:L627-L634. 

SI 6. Chen Q, Kosztin I, & Levin K (2000) Unusual thermodynamical and transport sig- 
natures of the BCS to Bose-Einstein crossover scenario below t,.. Phys Rev Lett 

517. Tallon JL, Loram JW, Cooper JR, Panagopoulos C, & Bernhard C (2003) Super- 
fluid density in cuprate high-T,. superconductors: A new paradigm. Phys Rev B 

518. Hashimoto K, et ai. (201 2) A sharp peak of the zero-temperature penetration depth 
at optimal composition in the Iron-based superconductor BaFe2(ASi_ j.P^)2. Sci- 
ence 334: 1554-1 557. 

519. Bianchi A, Movshovich R, Vekhter I, Paglluso PG, & Sarrao JL (2003) Avoided 
antiferromagnetic order and quantum critical point In CeColng. Phys Rev Lett 

520. Trunclk CJS et ai. (2012) Nodal quasiparticle dynamics In the heavy 
fermion superconductor CeColng revealed by precision microwave spectroscopy. 
arXiv:1 21 0.5571 

521 . Matsumoto Y, et ai. (201 1 ) Quantum criticality without tuning in the mixed valence 
compound 8-YbAIB4. Science 331:31 6-31 9. 

522. Gegenwart P, et al. (2002) Magnetic-field induced quantum critical point in 
YbRh2Sl2. PhysRev ieff 89:056402. 

523. Custers J, ef ai. (2003) The break-up of heavy electrons at a quantum critical point. 
Wafure 424:524-527. 

524. Paglione J, ef ai. (2003) Field-induced quantum critical point in CeColng . Phys Rev 
Le(f 91:246405. 

525. Izawa K, ef a/. (2001) Angular position of nodes In the superconducting gap of 
quasi-2D heavy-fermion superconductor CeColng . Phys Rev Lett 87:057002. 

8 I Footline Author