Unusually weak irradiation effects in anisotropic iron-based superconductor RbCa2Fe4As4F2

We report on the effects of 3.5 MeV proton irradiation in RbCa$_2$Fe$_4$As$_4$F$_2$, an iron-based superconductor with unusual properties in between those of the pnictides and of the cuprate high-temperature superconductors. We studied how structural disorder introduced by ion bombardment affects the critical temperature, superfluid density and gap values by combining a coplanar waveguide resonator technique, electric transport measurements and point-contact Andreev-reflection spectroscopy. We find an unusually weak dependence of the superconducting properties on the amount of disorder in this material when compared to other iron-based superconductors under comparable irradiation conditions. The nodal multigap state exhibited by pristine RbCa$_2$Fe$_4$As$_4$F$_2$ is also robust against proton irradiation, with a two-band $d-d$ model being the one that best fits the experimental data.


I. INTRODUCTION
Ion irradiation has proven to be an extremely useful tool for the investigation of fundamental properties of superconductors [1][2][3].In recent years, significant efforts have been devoted to the study of iron-based superconductors (IBSs) [4][5][6][7], hoping that understanding their pairing mechanism could lead to a comprehensive understanding of unconventional superconductivity [8].Although this high goal is still far from being achieved, several interesting findings were reported on these intriguing materials, contributing to deepen the knowledge of the effects of disorder on the superconducting condensate and related phenomena.In particular, besides the more standard result of enhancing the vortex pinning capability [7], an irradiation-disorderinduced transition was observed between s ± and s ++ states in two-band IBSs of the 122 family [3], that contributed to confirm the s ± state for the pristine compound, as predicted theoretically [1].
In the last few years, the attention moved to a newly discovered family of IBSs, the 12442 pnictide family [9][10][11], that shows some peculiar properties resembling those of the cuprates [12][13][14][15].The 12442-type ACa 2 Fe 4 As 4 F 2 (where A = K, Rb, Cs) compounds are obtained by the intergrowth of 1111-type CaFeAsF and 122-type AFe 2 As 2 .While in most IBSs superconductivity emerges after the suppression of antiferromagnetic order in the parent compound by carrier doping, the 12442 compounds are superconducting in their stoichiometric state, with a critical temperature ranging from 28 to 33 K, without the need for chemical substitution to achieve extra carrier doping.These compounds consist of alternate stacking of conducting Fe 2 As 2 layers and insulating Ca 2 F 2 blocks, and contrary to other IBSs, they have double FeAs layers between neighboring insulating layers, thus mimicking the case of double CuO 2 layers in the cuprates.As the cuprates (and contrary to 122 and 1144 IBSs), these compounds show a large anisotropy in the upper critical field and penetration depth [16].As for the pairing symmetry and gap structure, to date it is still controversial whether in 12442 there exist gap nodes or not, and, in the former case, whether nodal d-wave-like or accidental nodal s-wave gaps should be expected.Recently, a study of the gap structure of RbCa 2 Fe 4 As 4 F 2 (Rb-12442) led to the identification of two gaps with signatures of d-wave-like nodal behavior [16].Data are well described by a two-band d−d state with symmetry-imposed nodes, since they persist upon Ni doping.Within this framework, the need emerges to investigate the role of disorder of nonchemical origin, e.g.induced by ion irradiation.In fact, in this case it is possible to compare measurements on the very same crystals before and after irradiation, i.e. after the introduction of structural defects with defined shape and dimensions, and controlled density.This could help in completing the picture and give hints toward an overall interpretation of the superconductivity mechanism in the 12442 IBSs.
In this work, we report on the effects of disorder introduced via 3.5 MeV proton irradiation on the critical temperature, superfluid density and gap values of Rb-12442 single crystals, investigated by a combination of the complementary techniques of electric transport measurements, coplanar waveguide resonator (CPWR) measurements and point-contact Andreev-reflection spectroscopy (PCARS).Transport measurements allow the determination of the resistivity as a function of temperature, of the residual resistivity ratio RRR, and of the superconducting transition tempera-ture T ρ c , and thus provide information on how irradiation affects the scattering of electrons by defects and the onset of the superconducting state.CPWR measurements allow a direct determination of the London penetration depth λ L as a function of temperature, and consequently of the superfluid density ρ s and of the critical temperature T λ c at which it vanishes.The temperature dependence of ρ s provides indications on the symmetry of the superconducting gap(s), and in particular on whether unpaired quasiparticles persist down to the lowest temperatures, which is associated to an incomplete gapping of the Fermi surface typical for example of d-wave superconductivity.PCARS allows a direct determination of the amplitude and symmetry of the superconducting gap(s) and of the temperature at which they close, T A c , by measuring the differential conductance d I/d V of a small (point-like) contact between a normal metal (N) and the superconductor under study (S), as a function of the bias voltage V applied across the NS junction [17][18][19].The absence of an insulating barrier between the N and S banks enables Andreev reflection at the interface [17][18][19], which enhances the d I/d V for V ≤ ∆/e (∆ being the superconducting gap) and makes it sensitive to the coherence of the superconducting condensate -unlike tunnelling spectroscopy that probes the density of states of unpaired quasiparticles.Finally, the dependence of the critical temperature (T ρ c or T λ c ) on the disorder (which is directly related to the symmetry of the energy gap in single-band systems) can be compared to that observed in other, more widely-studied IBS systems.
The conclusion of this extensive study is that the critical temperature of 12442 decreases on increasing disorder in an unusually weak (if compared to other IBSs) and nonlinear manner, while the London penetration depth increases in an approximately linear way.The joint analysis of PCARS spectra and CPWR measurements supports the existence of nodes in at least one gap, and most probably indicates the persistence of the same gap symmetry observed in the pristine, non-irradiated compound.

A. Crystal growth and basic characterization
Single crystals of Rb-12442 were grown by the self-flux method [14,20,21] and their high quality was assessed through X-ray diffraction and energy-dispersive X-ray spectroscopy as detailed in Refs.[14,16].

B. Electrical resistivity measurements
The temperature dependence of the electrical resistivity was measured on selected single crystals by means of the van der Pauw method [22] after having electrically contacted them with thin gold wires and conducting silver paste.Each resistance measurement was performed by sourcing a small DC current ≈ 100 µA between two adjacent contacts and measuring the resulting voltage drop across the opposite pair of contacts.Common-mode offsets including thermoelectric voltages were removed using the current-reversal method.

C. CPWR measurements
The London penetration depth and superfluid density were measured by means of a coplanar waveguide resonator (CPWR) technique, particularly suitable to study small IBS crystals for irradiation experiments since it is minimally invasive and the same samples can be characterized several times before and after different irradiation steps [23].The measurement was carried out within a resonator perturbation approach.A coplanar waveguide resonator obtained by patterning an YBa 2 Cu 3 O 7−δ thin film deposited on MgO substrate was coupled to a vector network analyzer for the measurement of the scattering matrix at low input power.The transmitted power shows a resonance that can be fitted by a modified Lorentzian function, yielding the resonance frequency and quality factor.Measurements were carried out as a function of the temperature for the empty resonator and with the investigated samples coupled to it (i.e.placed on the central stripline, far from the edges).The presence of the sample causes shifts of the resonance frequency and of the quality factor that are due to the electromagnetic properties of the sample: from these shifts, after a self-consistent calibration procedure, the absolute value of the penetration depth and its temperature dependence can be assessed [24].

D. PCARS measurements
The structure of the superconducting gap was assessed via point-contact Andreev reflection spectroscopy (PCARS), by measuring the differential conductance (d I/dV ) of pointlike contacts made between a normal metal (N) and the superconducting sample (S).Unlike in tunnelling spectroscopy, in a direct N/S contact without insulating barrier, electrons impinging on the N/S interface are either transmitted from the metal to the superconductor -if their energy eV is higher than the gap ∆ in the S side -or experience Andreev reflection, i.e. are transmitted as a Cooper pair in S while a hole is reflected back in N [17][18][19].As a result, the differential conductance of the junction is enhanced for V smaller than ∆/e.In real junctions, a small potential barrier can be present at the interface, which makes a tunneling contribution be present in addition to the Andreev-reflection signal.The shape of the resulting d I/dV vs. V curve depends not only on the gap amplitude, but also on its symmetry in the reciprocal space.In particular, isotropic gaps give rise to conductance curves with a zero-bias minimum and two symmetric maxima at the gap edges, no matter what is the direction of current injection with respect to the crystallographic axes [18].Instead, anisotropic nodal gaps (especially in the case of a d-wave symmetry) can give rise to a variety of different curves depending on the angle of current injection [25].For example, zero-bias cusps or peaks are typical features that can be associated to a d-wave gap, in specific conditions of barrier height and current direction [19].In multigap systems, features such as peaks or shoulders are observed in the differential conductance, which are the hallmarks of the different gaps [19].
Point contacts were fabricated by using the so-called soft technique [18,19], i.e. by using a thin gold wire ( ≈ 12.7 µm) stretched over the platelet-like sample and touching its very thin and flat side surface in a single point.These contacts were (optionally) mechanically stabilized by small ( ≈ 50 µm) drop-cast droplets of silver paste.In either case, the "macroscopic" contact actually consists of several parallel nanoscopic N/S contacts.
The differential conductance was obtained by numerical derivation of the I − V characteristics of the contact, measured in the pseudo-four-probe configuration as detailed in Refs.[12,16].Within the present paper, the direction of (main) current injection was confined along the ab planes of the crystals by making the point contacts on the thin, but smooth and flat side surfaces of the crystals.For the d I/dV vs. V curves to contain spectroscopic information, the (maximum) energy of the injected electrons must be directly related to the bias voltage applied through the junction, i.e.E = eV .This requires that the resistance of the contact is much larger than the resistance of the normal bank, and that electrons do not lose energy while crossing the interface.The latter condition means that there is no Joule dissipation in the contact, and is fulfilled if the individual nanocontacts are in the ballistic or at most diffusive regime [18,26], i.e. their radius is smaller than the (inelastic) electronic mean free path.
Information about the gap number, amplitude and symmetry were extracted from the experimental d I/dV vs. V curves of spectroscopic contacts by fitting them with suitable models.Before fitting, the experimental curves at any temperature were normalized to the normal-state curve measured just above the superconducting transition -as customary in IBSs [18,27,28] -after this was rescaled by subtracting the contribution of the so-called spreading resistance R s , i.e. normal-state resistance of the crystal that appears, in series with the contact, in the proximity of the resistive transition and results in a horizontal stretching and a downward shift of the d I/d V spectra with respect to those at low temperature [29].As discussed in Refs.[12,16], for each contact we typically selected the value of R s that allows matching the tails of the normal and superconducting d I/dV spectra.Since this choice remains somewhat arbitrary, we used multiple reasonable values of R s to obtain the normalized spectra, and repeated each time the fitting procedure.The resulting values of the fit parameters were then averaged and their maximum differences taken as the uncertainty of the procedure.
The normalized d I/d V spectra were then fitted to the two-band version [18,19] of the appropriate anisotropic Blonder-Tinkham-Klapwijk model [30][31][32], as extensively discussed in Refs.[12,16].The reason why we used an effective two-band model even if we expect more than two gaps to be present in this compound is due to the necessity to limit the number of adjustable parameters.Indeed, for each gap i the model includes as adjustable parameters the gap amplitude ∆ i , the broadening parameter Γ i , and the barrier parameter Z i .The model also contains the relative weight of the large-gap band in the conductance w 1 and the direction α between the direction of (main) current injection and the antinodal line.The latter is crucial in reproducing zero-bias peaks in the experimental spectra arising from the interference effects occurring when electron-like and hole-like quasiparticles experience order parameters of opposite sign.Being based on the weak-coupling BCS theory, the model does not capture strong-coupling features such as electronboson-derived shoulders appearing at energies higher than the largest gap, which must therefore be excluded from the fitting procedure [19,27,29].

E. Ion irradiation
Ion irradiation was carried out at the CN facility of the Legnaro National Laboratories of INFN, at room temperature, in high vacuum, with a defocused 3.5 MeV proton beam incident parallel to the c axis of the samples.This energy ensures that ion implantation is avoided in all the investigated samples.Beam current was kept below 60 nA on a spot of 7 mm of diameter, to avoid excessive sample heating.
Irradiation with MeV-protons in IBSs is known to produce point-like defects (Frenkel pairs) and small clusters of nanometric size [33], that in the present case are quite uniformly distributed in the sample, since the Bragg peak is avoided.Such defects are efficient scattering centers that can strongly affect the superconducting state by suppressing the superfluid density.
In order to compare the effects of irradiation on the superconducting properties of our samples with other materials and irradiation conditions we computed the displacements per atom (dpa) through Monte Carlo simulations with the SRIM code [34].This approach has proven useful to study the effects of ion-irradiation-induced disorder because dpa is a good parameter to characterize the amount of scattering centers introduced in the sample [35].The dpa value is proportional to the proton fluence, that was evaluated through the total deposited charge obtained by integrating the measured beam current at the beam spot on the sample.

III. RESULTS
We first investigate the effects of ion irradiation on the electric transport properties of the Rb-12442 crystals.residual resistivity ratio, defined as RRR = ρ(300 K)/ρ(T on c ), is RRR = 20 ± 3 and larger than the literature values (12 − 16) [14,16,36]; whereas the resistivity at room T is ρ(300 K) = 290 ± 30 µΩ cm and lower than the previouslyreported values (≈ 1.1 mΩ cm) [14].Both features attest the high crystalline quality of our present samples.Here, T on c is defined as the onset of the superconducting transition, i.e. the temperature at which the resistivity starts to deviate from a linear extrapolation of the normal state (highlighted by the arrows in Figure 1B).The resistive superconducting transition temperature is defined as , where T α c is the temperature at which the resistivity reaches α% of ρ(T on c ), and is equal to T ρ c = 31.3± 0.2 K for the pristine crystal, in good agreement with the previous reports [12,14,16].
At the increase of the irradiation level, the resistivity increases in the entire temperature range up to ρ(300 K) = 470 ± 50 µΩ cm in the crystal at 5.3×10 −3 dpa, and at the same time both T ρ c and RRR are suppressed.As shown in Figure 1C, the suppression of T ρ c introduced by irradiation is significantly faster at lower dpa than at higher dpa values, and for both quantities it is minute with respect to those introduced by substitutional Ni doping.Specifically, we observe a maximum T ρ c suppression of 4.2% and a minimum RRR = 9 in our irradiated samples, against a complete disappearance of superconductivity and a minimum RRR ≈ 1 in Ni-doped samples.This indicates that the perturbations to the structural and electronic properties of Rb-12442 crystals introduced by irradiation are minimal with respect to those attained by substitutional doping.
A more extensive characterization of the superconducting properties of the Rb-12442 crystals upon increasing ion irradiation was carried out by analyzing the data from CPWR measurements.Figure 2A shows the temperature dependence of the London penetration depth λ L for the pristine crystal and for five irradiated crystals with increasing dpa values, from 1.7 up to 7.8×10 −3 dpa.It can clearly be seen that ion irradiation increases λ L in the entire T range and suppresses the onset temperature below which superconducting behaviour is observed.In particular, the superconducting critical temperature T λ c can be obtained as the value of T at which λ L diverges (or equivalently, at which the superfluid density ρ s = λ −2 L goes to 0). Figure 2B shows how this definition correlates extremely well to the midpoint of the resistive transition in transport measurements, for both pristine and irradiated samples.

Figure 2C tracks the evolution of T λ
c upon increasing irradiation level, which confirms both its small suppression as compared with those attained via substitutional doping (maximum T λ c suppression of 4.5% at 7.8×10 −3 dpa) and the progressive reduction in its suppression rate at higher dpa values.Figure 2D shows λ L measured at T = 10 K and the corresponding values of ρ s as a function of the dpa values: λ L (10 K) increases from 211 to 790 nm from the pristine crystal to the crystal irradiated at 7.8×10 −3 dpa, corresponding to a reduction in ρ s (10 K) from 22.4 to 7.8 µm −2 .This shows that ion irradiation is much more effective at tuning the magnetic penetration depth of the Rb-12442 crystals with respect to their superconducting transition temperatures, also when compared with Ni substitutional doping: for instance, 5% Ni-doped crystals -where T λ c was suppressed by ≈ 34% with respect to the pristine crystals -exhibited an increase of λ L (10 K) only up to ≈ 600 nm [16].Now we turn our attention to the effect of ion irradiation on the values of the superconducting energy gaps as measured by PCARS.Figures 3A, 3B and 3C show some ex- amples of low-temperature PCARS spectra of different point contacts (with different normal-state resistance) made on the sample with a disorder level of 5.3×10 −3 dpa.Most spectra (Figure 3A) show a zero-bias cusp or maximum, which is suggestive of a nodal gap [18,19,[30][31][32], and is perfectly consistent with the observations made on pristine Rb-12442 [12,16].Others, instead, show a broad maximum (Figure 3B) or even two symmetric maxima at low bias (Figure 3C) that are usually associated to an isotropic gap, but can be observed even in d-wave gap symmetry when the current is injected (mainly) along the antinodal direction [18,19,[30][31][32].
More detailed information about the energy gap can be extracted by normalizing the experimental spectra to the normal-state curve measured immediately above T c [18,27,28], and fitting the normalized spectra to a suitable model for Andreev reflection at a normal metal/superconductor interface [18,19].Figures 3D, 3E, and 3F show three different normalized spectra (black lines).None of these could be fitted by using a single gap, either in s-wave or d-wave symmetry, indicating that irradiation does not change the multiband nature of superconductivity in this compound.We then used an effective two-band model, in which the relative weight of the gaps is an adjustable parameter.Based on the results found in pristine, unirradiated Rb-12442 [12,16] and on the fact that some of the curves (Figure 3D) display a zero-bias cusp, we tried to fit the curves using either an isotropic s-wave gap and a nodal d-wave gap (s−d model in the following) or two d-wave gaps (d−d model).As discussed in Ref. [12], the s−d model actually mimics the nodal s ± and the d x 2 − y 2 symmetries, whereas the d−d model mimics the d x y symmetry, all compatible with the tetragonal structure of Rb-12442 [8].In both the models, the partial conductance associated to each gap is modelled as in Ref. [32].The fitting parameters are then, in addition to the gap amplitudes ∆   Z 2 , the broadening parameters Γ 1 and Γ 2 , and the angle α between the direction of (main) current injection and the antinodal direction of the d-wave gap.In the case of the d−d model, we assume α 1 = α 2 , that implies that both the gaps have nodes in the same direction.
It turns out that both the models allow a good fit of the curves, as already pointed out in the case of pristine Rb-12442 and of the Ni-doped compound [12,16].The fact that irradiation does not change the gap structure is compatible both with the small effect on the critical temperature, and with the observed persistence of the same structure even in Ni-doped samples [12,16] where the dopants introduce disorder and give rise to a much more significant T c suppression (about 10 K).It is worth noting that, in order to reproduce the shape of the curves in panels E and F, the current injection has to be assumed to occur along the antinodal direction of the d-wave gap (i.e.α ≃ 0).Just as an example, Figures 3D, 3E and 3F show the results of the d−d fit of three different curves measured at the lowest temperature.The fit is very good in the central part of the curves, while for higher energies (starting at about 10 meV) it fails due to an excess conductance probably due to the coupling of charge carriers with some bosonic mode involved in the superconducting pairing, as always observed in IBSs [19,27,29,37].The amplitudes of the superconducting gaps vary from one contact to another, probably due to the fact that the compound features more than two gaps, while we are using here a minimal (effective) two-gap model in order to keep the number of adjustable parameters to a minimum; moreover, the presence of the bosonic structures prevents an accurate determination of the larger gap, presumably causing its overestimation.The same hap-pened in pristine Rb-12442, where the gap amplitudes were spread over a rather wide band of energies.In particular we found, for the d−d model, ∆ d1 ∈ [4.8, 7.3] meV and ∆ d2 ∈ [1.3, 3.0] meV.The gap amplitudes that result from the d−d fit of the curves in irradiated Rb-12442 are still generally compatible with the same distribution of gap values, meaning that no significant change in the gap amplitudes has been produced by irradiation.The same happens in the case of the s−d fit.
Figure 3G reports, as an example, the temperature dependence of the conductance curves of a point-contact whose normal-state resistance is 48.7 Ω. Above 31 K (black curve) the conductance curves are smooth and do not show any feature at zero bias.This means that, at this temperature, the transition to the normal state is complete (in accordance with the green curve in Figure 1B).We thus chose this curve as the normal-state conductance and used it to normalize all the others at lower temperatures, with no need of any further rescaling or shift.The resulting normalized curves are shown in Figure 3H (symbols), vertically offset for clarity (apart from the lowest-temperature curve).The best-fitting curves, obtained within the d−d model, are shown as well as solid lines.The fit follows rather well the temperature evolution of the experimental curves; up to 23 K, it was possible to fit them by keeping all the parameters of the model (apart from the gaps) fixed to their lowest-temperature values, which is always a good indication of both the reliability of the curves (and of the normalization) and of the ability of the model to grasp the physics of the system.In particular, we used w = 0.9, Γ 1 = 0.45 meV, Z 1 = 0.14, Γ 2 = 0.25 meV, Z 2 = 0.32.Above 23 K, instead, the conductance curves start to shrink, becoming too narrow with respect to what predicted by the model; if one wants to reproduce this shrinking, the values of the barrier parameters and of the broadening parameters must be reduced.The possible reason for this behaviour can be understood by plotting the gap amplitudes resulting from the fit as a function of the temperature, as in Figure 3I.This plot shows that the small gap amplitude is rather robust against the model used (either s−d or d−d) which is natural since in either case the small gap has the d-wave symmetry.The amplitude of the larger gap instead depends on whether it is s-wave or d-wave, which is perfectly natural due to the different shape of the Andreev-reflection features that these two symmetries produce.The most relevant thing is that, no matter what the model is, the gap amplitudes follow a very clear trend which can be approximated by the solid and dashed curves, whose expression is where ∆ i is the gap of the i-th band and T A c = 29 K.This indicates that the behaviour is not exactly BCS, and that the Andreev critical temperature T A c (i.e. the temperature at which the Andreev signal disappears, and at which the superconducting bank of the junction becomes normal) is about 29 K, and thus smaller than T c = 31 K at which the narrow zero-bias maximum in the experimental curves dis-appears (see Figure 3G).As is clear from Figure 1B, T A c falls in the proximity of the superconducting transition, but still in a region where the resistivity is zero.This is likely to be Temperature dependence of the zero-bias resistance extracted from the curves in Figure 3G (blue hollow circles, left scale) compared to the resistivity of the same crystal (solid green line, right scale).The inset shows the zero-bias conductance (filled blue circles), that follows the same temperature dependence of the gaps up to 29 K (solid black line).due to the non perfectly ballistic nature of the contact, so that a very small heating occurs in the contact region giving rise to a small mismatch between the "bulk" T ρ c and T A c .The same effect explains the narrow zero-bias peak that emerges at high temperature in the conductance curves: while at zero bias the temperature of the contact is equal to that of the bulk, on increasing the bias voltage the current gives rise to heating and, being close to the normal transition, this produces an increase in resistance (and thus a decrease in conductance).Such an effect ceases when the temperature at which the curve is acquired is already in a region where the normal transition (as measured by transport) is complete, since here the resistivity has a much weaker dependence on temperature.Figure 4 displays the temperature dependence of the zero-bias resistance (main panel) and conductance (inset).It is clear that the high-temperature behaviour of both these quantities is no longer due to the temperature dependence of the superconducting gap, but is related instead to the temperature dependence of the crystal's resistivity.

IV. DISCUSSION
The response to proton-irradiation-induced disorder of Rb-12442 is qualitatively as expected for the IBSs, showing a suppression of the critical temperature and of the superfluid density due to the enhanced scattering of carriers in a multiband system.However, there is a relevant quantitative difference: below about 4 − 5×10 −3 dpa, the T c suppression rate with dpa in optimally-doped (Ba,K)Fe 2 As 2 , filled blue circles), from dc resistivity (T ρ c , hollow violet circles) and from PCARS measurements (T A c , filled green diamonds -the smaller value is due to localized contact heating as explained in the text).The dashed line is a guide for the eye.Data are normalized to the corresponding transition temperature at 0 dpa (T c,0 ) and compared to the trend (black and red lines) showed by two optimally-doped systems of the 122 IBS family (K and Co-doped BaFe 2 As 2 ), for which the scaling laws with dpa were reported [35].(B) Proton irradiation-induced increase of the London penetration depth λ L at T = 10 K for Rb-12442 samples (symbols, dashed line is a guide to the eye), compared to the trends shown by the same two IBS Ba-122 systems shown in (A) [35].
is about twice than for Rb-12442 samples, despite the similar pristine critical temperature of the two compositions.This is reported in Figure 5A, where the data obtained in this work on Rb-12442 are compared to the typical behavior exhibited by the 122 family of IBSs when exposed to the same irradiation process [35].It is also evident that the T c vs dpa curve of Rb-12442 shows an upward curvature (or alternatively, that it can be divided in two separate linear trends with different slopes), implying an even lower dpa dependence at higher disorder levels.This behavior is quite unexpected, and could -in principle -be the result of an important change in the symmetry of the superconducting gap of the system.However, the fact that the experimental PCARS spectra measured on a sample with a disorder level of 5.3×10 −3 dpa can be satisfactorily fitted using models employing the same gap symmetry describing pristine and Ni-doped Rb-12442 [12,16] makes this scenario highly unlikely.Similarly to T c , also the low-temperature value of the London penetration depth, as shown in Figure 5B, has a much smaller variation with dpa when compared to other IBSs.In this case, the modification with dpa is approximately linear in the whole investigated range.
These two observations indicate that the condensate in this system is affected much less than similar systems by the exposition to the same proton irradiation conditions.Since the composition, density and chemical nature of the 12442 and 122 families of IBSs is similar, a similar defects morphology is induced in the samples.Therefore, either the recombination rate of the created defects is much higher in the 12442 family (possibly due to the more layered and anisotropic structure), or the superconducting state itself is more radiation-hard.The latter option could be related to the other peculiarities observed in this system, such as the observation of nodal behavior [16,[38][39][40] and its vicinity to the disappearance of the electron-like pockets [12].
A further confirmation of the robustness of nodal behavior in Rb-12442 against the introduction of disorder via proton irradiation can be obtained from the temperature dependence of the superfluid density.Figure 6A shows ρ s = λ −2 L as a function of T determined from the λ L (T ) curves shown in Figure 2A: beyond the suppression of the superfluid density in the entire temperature range with increasing dpa, it can be seen that all the ρ s (T ) curves exhibit a smooth, kinkless behavior in the experimentally-accessible temperature range, and a clear linear dependence at low temperatures.Both features are consistent with the previous reports in pristine and Ni-doped Rb-12442 crystals [12,16] and are typically associated with multigap superconductivity with mainly interband coupling [41] and with the existence of at least one nodal gap [38,39] respectively.This interpretation can be further corroborated by combining the results obtained independently by means of the CPWR and PCARS measurements.As a matter of fact, the superfluid density can also be computed from the zerotemperature London penetration depth, λ L (0), and from the temperature-dependence of the gap values extracted from PCARS, ∆ i (T ), as [42]: where i identifies the band, ∆ i (φ, T ) is the angle-dependent superconducting gap function, f = [1 + exp(E/k B T )] −1 is the Fermi function and w i is the mixing weight of the i-th gap contribution (constrained by w 1 + w 2 = 1).For this latter parameter we use the same value previously employed for pristine undoped Rb-12442 [16] and consider it as a fixed parameter.One can factorize the angular and temperature dependencies of the gap as ∆ i (φ, T ) = ∆ 0,i f φ (φ) f T (T ), where the angular function for d-wave (nodal) superconductors is f φ (φ) = cos(2φ) and the gap amplitude is the one measured by PCARS.To obtain the λ L (0) value for the disorder level of the sample measured by PCARS we linearly interpolate the data shown in Figure 5B.
The comparison between the superfluid data from PCARS obtained within this approach and the data from the CPWR technique is shown in Figure 6B, for the disorder level of 5.3×10 −3 dpa.Violet (green) crosses indicate the values of ρ s calculated by starting from the PCARS gaps extracted from the s−d (d−d) fit, while the black dotted line represents the ρ s (T ) curve deduced from the CPWR analysis for a disorder level of 5.3×10 −3 dpa.The PCARS d−d case seems to yield the best qualitative match with the CPWR trend since it reproduces the linear behavior at low temperatures (even though with a lower slope), contrarily to what is found for the s−d case, where the s-wave component forces the curve to saturate.
In summary, we reported on the unusually weak irradiation effects on the superconducting properties of the IBS Rb-12442.The analysis was performed by CPWR and PCARS on single crystals before and after 3.5-MeV proton irradiation, at different damage levels up to 7.8×10 −3 dpa.Both the critical temperature and the superfluid density are shown to decrease with disorder, but at a much lower rate than for other IBS systems with similar pristine values of T c and ρ s .Moreover, the critical temperature degradation shows also a qualitatively different trend with respect to the other mentioned IBSs, i.e. a change of linear slope or possibly an upward curvature, implying an even smaller disorder dependence for higher irradiation doses.A comparison between results from the two measurement techniques allowed us to elaborate about the order parameter symmetry, in an effective two-gap approach.The evidence of a nodal behavior is confirmed also for the irradiated crystal and on the basis of the comparison we are inclined to select the d−d scenario rather than the s−d one, even though an ultimate evidence is still needed, which encourages further investigation.

FIG. 1 .
FIG. 1. (A) Resistivity ρ as a function of temperature T for three Rb-12442 crystals in the pristine state (corresponding to dpa = 0) and at increasing proton irradiation levels.(B) Magnification of the curves shown in (A) in the vicinity of the superconducting transition.The arrows highlight the onset temperature of the resistive transition T on c (see text for details).(C) Superconducting transition temperature T ρ c (filled violet circles, right scale) and residual resistivity ratios RRR (hollow green circles, left scale) of the same crystals as a function of the disorder level measured in displacements per atom (dpa).Dashed lines are guides to the eye.

FIG. 2 .
FIG. 2. (A) Temperature (T ) dependence of the London penetration depth λ L determined from CPWR measurements for pristine and irradiated crystals for increasing proton irradiation levels, measured as displacements per atom (dpa).(B) Comparison of the T dependencies of λ L (hollow circles, left scale) and of the electrical resistivity ρ (solid lines, right scale) measured in the same two samples (a pristine crystal and an irradiated crystals at 2.8×10 −3 dpa).The superconducting critical temperature determined from CPWR measurements, T λ c , (defined as the T at which λ L diverges) correlates well with the midpoint of the resistive transition.(C) T λ c and (D) λ L (hollow green circles, left scale) and ρ s (filled violet circles, right scales) measured at 10 K determined from the curves shown in panel (A) as a function of the disorder level measured in dpa.Dashed lines are guides to the eye.
s and ∆ d (in the s−d model) or ∆ d1 and ∆ d2 (in the d−d model), the weight w 1 of the large gap (such that w 2 = 1 − w 1 ), the barrier parameters Z 1 and

FIG. 3 .
FIG. 3. (A),(B),(C) Some examples of as-measured conductance curves (d I/dV vs V ) of ab-plane contacts made on the irradiated crystal for which the disorder level is 5.3×10 −3 dpa.The normal-state resistance of each contact is indicated in the labels.(D),(E),(F) Examples of normalized conductance curves (black lines) with the relevant d−d fit (blue lines).The labels indicate the temperature at which the spectra were acquired, the value of the weight w 1 given by the fit, and the amplitude of the two gaps.(G) Temperature dependence of the raw conductance spectra of the contact with R N = 48.7 Ω.The black curve (at T = 31 K) is the first that does not show zero-bias features and corresponds to the attainment of the bulk normal state.(H) The same conductance curves of panel (G) after normalization (symbols) and the relevant best-fitting curves (lines) calculated within the d−d model.(I) Temperature dependence of the energy gaps extracted from the fit of the curves in panel (G) using the s−d model (filled circles) and the d−d model (hollow circles).Solid lines are functions of the form of Equation (1).
FIG. 4.Temperature dependence of the zero-bias resistance extracted from the curves in Figure3G(blue hollow circles, left scale) compared to the resistivity of the same crystal (solid green line, right scale).The inset shows the zero-bias conductance (filled blue circles), that follows the same temperature dependence of the gaps up to 29 K (solid black line).

FIG. 5 .
FIG. 5. (A)Proton-irradiation-induced suppression of the superconducting critical temperature T c for Rb-12442 samples, determined from CPWR analysis (T λ c , filled blue circles), from dc resistivity (T ρ c , hollow violet circles) and from PCARS measurements (T A c , filled green diamonds -the smaller value is due to localized contact heating as explained in the text).The dashed line is a guide for the eye.Data are normalized to the corresponding transition temperature at 0 dpa (T c,0 ) and compared to the trend (black and red lines) showed by two optimally-doped systems of the 122 IBS family (K and Co-doped BaFe 2 As 2 ), for which the scaling laws with dpa were reported[35].(B) Proton irradiation-induced increase of the London penetration depth λ L at T = 10 K for Rb-12442 samples (symbols, dashed line is a guide to the eye), compared to the trends shown by the same two IBS Ba-122 systems shown in (A)[35].

FIG. 6 .
FIG. 6. (A) Superfluid density (ρ = λ −2 L ) as a function of temperature T for increasing values of proton irradiation-induced disorder measured in displacements per atom (dpa), from CPWR measurements.(B) Data calculated from the gap values measured by PCARS (crosses) are compared to the temperature dependence of ρ s deduced starting from CPWR data for the disorder level of 5.3×10 −3 dpa.