Spectroscopic evidence of mixed angular momentum symmetry in non-centrosymmetric Ru\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_7$$\end{document}7B\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document}3

Superconducting crystals with a lack of inversion symmetry can potentially host unconventional pairing. However, till today, no direct conclusive experimental evidence of such unconventional order parameters in non-centrosymmetric superconductors has been reported. In this paper, through direct measurement of the superconducting energy gap by scanning tunnelling spectroscopy, we report the existence of both s-wave (singlet) and p-wave (triplet) pairing symmetries in non-centrosymmetric Ru\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_7$$\end{document}7B\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document}3. Our temperature and magnetic field-dependent studies also indicate that the relative amplitudes of the singlet and triplet components change differently with temperature.

www.nature.com/scientificreports/ two parts exhibited significantly different sensitivity to the applied magnetic field and led to two dramatically different field scales for the upper critical fields ( H c2 ∼ 1.1 T and 5 T respectively). In this context, Fang et al. discussed the possibility of a mixed angular momentum symmetry of the superconducting order parameter in Ru 7 B 3 . In addition, though it was ignored by the authors, a possible signature of unconventional pairing was also present in the specific heat data as presented in 23 . More recently, Cameron et al. performed small-angle neutron scattering 26 on Ru 7 B 3 and reported that the orientation of the vortex lattice in Ru 7 B 3 strongly depends on the history of the applied magnetic field, thereby indicating the possibility of broken time-reversal symmetry in the order parameter. The superconducting order parameter is an energy-resolved quantity. Hence, a direct energy-resolved measurement technique, like scanning tunneling spectroscopy, is essential to understand its true nature. In RuB 2 , which is a centrosymmetric superconducting phase of the same elements Ruthenium and Boron, recently such tunneling spectroscopic data revealed unusual multi-band signatures in the order parameter 27 . Motivate by that, in order to probe the true order parameter symmetry of the non-centrosymmetric phase Ru 7 B 3 also, we carried out a detailed temperature and magnetic field-dependent tunneling measurements on Ru 7 B 3 . The analysis of such data reveals the spectroscopic signatures of an order parameter with mixed angular momentum symmetry.

Experimental methods
The single crystals used for our measurements showed a bulk superconducting transition at 2.6 K (Fig. 1b). The scanning tunneling microscopy (STM) and spectroscopy (STS) experiments were performed in a Unisoku system with an RHK R9 controller, inside an ultra-high vacuum (UHV) cryostat at ∼ 10 −10 mbar pressure. The lowest temperature down to which the measurements were performed was 300 mK. The STM is also equipped with a superconducting solenoid capable of producing a magnetic field up to 11 T. Since the single crystals could not be cleaved using the standard cleaving technique (optimized for layered materials only), we cleaned the surface by reversed sputtering for 30 min with Argon (Ar) ion in-situ inside an integrated UHV preparation chamber. Following that, we immediately transferred the sample to the scanning stage at low temperature. The Tungsten (W) tip, which was prepared outside by electrochemical etching, was also cleaned in UHV by bombarding it with a high-energy electron-beam. This process helped us probe the pristine surface of Ru 7 B 3 . In the inset of Fig. 1b, we show an STM topographic image showing the distinctly visible crystallites with average grain size ∼ 3 nm.

Results and discussion
In Fig. 1a we show six representative tunneling spectra ( S 1 -S 6 ) captured at randomly chosen points on the surface of Ru 7 B 3 at ∼ 310 mK. A visual inspection reveals that, the spectra can be distinctively divided into two categories based on their overall shapes. The first type ( S 1 , S 2 , and S 3 ) shows coherence peaks around ± 0.30 mV and shallow www.nature.com/scientificreports/ 'V'-shaped dip in between the peaks. The spectra ( S 4 , S 5 , and S 6 ) of the second type show coherence peaks at ± 0.47 mV, and they have a higher curvature and depth (dI/dV ∼ 0 at V=0) below that. We have also analyzed these spectra within a single band 's-wave' model 1 using Dyne's formula 28 : The tunnelling current is given by Here, N s (E) and N n (E) are the density of states (DOS) of the superconducting sample and the normal metallic tip, respectively and f(E) is the Fermi-Dirac distribution function. Ŵ is the Dyne's parameter that takes care of the broadening of DOS. The theoretical plots thus generated are shown as black lines on the experimental data points in Fig. 1a. The spectra S 1 , S 2 and S 3 give of the order of 0.31 meV, 0.29 meV, and 0.24 meV, respectively while spectra S 4 , S 5 , and S 6 provide the values 0.48 meV, 0.42 meV, and 0.49 meV respectively from the same analysis. To note, while the first group of spectra ( S 1 , S 2 , and S 3 ) shows reasonably good fitting with the single gap 's-wave' model (albeit with large Ŵ ), the second group ( S 4 , S 5 , and S 6 ) exhibits a significant departure from that.
The above-mentioned discrepancy between the experimental spectra and the spectra generated theoretically within a single-band 's-wave' model, prompted us to consider other possible symmetries of the order parameter. To perform such an analysis, we modified Dyne's equation by introducing a more general expression 29 than an isotropic . The modified Dyne's equation reads as N s (E, θ) ∝ Re Here, θ is the polar angle (w.r.t. (001)) and the integer n can be 0, 1 or 2 for s, p, and d wave symmetries, respectively. The expression for tunneling current is also modified to In Fig. 1c we show the experimental spectrum S 4 along with theoretical plots considering isotropic 's-wave' (red line) and anisotropic 'p-wave' (blue line). It is clear that the spectrum, especially the 'V'-shaped part of that between the coherence peaks, is better described by the 'p-wave' symmetry. It is also interesting to note that the extracted value of (0.47 meV) for such fit does not differ much from the best 's-wave' fit (0.48 meV).
Such 'V'-shaped STS spectra are often seen for superconductors with possible unconventional symmetries and are well described by the Tanaka-Kashiwaya model 29,30 we used here. For example, considering various possible symmetries under the framework of an "anisotropic s-wave' model in MgB 2 , Seneor et al. 31 explained spectroscopic signatures, e.g., very high coherence peaks, 'V'-shaped valley, and also zero-biased coherence peak (ZBCP). In anisotropic Bi 2 Sr 2 CaCu 2 O 8 , Ichimura et al. 32 reported two types of spectra, one with a 'V'-shape and another with a flat bottom. The authors explained those with a model of mixed ' s + d-wave' symmetry where the former isotropic component is the dominant one. Millo et al. 33 reported both 'V'-shaped spectra and spectra with ZBCP in SmFeAsO 0.85 and adduced such shapes with the 'd-wave' order parameter. To note, it was also reported there that some of the STM spectra could also be fitted well within a pure 's-wave' model but with significantly smaller and relatively higher Ŵ 33 -a situation similar to our group I spectra. From our analysis, an order parameter with either 'p-wave' or 'd-wave' symmetry can reproduce a 'V'-shaped spectra as mentioned in Group II. However, as reported by Cameron et al. 26 , the spontaneous magnetic field present in Ru 7 B 3 breaks the time-reversal symmetry in the order parameter and supports the 'p-wave' gap as a more favorable possibility. Now we focus on the spectrum S 1 , which belongs to group I and fits reasonably well with a single 's-wave' gap (and with relatively large Ŵ ). A closer inspection, however, reveals that there is a small discrepancy between the experimental data and the 's-wave' model spectrum. We investigated the evolution of this discrepancy with temperature. Temperature dependence of approximately followed the BCS prediction ( Fig. 2c) with 0 = 0.3 meV. The broadening parameter Ŵ did not change much within this range. The departure of the experimental spectrum from the s-wave model rapidly decreased with increasing temperature, and at around 750 mK, the discrepancy almost disappeared (Fig. 2a). To illustrate this effect clearly, we show two spectra, one at 338 mK and another at 961 mK, along with their theoretical (s-wave) fits in Fig. 2b. This observation indicates the possibility of a mixed angular momentum symmetry in the order parameter where the amplitudes corresponding to different l vary differently with temperature. This also explains why Fang et al. did not find any signature of unconventionality in their lower critical field ( H c1 ) studies 24 down to 1.2 K, which is well above the temperature window where we could notice a deviation from the s-wave behaviour in our data. We have also performed magnetic field dependence of another spectrum ( S 3 ) from the same group I (Fig. 3a). The key observation is that while the spectra deviate appreciably from 's-wave' theoretical curves at low fields, above 10 kG the spectra resembled far more closely to the 's-wave' predictions. For demonstrating this effect more clearly, we show the spectra and respective 's-wave' fits both at zero field and at 16 kG, in Fig. 3b. The evolution of the extracted gap ( ) and broadening parameter ( Ŵ ) with the applied magnetic field are shown in Fig. 3c.
In Fig. 4a, we demonstrated a detailed temperature dependent study of spectra S 6 , which belong to group II. At ∼ 740 mK, a small zero-bias conductance peak (ZBCP) appeared. It became more and more pronounced up to ∼ 1370 mK and then started fading until it, along with other prominent spectral features, gradually disappeared at 1.9 K. To note, the emergence of such a ZBCP can be attributed to a 'p-wave' component in the order parameter where the interface normal and the lobe-direction of the 'p-wave' maintain an acute angle between them 34,35 . Since the surface of our Ru 7 B 3 has crystallites with arbitrary orientations, and the tip is engaged randomly at different points, this condition can naturally be satisfied sometimes.
Based on the discussion above, if Ru 7 B 3 has a mixed angular momentum symmetry in its order parameter, then we would expect the unconventional component to be present in group-I spectra too. To understand this aspect in detail, we used an ' s + p-wave' model to analyze the spectra S 1 . In this model, the effective gap is given by s+p = s + � p Cosθ . In Fig. 4b, we present the spectrum S 1 (black circles) at 338 mK with numerically generated spectra considering ' s + p-wave' symmetry (blue line). Pure 's-wave' fit (red line) is also shown for comparison. It is clear that the mixed angular momentum symmetry provides a better description of the data. Temperature evolutions for the amplitudes of the two components s and p extracted from the above ' s + p-wave' fittings for the whole spectra S 1 are presented in Fig. 4c. It is not a surprise that, while the conventional s follows a smooth BCS like dependence up to 2 K, the smaller p sharply drops and goes beyond our measurement resolution at 0.9 K.   www.nature.com/scientificreports/ It should be noted here that the superconducting gap is the manifestation of a phase-coherent macroscopic condensate. Therefore, the measured should ideally be unique irrespective of the measurement techniques.
Since the system is metallic in nature 23,24 , any role of special surface states can be ruled out. The absence of multiple distinct pairs of coherence peaks in any spectra makes the possibility of a multi-band superconductivity or a proximity-induced superconductivity (PIS) in the bottom of the STM tip unlikely. In a different context, a similar variety of transport spectra is possible if some superconducting nanoislands are deposited and separated spatially from each other on an otherwise non-superconducting metal or insulator 36,37 . PIS is reported on such non-superconducting substances in the neighborhood of those superconducting islands, which eventually gets vanished while moving away from the islands. Nevertheless, in our case, the whole surface under the tip is of a pure superconducting Ru 7 B 3 single-crystal, unlike such geometry, and we did not find an entirely gap-less flat spectrum anywhere on the surface. Also, as reported by Cherkez et al. 37 , spectra probed on the superconducting islands and away from the islands, both could be well fitted with the BCS model with just different sets of parameter values. However, in our case, a pure BCS model fails to explain the spectra, and an anisotropic component, whether relatively small ( for group I ) or large ( for group II ) is essential. In that way, we can safely exclude the possibility of any proximity effect in our measurements. The T c and the H c that we measure for all our spectra match well irrespective of whether they fall under Group I or II. Therefore, it can be concluded that all the spectra falling under two distinct groups differ from each other based on which component of the order parameter contributes predominantly for a particular crystallite that the measurement is performed on.

Conclusion
In summary, we performed scanning tunneling spectroscopy on single-crystal Ru 7 B 3 and recorded several spectra, which can be broadly categorized into two groups, group-I and II. Group-I consists of spectra that are shallow in shape. They show overall good agreement with single gap 's-wave' symmetry except at very low temperatures where they deviate from such a pure s-wave model. Smaller superconducting gaps and larger broadening parameters are characteristics of these spectra. Group-II spectra are broader in shape and show sharper coherence peaks. They significantly deviate from the predictions of the 's-wave' model, but a theoretical model with 'p-wave' symmetry shows a better agreement. The superconducting gaps are larger, and the broadening parameters are very small for such spectra compared to those belonging to the first group. The temperature dependences of the spectra belonging to both groups indicate the presence of a mixed ' s + p-wave' symmetry in the order parameter where the two components have different temperature dependence.