# Full-Gap Superconductivity Robust against Disorder in Heavy-Fermion CeCuSi

###### Abstract

A key aspect of unconventional pairing by the antiferromagnetic spin-fluctuation mechanism is that the superconducting energy gap must have opposite sign on different parts of the Fermi surface. Recent observations of non-nodal gap structure in the heavy-fermion superconductor CeCuSi were then very surprising, given that this material has long been considered a prototypical example of a superconductor where the Cooper pairing is magnetically mediated. Here we present a study of the effect of controlled point defects, introduced by electron irradiation, on the temperature-dependent magnetic penetration depth in CeCuSi. We find that the fully-gapped state is robust against disorder, demonstrating that low-energy bound states, expected for sign-changing gap structures, are not induced by nonmagnetic impurities. This provides bulk evidence for -wave superconductivity without sign reversal.

Theories of unconventional superconductivity, where the pairing is not mediated by phonons, have been developed over the past decades to explain superconductivity in strongly correlated materials such as heavy fermions and high- cuprates. These theories have been challenged by recent and surprising results on the heavy-fermion superconductor CeCuSi Yamashita () which was the first discovered heavy fermion superconductor Steglich (), and as such, the first candidate for an unconventional superconducting state. The fact that superconductivity in CeCuSi emerges near a quantum critical point of antiferromagnetic order has led to the almost universally held conclusion that its superconductivity is unconventional with Cooper pairing mediated by spin fluctuations. The essence of this mechanism is that the momentum () dependent repulsive interactions can effectively pair the electrons as long as the superconducting gap changes sign in -space. Depending on the structure of the Fermi surface and spin-fluctuations, this can lead to different sign-changing gap structures such as -wave in cuprates or -wave in iron pnictides. For CeCuSi, early experiments such as nuclear quadrupole resonance relaxation rate Ishida (); Fujiwara () and specific heat Arndt () suggested -wave superconductivity with line nodes in . Inelastic neutron scattering measurements have shown an enhancement of magnetic spectral weight at around , which has been interpreted as a spin resonance expected for a sign-changing also consistent with -wave symmetry Stockert (); Eremin ().

In contrast to this, recent experiments which have combined specific heat Kittaka (), penetration depth, and thermal conductivity measured down to very low temperatures have shown that gap nodes do not exist at any point on the Fermi surface of CeCuSi Yamashita (). This nodeless structure might still be explained by a spin-fluctuation mechanism if the points in -space where the gap changes sign do not coincide with the Fermi surface sheets, as is the case for most iron-pnictide materials. Specific calculations for CeCuSi have shown that spin-fluctuations can lead to an -type structure but the closeness and corrugations of the Fermi-surface sheets means that accidental nodes are unavoidable Ikeda (). However, a nodeless state cannot be ruled out by this alone because the experimental Fermi surface of CeCuSi has not been fully determined. Hence, experiments to specifically probe for the presence or absence of a nodeless sign changing gap structure are needed.

One such experiment is the effect of impurity scattering on . It was demonstrated that for CeCuSi, increasing impurity scattering leads only to a very weak reduction in Yamashita (); Adrian () which would appear to be inconsistent with a sign-changing . However, given the quantitative nature of this argument further experimental confirmation is needed. Previous attempts at phase-sensitive measurements focused on the Josephson effect between CeCuSi and Al Poppe (); Sumiyama (). Although a finite Josephson current and conventional Fraunhofer pattern were observed, the results are not conclusive because the polycrystalline nature of the samples used complicates the interpretation. Moreover, in such measurements the possibility of a surface-induced -wave component cannot be ruled out Sumiyama (). Thus bulk measurements that are sensitive to a possible sign change in of CeCuSi are desirable.

Here we report on systematic measurements of the temperature dependence of the magnetic penetration depth in single crystals of CeCuSi, where impurity scattering has been introduced in a controlled way by electron irradiation. The irradiation creates point-like defects that act as nonmagnetic scatterers, which in the case of sign-changing order parameters should induce Andreev bound states at low energies. In the case, therefore, the fully-gapped state is expected to change to a gapless state with low-lying quasiparticle excitations Wang (). Indeed, in some iron-based superconductors a change in low-temperature from an exponential dependence to a dependence with increasing defect level is found which supports the model of an gap structure Mizukami (). In CeCuSi, we observe essentially no change in the low-temperature behavior of over a wide range of impurity scattering rates, indicating the robustness of the fully-gapped state against disorder. This provides strong evidence that the gap structure of CeCuSi is non-sign changing -wave state ().

High-quality single crystals of CeCuSi were synthesized by the self flux method and characterized by x-ray diffraction Seiro (). The crystals were cut into samples with typical dimension about m (in the plane) and thickness about 50 m (along the axis). To introduce spatially homogeneous defects in a controllable way, we employed 2.5-MeV electron beam irradiation at the SIRIUS Pelletron linear accelerator operated by the Laboratorie des Solides Irradiés (LSI) at École Polytechnique. This incident energy is sufficient to form vacancy-interstitial (Frenkel) pairs, which act as point-like defects. The attenuation distance of these irradiation electrons in CeCuSi is about 2.7 mm, which is much longer than our sample thickness. During the irradiation, the sample was kept at K by using a liquid hydrogen bath, which is important to prevent defect migration and clustering. The temperature dependent changes in the in-plane magnetic penetration depth were measured using the tunnel diode oscillator (TDO) technique operating at 14 MHz carrington99 () in a dilution refrigerator down to mK. The shift of the resonant frequency of the oscillator including the sample directly reflects the change in the magnetic penetration depth, . The geometric factor is determined from the geometry of the sample and the coil prozorov00 ().

Figure 1(a) shows the temperature dependence of the normalized frequency shifts for the pristine and irradiated samples with irradiation doses of 2.2, 3.7, and 4.8 C/cm. The superconducting transition remains sharp after irradiation, indicating that the point-like defects are introduced uniformly. In Fig. 1(b) we show the dose dependence of determined by the onset of the diamagnetic signal and compare this to the evolution of the and residual resistivity reported previously Yamashita (). Here the vertical axis is normalized by K, which is estimated from the linear extrapolation to the zero defect limit. These two independent results measured in different crystals are in good agreement, demonstrating that all irradiated samples are homogeneous.

Figure 2 depicts the variations of penetration depth as a function of for pristine and irradiated crystals, after normalization by their values of . We find no significant change in the temperature dependence of at low temperatures, and all the curves almost collapse to a single one. This indicates that the introduced defects have essentially no effect on the low-energy quasiparticle excitations. The fact that we do not observe any evidence for a Curie-like upturn in down to the lowest temperature of mK even in most irradiated sample implies that introduced point-like defects are nonmagnetic in nature. Any magnetic impurities would result in a Curie-like upturn in the normal-state susceptibility, which would lead to an additional contribution to , with . and are respectively the density and the effective moment of the magnetic impurity Cooper (); Malone (); Mizukami (). We estimate about 2 vacancies per 1000 Ce atoms are formed per 1 C/cm electron irradiation, so the absence of the upturn in the 4.8 C/cm sample gives an upper limit of about per defect. This is much smaller than the moment of for a free Ce ion with total momentum .

Having established the nonmagnetic nature of the defects, we analyze the low-temperature data to make a more quantitative analysis of the changes in the gap structure. We use two procedures: one is a fit to the power-law dependence with a variable exponent and the other is a fit to the exponential dependence with a variable minimum gap . In both cases, we examine the changes in the fitting parameters as a function of the upper bound of fitted temperature range, . The obtained results for and are shown in Figs. 3(a) and (b), respectively. For the power-law fitting procedure we find that all the data before and after irradiation give high exponent values , which is far beyond the dirty-limit line-node case of exponent. This high-power dependence is practically indistinguishable from an exponential temperature dependence. The obtained gap values in the exponential fits are smaller than the BCS value of , suggesting that the gap has strong dependence with a large difference between minimum and maximum values.

Our principal finding is the robustness of the fully-gapped superconductivity against disorder in CeCuSi. This is most clearly demonstrated by plotting the exponent in the power-law fit as a function of pair-breaking parameter in Fig. 3(c), in which we compare with the typical results for -wave CeLaCoIn Kim () and for -wave BaFe(AsP) Mizukami (). The impurity scattering time is calculated with . In La-substituted CeCoIn the exponent increases with impurity scattering and saturates at , which is consistent with the gapless state expected theoretically in the dirty -wave superconductors Hirschfeld () and established experimentally for Zn-substituted YBaCuO Bonn (). In optimal BaFe(AsP), which has a -linear behavior due to the accidental line nodes in the clean limit Hashimoto (), the exponent initially shows a large increase from to , indicating the lifting of nodes by the impurity-induced averaging effect of the dependence, which occurs only when the nodes are not symmetry protected. Further irradiation yields a decrease of toward the gapless value of 2, demonstrating the creation of the low-energy states that are expected only for sign-changing cases. These results established a nodal -wave state in this iron pnictide. Thus in both -wave and -wave cases, a gapless state with the exponent appears for pair-breaking parameter of the order of unity. In stark contrast, our data for CeCuSi reveal that the exponent remains high () even when exceeds 20, which clearly indicates the absence of impurity-induced low-energy states, evidencing no sign change in .

The minimum superconducting gap size normalized by shows no appreciable change against the relative suppression of the transition temperature as shown in Fig. 3(d). At first glance this appears counterintuitive because the gap averaging effect due to impurity scattering might be expected to lead to an increase of the minimum . However similar behavior is also observed in the protypical two gap superconductor MgB, where following an initial increase, remains unchanged in a wide suppression range up to about half of Putti (); Daghero ().

The anisotropic nature of the gap structure in CeCuSi can be seen in the full temperature dependence of the normalized superfluid density . To calculate , we need the value of for each sample, which cannot be determined directly by using the TDO technique. So instead we have estimated from the lower critical field measured by micro Hall-probe magnetometry as described in Ref. SM (). To minimize errors due to geometrical demagnetization factors we measured the same sample both before and after irradiation with a dose of 1.9 C/cm, which reduced from 0.64 K to 0.52 K. For the irradiated sample we found mT at 100 mK compared to mT in the pristine sample Yamashita (). From this we estimate that is increased from nm for the unirradiated sample to nm for the irradiated one. An increase in upon irradiation is expected because the effective penetration depth depends on the mean free path of quasiparticles. For the unirradiated sample we estimated the in-plane mean free paths and coherence lengths to be nm and nm respectively so the sample is between the clean and dirty limits. Then from the change in we would expect to decrease by a factor 2 for this irradiation level, pushing the sample closer to the dirty limit and thus increasing . In the dirty limit, is expected to be proportional to , which appears to hold as shown in the inset of Fig. 4. From this relation we estimate nm for the most irradiated sample (4.8 C/cm). Figure 4 displays the extracted curves before and after irradiation, which again show the robustness of flat temperature dependence at low temperatures indicating the absence of the low-energy states. The multi-gap or strong -dependent nature of manifests itself in the concave curvature near . Contrary to the case of MgB Fletcher (), however, a simple two-gap model does not fit the data very well. Possible reasons for this include significant interband scattering and largely varying for each band. For the irradiated sample, the concave curvature of near is less pronounced and the curve becomes closer to the single-gap -wave one, which is consistent with the reduced anisotropy of by impurity scattering.

In summary, systematic measurements of magnetic penetration depth in electron-irradiated single crystals of CeCuSi show that nonmagnetic impurity scattering does not induce any low-energy quasiparticle excitations. This provides bulk evidence for the absence of a sign change in the gap function in the superconducting state of this heavy-fermion superconductor. The -wave state inferred in this study is generally a manifestation of on-site attractive interactions, but how this can overcome the strong Coulomb repulsion in such a strongly correlated electron system calls for new theoretical approaches beyond the wide-spread spin-fluctuation based unconventional mechanism of superconductivity. Very recent calculations show that in the vicinity of magnetic quantum critical point, the orbital fluctuations may lead to -wave superconductivity Tazai (). Indeed, the importance of orbital degrees of freedom has been pointed out in several aspects for some Ce-based materials including CeSi, where is a transition metal element Hattori (); Ren (); Dong () . Thus the relationship between the orbital effects and gap symmetry in heavy-fermion superconductors deserves further studies.

We thank H. Ikeda and H. Kontani for fruitful discussion. We also thank the SIRIUS team, O. Cavani, B. Boizot, V. Metayer, and J. Losco, for running electron irradiation at Laboratorie des Solides Irradiés (LSI) in École Polytechnique [supported by EMIR network, proposal Nos. 16-0398, 16-9513, and 17-1353]. This work was supported by Grants-in-Aid for Scientific Research(KAKENHI) (No. 25220710, No. 15H02106) and Grants-in-Aid for Scientific Research on Innovative Areas “Topological Materials Science” (No. 15H05852) from Japan Society for the Promotion of Science (JSPS), and by the UK Engineering and Physical Sciences Research Council (grant no.EP/L025736/1 and EP/L015544/1).

## References

- (1) T. Yamashita, T. Takenaka, Y. Tokiwa, J. A. Wilcox, Y. Mizukami, D. Terazawa, Y. Kasahara, S. Kittaka, T. Sakakibara, M. Kończykowski, S. Seiro, H. S. Jeevan, C. Geibel, C. Putzke, T. Onishi, H. Ikeda, A. Carrington, T. Shibauchi, and Y. Matsuda, Sci. Adv. 3, e1601667 (2017).
- (2) F. Steglich, J. Aarts, C. D. Bredl, W. Lieke, D. Meschede, W. Franz, and H. Schäfer, \bibfieldjournal Phys. Rev. Lett. 43, 1892 (1979).
- (3) K. Ishida, Y. Kawasaki, K. Tabuchi, K. Kashima, Y. Kitaoka, K. Asayama,C. Geibel, and F. Steglich, \bibfieldjournal Phys. Rev. Lett. 82, 5353 (1999).
- (4) K. Fujiwara, Y. Hata, K. Kobayashi, K. Miyoshi, J. Takeuchi, Y. Shimaoka, H. Kotegawa, T. C. Kobayashi, C. Geibel, and F. Steglich, \bibfieldjournal J. Phys. Soc. Jpn. 77, 123711 (2008).
- (5) J. Arndt, O. Stockert, K. Schmalzl, E. Faulhaber, H. S. Jeevan, C. Geibel, W. Schmidt, M. Loewenhaupt, and F. Steglich, \bibfieldjournal Phys. Rev. Lett. 106, 246401 (2011).
- (6) O. Stockert, J. Arndt, E. Faulhaber, C. Geibel, H. Jeevan, S. Kirchner, M. Loewenhaupt, K. Schmalzl, W. Schmidt, Q. Si, and F. Steglich, \bibfieldjournal Nat. Phys. 7, 119 (2011).
- (7) I. Eremin, G. Zwicknagl, P. Thalmeier, and P. Fulde, \bibfieldjournal Phys. Rev. Lett. 101, 187001 (2008).
- (8) S. Kittaka, Y. Aoki, Y. Shimura, T. Sakakibara, S. Seiro, C. Geibel, F. Steglich, H. Ikeda, and K. Machida, \bibfieldjournal Phys. Rev. Lett. 112, 067002 (2014).
- (9) H. Ikeda, M.-T. Suzuki, and R. Arita, \bibfieldjournal Phys. Rev. Lett. 114, 147003 (2015).
- (10) G. Adrian and H. Adrian, \bibfieldjournal Europhys. Lett. 3, 819 (1987).
- (11) U. Poppe, \bibfieldjournal J. Magn. Magn. Mater. 52, 157 (1985).
- (12) A. Sumiyama, N. Miyakawa, Y. Ushida, G. Motoyama, A. Yamaguchi, and Y. Oda, \bibfieldjournal J. Phys: Conf. Ser. 273, 012086 (2011).
- (13) Y. Wang, A. Kreisel, P. J. Hirschfeld, and V. Mishra, \bibfieldjournal Phys. Rev. B 87, 094504 (2013).
- (14) Y. Mizukami, M. Kończykowski, Y. Kawamoto, S. Kurata, S. Kasahara, K. Hashimoto, V. Mishra, A. Kreisel, Y. Wang, P. J. Hirschfeld, Y. Matsuda, and T. Shibauchi, \bibfieldjournal Nat. Commun. 5, 5657 (2014).
- (15) S. Seiro, M. Deppe, H. Jeevan, U. Burkhardt, and C. Geibel, \bibfieldjournal Phys. Status Solidi B 247, 614 (2010).
- (16) A. Carrington, R. W. Giannetta, J. T. Kim, and J. Giapintzakis \bibfieldjournal Phys. Rev. B 59, R14173 (1999).
- (17) R. Prozorov, R. W. Giannetta, A. Carrington, and F. M. Araujo-Moreira \bibfieldjournal Phys. Rev. B 62, 115 (2000).
- (18) J. R. Cooper, \bibfieldjournal Phys. Rev. B 54, R3753 (1996).
- (19) L. Malone, J. D. Fletcher, A. Serafin, A. Carrington, N. D. Zhigadlo, Z. Bukowski, S. Katrych, and J. Karpinski, \bibfieldjournal Phys. Rev. B 79, 140501 (2009).
- (20) H. Kim, M. A. Tanatar, R. Flint, C. Petrovic, R. Hu, B. D. White, I. K. Lum, M. B. Maple, and R. Prozorov, \bibfieldjournal Phys. Rev. Lett. 114, 027003 (2015).
- (21) E. E. M. Chia, D. J. Van Harlingen, M. B. Salamon, B. D. Yanoff, I. Bonalde, and J. L. Sarrao, \bibfieldjournal Phys. Rev. B 67, 014527 (2003).
- (22) C. Petrovic, S. L. Bud’ko, V. G. Kogan, and P. C. Canfield, \bibfieldjournal Phys. Rev. B 66, 054534 (2002).
- (23) M. Putti, M. Affronte, C. Ferdeghini, P. Manfrinetti, C. Tarantini, and E. Lehmann, \bibfieldjournal Phys. Rev. Lett. 96, 077003 (2006).
- (24) D. Daghero, A. Calzolari, G. A. Ummarino, M. Tortello, R. S. Gonnelli, V. A. Stepanov, C. Tarantini, P. Manfrinetti, and E. Lehmann, \bibfieldjournal Phys. Rev. B 74, 174519 (2006).
- (25) P. J. Hirschfeld and N. Goldenfeld, \bibfieldjournal Phys. Rev. B 48, 4219 (1993).
- (26) D. A. Bonn, S. Kamal, K. Zhang, R. Liang, D. J. Baar, E. Klein, and W. N. Hardy, \bibfieldjournal Phys. Rev. B 50, 4051 (1994).
- (27) K. Hashimoto, M. Yamashita, S. Kasahara, Y. Senshu, N. Nakata, S. Tonegawa, K. Ikada, A. Serafin, A. Carrington, T. Terashima, H. Ikeda, T. Shibauchi, and Y. Matsuda, Phys. Rev. B 81, 220501(R) (2010).
- (28) See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.*** for details of the measurements, which includes Refs. Brandt (); Putzke ().
- (29) E. H. Brandt, Phys. Rev. B 60, 11939 (1999).
- (30) C. Putzke, P. Walmsley, J.D. Fletcher, L. Malone, D. Vignolles, C. Proust, S. Badoux, P. See, H. E. Beere, D. A. Ritchie, S. Kasahara, Y. Mizukami, T. Shibauchi, Y. Matsuda, and A. Carrington, Nat. Commun. 5, 5679 (2014).
- (31) J. D. Fletcher, A. Carrington, O. J. Taylor, S. M. Kazakov, and J. Karpinski, Phys. Rev. Lett. 95, 097005 (2005).
- (32) R. Tazai, Y. Yamakawa, M. Tsuchiizu, and H. Kontani, J. Phys. Soc. Jpn. 86, 073703 (2017).
- (33) K. Hattori, J. Phys. Soc. Jpn. 79, 114717 (2010).
- (34) Z. Ren, L. V. Pourovskii, G. Giriat, G. Lapertot, A. Georges, and D. Jaccard, Phys. Rev. X 4, 031055 (2014).
- (35) R. Dong, X. Wan, X. Dai, and S. Y. Savrasov, Phys. Rev. B 89, 165122 (2014).

## Appendix A Supplemental Material

### a.1 Lower Critical Field

In order to estimate the absolute value of the penetration depth in the pristine and irradiated samples, we employ a technique using a micro-Hall array to determine the lower critical field of the samples. The sample is placed on top of an array of Hall sensors, where each Hall sensor measures the perpendicular component of the magnetic induction through the active area as a function of an applied field . The sample is cooled in zero-field, with the vacuum can of the refrigerator shielded in mu-metal to reduce the effect of the earth’s field. A small coil inside the can provides the magnetic field.

At a given temperature, the applied field, oriented perpendicular to the Hall sensors and parallel to the -axis of the sample, is increased from zero to the maximum value (positive or negative). At low field, there is an increase in with due to the incomplete shielding of the Hall sensors by the sample. At a well defined field, the magnetic induction increases sharply, indicating the field at which flux has entered the sample. An example sweep is shown in Fig. S1(a). A linear background is subtracted from the data to account for the incomplete shielding below , and the values of for each sweep direction are indicated. After each sweep, the sample is warmed above and cooled again in zero-field. The overall temperature dependence is given in Fig. S1(b).

The field at which flux enters the sample is not equal to the lower critical field , but it is the lower critical field reduced by a factor related to the demagnetising effects due to the geometry of the sample. The lower critical field is related to the penetration field by the relation given by Brandt for a strip Brandt2 ():

(S1) |

where we take as the shorter of the two in-plane dimensions and is the thickness in the -axis direction. Although this expression was calculated for a strip of infinite third dimension, it has been found to provide a good description of samples even if the two in-plane dimensions are approximately equal (see for example Ref. Putzke2 () where consistent results were found when the aspect ratio of a sample was varied by cutting). For the sample with an irradiation dose of 1.9 C cm (dimensions mm) we find a value of mT at mK. The penetration depth is then determined from the Ginzburg-Landau equation

(S2) |

where is the flux quantum and is the coherence length. Using nm found for the pristine sample Yamashita2 (), gives a value of nm at mK which within the uncertainty we take as being equal to the zero temperature value.

## References

- (S1) E. H. Brandt, Phys. Rev. B 60, 11939 (1999).
- (S2) C. Putzke, P. Walmsley, J.D. Fletcher, L. Malone, D. Vignolles, C. Proust, S. Badoux, P. See, H. E. Beere, D. A. Ritchie, S. Kasahara, Y. Mizukami, T. Shibauchi, Y. Matsuda, and A. Carrington, Nat. Commun. 5, 5679 (2014).
- (S3) T. Yamashita, T. Takenaka, Y. Tokiwa, J. A. Wilcox, Y. Mizukami, D. Terazawa, Y. Kasahara, S. Kittaka, T. Sakakibara, M. Kończykowski, S. Seiro, H. S. Jeevan, C. Geibel, C. Putzke, T. Onishi, H. Ikeda, A. Carrington, T. Shibauchi, and Y. Matsuda, Sci. Adv. 3, e1601667 (2017).