Introduction

Unconventional superconductivity (SC) often is observed in close proximity to a magnetically ordered phase, where the SC transition temperature Tc forms a dome against a non-thermal control parameter, such as the external pressure, chemical substitution, or magnetic field1,2,3,4,5,6. At an optimal value of the tuning parameter, where Tc is the highest, normal state properties do not follow predictions for Landau–Fermi liquids: the electrical resistivity (ρ) does not exhibit a T2 dependence, and the electronic specific heat coefficient (γ=C/T) does not saturate, but rather diverges with decreasing temperature1, 2, 7. These non-Fermi liquid (NFL) behaviors arise from incoherent critical fluctuations associated with a quantum critical point (QCP) hidden inside the SC dome of heavy fermion compounds and some Fe-based superconductors, such as BaFe2(As1−xPx)21, 2, 4, 6, 8. Because the zero-temperature quantum phase transition is typically not accessible without destroying SC, the role of critical magnetic fluctuations on properties of unconventional superconductors has yet to be explored in depth.

The critical current (Ic), which limits the current capacity of a zero-resistance state, is characteristically taken to depend on the strength of vortex pinning, which, in turn, is determined by the geometry and distribution of microstructural defects9,10,11. Because application of pressure should not lead to the creation of different or additional defects or to a substantial change in sample dimensions, Ic in relation to Tc should be at most weakly pressure-dependent. A substantial variation in Ic(P) or Ic/Tc(P), then, logically, should be attributed to intrinsic changes in the superconducting state itself. For example, the zero-field critical current density Jc (equal to Ic/A, where A is the sample cross sectional area perpendicular to current) of the hole-doped high-Tc cuprate superconductor Y0.8Ca0.2Ba2Cu3Oy has a sharp peak that is centered on a critical hole-doping where the pseudogap boundary line projects to zero temperature, and that is attributed in model calculations to changes in the superfluid density12, 13. These results indicate that Ic measurements may provide an opportunity to explore the relationship between unconventional SC and any QCP that is hidden beneath the SC dome.

Here we report a peak in the zero-field critical current, Ic(0), at a critical pressure Pc in pure CeRhIn5 (Rh115) and 4.4% Sn-doped CeRhIn5 (SnRh115), where their respective antiferromagnetic boundary TN(P) extrapolates to T = 0 K inside a dome of pressure-induced SC. The temperature dependence of Ic(0)s for pure Rh115 and SnRh115 under pressure is similar to that of superconducting CeCoIn5, which is close to quantum criticality at ambient pressure. Normalized values of Ic(T, P) follow a common universal curve for each material, suggesting an intrinsic, fundamental connection to quantum criticality. Supporting this conclusion, the magnetic field dependence of the flux-pinning force (Fp = Ic × μ0H), normalized to its maximum value, also forms a pressure-invariant universal curve for each compound. As will be discussed, these discoveries demonstrate that the pressure evolution of zero-field Ic is determined mainly by quantum critical fluctuations, and that the peak in Ic is a direct link to the hidden QCP.

Results

Temperature–pressure phase diagrams

Figure 1a and b presents a contour plot of the zero-field Ic(P, T) in the SC phase and the in-plane resistivity ρab(P, T) in the normal state for pure CeRhIn5 (Rh115) and 4.4% Sn-doped CeRhSn0.22In4.78 (SnRh115) single crystals. The dependence on pressure of the in-plane resistivity and current–voltage curves upon which Fig. 1 is based is displayed in Supplementary Figs. 1 and 2, respectively. The quantum critical region veiled by the superconducting phase is fully exposed by the pressure dependence of zero-field Ic(T). A sharp peak in the value of Ic(T) is clearly observed for pressures around the QCP at Pc, where a large enhancement in the resistivity is accompanied by strong quantum fluctuations3,4,14. In addition, Ic(P) abruptly increases at pressures around Pc*, the critical pressure where coexisting phases of magnetism and SC evolve into a single SC state. In undoped Rh115, large differences between Tcs measured by heat capacity (C) and resistivity (ρ) at pressures below Pc* are ascribed to textured SC originating from an incommensurate long-range magnetic order15,16,17,18.

Fig. 1
figure 1

Temperature–pressure phase diagrams of CeRhIn5 and CeRhSn0.22In4.78 single crystals. In the superconducting state below Tc(P), false colors denote the magnitude of the zero-field critical current Ic(P, T). At temperatures above Tc(P), false colors reflect the magnitude of the in-plane resistivity ρab(P, T). a CeRhIn5 (Rh115) and b CeRhSn0.22In4.78 (SnRh115). For both materials, ρab(P, T) is enhanced around the quantum critical point Pc due to pronounced incoherent inelastic scattering. Similarly, the zero-field Ic(P, T) is largest at Pc, where the QCP is expected, as indicated by the arrow. In both a and b the vertical hashed rectangle is at Pc*, the pressure that separates a phase of coexisting superconductivity and magnetism from a purely SC phase for P > Pc*. Open squares in both a and b represent the antiferromagnetic transition temperature (TN). SC transition temperature (Tc) of Rh115 is evaluated from specific heat (Tc,C) and resistivity (Tc,ρ) measurements, and Tc of SnRh115 is determined as Tc onset (Tc,on) and 50% (Tc,50%) of the normal state resistivity value at Tc,on. AFM, SC, and NFL stand for antiferromagnetic, superconducting, and non-Fermi liquid regions, respectively

Temperature dependences of the zero-field critical current

The antiferromagnetic transition temperature (TN ~ 3.8 K) in pure Rh115 is suppressed by Sn doping, which induces a shift of its extrapolated T = 0 K antiferromagnetic transition, and pressure-induced superconductivity emanates from the tuned QCP4,5 (see Supplementary Fig. 3). Figure 2a and b shows the temperature dependence of the zero-field critical current, Ic(0), for Rh115 and SnRh115 at several pressures, respectively. Here, Ic is determined by using the voltage criterion of 0.1 μV (see Supplementary Fig. 2). Analysis of the flux-pinning force, Fp = Ic × μ0H, shows that the normalized flux-pinning force follows a power-law dependence on magnetic field, fp(h) hp(1−h)q, and is peaked around hpeak ≈ 0.6, which is characteristic of type-II superconductors with weak pinning (see Supplementary Fig. 4). Here, the normalized pinning force is fp = Fp/Fp,max and the reduced field is h = H/Hirr, where Fp,max is the maximum flux-pinning force and Hirr is the irreversible field. The dependence on temperature of the critical current has been widely explained by Ic(t)/Ic(0) = (1−t2)α(1 + t2)β for type-II superconductors, such as high-Tc cuprates, Fe-based superconductors, and MgB219,20,21,22, where t = T/Tc is the reduced temperature. When Tc variations surrounding defects are important (δTc-pinning), α = 7/6 and β = 5/6, but α = 5/2 and β = −1/2 for δl-pinning that arises from spatial variations in the charge-carrier mean free path (l) near a lattice defect10, 19, 23. These functional forms are shown by the dotted and dashed lines for δTc-pinning to δl-pinning in Fig. 2c, respectively. A crossover of the mechanism from δTc-pinning to δl-pinning has been often reported by introducing additional defects via chemical substitution or heavy ion irradiation, indicating that δTc-pinning is preferred in clean crystals19,20,21.

Fig. 2
figure 2

Temperature dependences of the zero-field critical current for Ce-based heavy fermion materials under pressure. a Temperature dependence of the zero-field critical current, Ic(0), for CeRhIn5 at various pressures. b Zero-field Ic for CeRhSn0.22In4.78 at various pressures. c Reduced temperature (t = T/Tc,0) dependence of Ic, Ic(t), normalized by its value at t = 0.4 for Rh115 and SnRh115 at representative pressures and for CeCoIn5 at ambient pressure. The normalized values of Ic(P, t) for all crystals can be described by a single curve, Ic(t) (1 − t2)5/6(1 + t2)2/3, indicating universal behavior of Ic(t) with respect to pressure in the CeMIn5 (M = Co, Rh) materials. Dotted and dashed curves are for δTc-pinning and δl-pinning, respectively, as discussed in the text

The values of Ic(t) for SnRh115 at pressures around Pc are fitted together with those for Rh115 and CeCoIn5 in Fig. 2c. The temperature dependence of Ic(0) for CeCoIn5 at ambient pressure is measured to compare it with that of CeRhIn5, because Rh115 is believed to have a SC pairing mechanism similar to that in CeCoIn5. The values of Ic(0) for all samples can be expressed well by one curve with the relation Ic(t) (1−t2)5/6(1 + t2)2/3, which is distinct from that for Ic(t) controlled by either δTc-pinning or δl-pinning. This universal curve underscores that the origin of the zero-field Ic is the same for each compound, and that it does not change under pressure for these Ce-based quantum critical materials. The fact that external pressure does not create new defects inside the crystals suggests that the pressure evolution of Ic should be related to the pressure dependence of the SC coupling strength.

Discussion

Figure 3a presents the pressure dependences of Ic(0) and Tc,0 for SnRh115, which are similar to each other. However, their relative fractional variations in Ic and Tc, γI ≡ Ic,0(P)/Ic,0(Pc) × 100 and γT ≡ Tc,0(P)/Tc,0(Pc) × 100, where Ic,0(Pc) is Ic extrapolated to zero temperature at Pc and Tc,0(Pc) is the SC transition temperature at Pc, are much different, as shown in Fig. 3b: at 1.0 GPa, the critical current is 45% of the maximum value, and Tc,0 is 67% of its maximal value. The stronger pressure dependence of Ic relative to that of Tc,0 is clearly visible in ratio Ic,0/Tc,0 for SnRh115, as presented in Fig. 3c, d. An abrupt enhancement in Ic,0/Tc,0 is observed at Pc*, and the peak in the pressure dependence of Ic,0/Tc,0 is achieved at Pc. The dc conductivity (=σdc) at Tc onset is shown as a function of pressure in the right ordinate of Fig. 3c, where a minimum value appears near Pc (see Supplementary Fig. 5). The anti-correlation between Ic,0/Tc,0 and σdc in these Ce-based quantum critical compounds may be related with the presence of the hidden QCP at Pc, because the associated critical quantum fluctuations not only act as the SC pairing glue, but also strongly enhance incoherent electron scattering, thus leading to a minimum in σdc at Pc24, 25. Homes’ scaling relation26,27,28 states that the superfluid density ns is proportional to σdcTc in many correlated superconductors and, consequently, that the ratio ns/Tc should be proportional to σdc. The fact that σdc is the minimum at Pc, where Ic,0/Tc,0 is the maximum in these Ce-based compounds, suggests a violation of Homes’ scaling if the strength of the condensate ns is proportional to the critical current Ic,0. Pressure-dependent optical conductivity and/or penetration depth experiments that directly measure ns will be important to provide a stringent test for the validity of Homes’ law in quantum critical superconductors.

Fig. 3
figure 3

Pressure evolution of the zero-field critical current in Sn-doped CeRhIn5. a Pressure dependences of Tc,0 and Ic,0 for SnRh115, where Ic,0 is the value of Ic obtained from an extrapolation of data in Fig.  2b zero Kelvin. b Fractional variations in Ic,0 and Tc,0 for SnRh115 under pressure. The fractions are defined as γI (at 0 K) ≡ Ic,0(P)/Ic,0(Pc) × 100 and γT ≡ Tc,0(P)/Tc,0(Pc) × 100, where Ic,0(Pc) is Ic extrapolated to zero temperature at Pc and Tc,0(Pc) is the superconducting transition temperature at Pc. Values of γI are plotted as a function of pressure for measured or estimated Ic(T, P) at 0 K (squares) and 0.3 K (circles). c The ratio between the critical current and SC transition temperature, Ic,0/Tc,0, plotted together with the dc conductivity at Tc onset, σdc, as a function of the pressure difference PPc, where Pc = 1.35 GPa is the QCP. d A contour plot of Ic,0/Tc,0 displayed in the temperature (T) and pressure (PPc) plane. The ratio Ic,0/Tc,0 forms a dome centered around the quantum critical point Pc and its values decrease with distance from Pc

Our study demonstrates that the critical current, a fundamental superconducting parameter, is a powerful tool for investigating the presence of a hidden QCP inside the superconducting dome without destroying the superconducting phase. The dependence on temperature of the zero-field Ic for both pure Rh115 and Sn-doped Rh115 exhibits the same functional form under pressure, underscoring that the peak at Pc in the pressure dependence of Ic arises from an enhanced fluctuations around the hidden QCP. Even though these results are specific to the Ce115 heavy-fermion materials, the prediction of similar results for the hole-doping dependence of the critical current density Jc(x) in high-Tc cuprates29 suggests a universal behavior of Jc among unconventional superconductors. These discoveries should stimulate more theoretical and experimental effort to understand the intimate link between quantum criticality and the origin of unconventional superconductivity in various families of correlated electronic systems.

Methods

Measurement outline

CeRhIn5, Sn-doped CeRhIn5, and CeCoIn5 single crystals were synthesized by the indium (In) self-flux method30,31,32. Pressure was generated in a hybrid clamp-type pressure cell with Daphne 7373 as the pressure-transmitting medium, and the pressure was determined by monitoring the shift in the value of Tc for lead (Pb). Measurements of current–voltage (IV) characteristics under pressure were performed in a Heliox VL system (Oxford Instruments) with a vector magnet (y = 5 T and z = 9 T, American Magnetics Inc.) and in a Physical Property Measurement System (PPMS 9 T, Quantum Design), where the current was provided by a Keithley 6221 unit and the voltage was measured with a Keithley 2182A nanovoltmeter.

Measurement details

Measurements of I–V characteristics were performed in a pulsed mode to minimize Joule heating developed at Ohmic contacts to the samples and copper (Cu) wires between the pressure cell and the connector. The duration of the pulsed current was 10–11 ms, and the repetition rate was one pulse every 2 s, which was sufficient to eliminate Joule heat in the samples33, 34. A standard four-probe method was used to determine IV, and good Ohmic contact to samples was achieved by using silver epoxy. The critical current was based on a 10−7 V criterion35, which was averaged over three measurements. The dimensions of the measured crystals were 920×330×20, 650×200×22, and 1100×200×47 μm3 for CeRhIn5 (Rh115), CeRhSn0.22In4.78 (SnRh115), and CeCoIn5, respectively. The magnetic-field dependence of the critical current was measured at several pressures and the flux-pinning force (Fp) was estimated from the relation Fp = Ic × μ0H 36,37,38

Data availability

The data sets generated and/or analyzed in this study are available from the corresponding author on reasonable request.