# Superconducting and normal state properties of noncentrosymmetric superconductor NbOs investigated by muon spin relaxation and rotation

###### Abstract

Noncentrosymmetric superconductors with -manganese structure has attracted much attention recently, after the discovery of time-reversal symmetry breaking in all the members of ReX (X = Ti, Hf, Zr) family. Similar to ReX, NbOs also adopts - structure and found to be superconducting with critical temperature 2.7 K. The results of the resistivity, magnetization, specific heat and muon-spin relaxation/rotation measurements show that NbOs is a weakly coupled type-II superconductor. Interestingly, the zero-field muon experiments indicate that the time-reversal symmetry is preserved in the superconducting state. The low-temperature transverse-field muon measurements and the specific heat data evidence an conventional isotropic fully gapped superconductivity. However, the calculated electronic properties in this material show that the NbOs is positioned close to the band of unconventionality of the Uemura plot, indicating that NbOs potentially borders an unconventional superconducting ground state.

## I Introduction

The discovery of exotic superconducting properties in the heavy fermion noncentrosymmetric superconductor (NCS) CePtSi EB (), sparked renewed research interest both from both experimental and theoretical perspectives, to understand the role of structural asymmetry in superconductivity. Theoretically, in the Bardeen-Cooper-Schrieffer (BCS) superconductors, the pairing state is protected by the crystal inversion symmetry. When a superconductor has an inversion center, the pairing states can be unambiguously classified into degenerated states of spin-singlet (even parity) and spin-triplet (odd parity). However, this scenario changes in noncentrosymmetric superconductors, where the lack of inversion symmetry induces odd parity antisymmetric spin-orbital coupling (ASOC). A strong ASOC results in mixed parity states and develops complicated spin structures EB1 (); LP (); PAF (). Such an admixture of spin states reflects in a variety of unique properties, such as: anomalous upper critical field values EB2 (); NK (); IS (), line nodes IB (); EB3 (), and magnetoelectric effects VM (); SF () etc. Besides, recently it was proposed that some NCS can also be a potential candidate for topological superconductivity, of which known examples are PbTaSe GB1 (); SYG () and BiPd ZSM (); MN (), which shows topological surface states and potentially lead to signatures of Majorana fermions.

Since CePtSi, several more heavy fermion superconductors have been found to exhibit a wide variety of unusual superconducting properties EB1 (). But it was soon realized that heavy fermion noncentrosymmetric superconductors are not the simplest of systems to study. In these superconductors, additional complications come due to the presence of strong electronic correlation effects and quantum criticality, which often hinders the research aimed to understand the interplay between crystal symmetry and superconductivity. As a result, many new noncentrosymmetric superconductors with weak correlation have been targeted to look for singlet-triplet mixing. Some of these compounds are LaNiC ADH (); ADH1 (), LiPdB KTP (); PBT (); PBT1 (); MNY (); HTK (); HQY (), LiPtB PBT1 (); MNY (); HTK (); HQY (); MNY1 (); HTM1 (), RuB LFH (), NbRe ABK (), ReW PKB (), YC JCM (), MoAlC MLC (); MLC1 (). A particularly interesting example is LiB (M = Pd, Pt) HTK (); HQY (); MNY1 (); HTM1 (); TMH (); RKI (); MMD (); HQY1 (); STT (); ASY (); DCP (); GED (). LiPdB exhibits a isotropic fully gapped superconductivity MNY (); HTK (); HQY (), whereas the gap of LiPtB has line nodes HQY (); MNY1 (); HTM1 (). Here the substitution of Pd with Pt element strengthens the ASOC and enhances the relative pairing mixing ratio with a dominant spin-triplet component HQY (). Hence, the change from fully gapped to nodal superconductivity is directly dependent on the magnitude of ASOC. Indeed, NMR experiments also conclude that the Cooper pair included about 60 of the spin-triplet pairing in LiPtB.
Accordingly, many noncentrosymmetric superconductors with wide variety of SOC have been studied, yet most of them such as LaMSi (M = Pd, Pt)MSA3 (), LaMSi (M = Rh, Ir)VKA1 (); VKA2 (), BiPd KMS (); XBY (), NbRe ABK (), ReW PKB () found to be compatible with fully gapped superconductivity. Thus, the importance of ASOC strength in the degree of mixing ratio is not fully understood, and therefore requires to find more examples where the effects of modified spin-orbit coupling on the singlet-triplet mixing can be studied directly.

Another noteworthy feature of unconventional superconductivity is time-reversal symmetry (TRS) breaking, which is an extremely rare phenomenon has been found only in few noncentrosymmetric superconductors, such as: LaNiC ADH (), LaIr JAT (), SrPtAs PKB1 () and ReX (X = Ti,Hf,Zr) RPS (); DSJ (); DSJ1 (). It is the latter family of compounds that are particularly interesting owing to the frequent occurrence of TRS breaking among all the members. In ReX, the TRS breaking effects seem very similar irrespective of the substitution at X-sites, suggesting a negligible effect on the ASOC and hints towards its common origin. Interestingly, the band-structure calculations of the Re-based binary alloys r1 (); r2 (), indicate that the density of states (DOS) at the Fermi level to be dominated by the d-bands of Re and could be the critical factor leading to TRS breaking. Indeed, recent observation of TRS breaking in pure Re metal (centrosymmetric), strongly suggests that the sizable Re spin-orbit coupling to be the origin of TRS breaking in the ReX family TMC ().
However, to strengthen the above conclusion one need to study other Re-free materials, whilst still having the same - type structure. In this context, we report a detailed investigation of the superconductivity of NCS NbOs, isostructural to the ReX family. NbOs exhibits bulk superconductivity around 2.7 K BT (); BT1 (), characterized via electrical resistivity, specific heat, magnetic susceptibility, and muon spectroscopy. The principal goal of the present work is to study the superconducting properties of NbOs and search for TRS breaking in a Re-free material with - structure.

## Ii Experimental Details

The polycrystalline samples of NbOs were prepared by arc melting a stoichiometric mixture of Nb (99.95 ) and Os (99.95 ) under high purity argon gas atmosphere on a water cooled copper hearth. The as-cast ingot was flipped and remelted several times to ensure the phase homogeneity. The observed mass loss was negligible. The sample was then sealed inside an evacuated quartz tube and annealed at 800 C for one week. To verify the phase purity we performed room temperature (RT) powder x-ray diffraction (XRD) using a X’pert PANalytical diffractometer (Cu-K radiation, = 1.540598 Å). The superconducting properties of NbOs were measured using magnetization , ac susceptibility , electrical resistivity , specific heat , and muon relaxation/rotation (SR) measurements. The dc magnetization and ac susceptibility data were collected using a Quantum Design superconducting quantum interference device (SQUID). The electrical resistivity and specific heat measurements were performed on a Quantum Design physical property measurement system (PPMS). The SR measurements were conducted at the ISIS Neutron and Muon facility, in STFC Rutherford Appleton Laboratory, United Kingdom using the MUSR spectrometer. The powdered NbOs sample was mounted on a silver holder and placed in a sorption cryostat, which we operated in the temperature range 0.3 K - 3.0 K. Zero-field muon spin relaxation (ZF-SR) and the transverse-field muon spin rotation (TF-SR) measurements were performed on the MuSR beam line at the ISIS pulsed muon source. A full description of the SR technique may be found in Ref. SLL (). In ZF-SR, the contribution from the stray fields at the sample position due to neighbouring instruments and the Earth’s magnetic field is cancelled to within 1.0 T using three sets of orthogonal coils. TF-SR measurements are performed to measure the temperature dependence of the magnetic-penetration depth (T). (T) is proportional to the superfluid density, and can provide information on the symmetry of the gap structure.

## Iii Results and Discussion

### a Normal and superconducting state properties

Figure 1 shows the room temperature powder X-ray diffraction pattern. Rietveld refinement of the RT XRD pattern shows that our sample is in the single phase, crystallizing with a cubic - structure (space group ) with lattice cell parameter a = 9.654(3) Å. A schematic view of the crystal structure of NbOs is shown in the inset of Fig. 1. Refined lattice parameters are shown in Table 1 and are in good agreement with the published literature BT (). In the crystal structure except for the Nb(1) site (2a), no other sites have an inversion center. The sample used here is found to have a refined stoichiometry of NbOs(3), close to the target composition NbOs, with the possibility of some site mixing between Nb and Os sites CT ().

Structure Cubic

Space group

Lattice parameters

a(Å) 9.654(3)

Atomic Coordinates | ||||
---|---|---|---|---|

Atom | Wyckoff position | x | y | z |

Nb1 | 2a | 0 | 0 | 0 |

Os1 | 8c | 0.322(1) | 0.322(1) | 0.322(1) |

Nb2 | 8c | 0.322(1) | 0.322(1) | 0.322(1) |

Os2 | 24g | 0.3536(6) | 0.3536(6) | 0.0371(8) |

Nb3 | 24g | 0.3536(6) | 0.3536(6) | 0.0371(8) |

Os3 | 24g | 0.0912(3) | 0.0912(3) | 0.2828(4) |

Temperature dependence of the electrical resistivity in the temperature range 1.8 K to 300 K in zero applied magnetic field is shown in Fig. 2 (a). The residual resistivity ratio (RRR) is 1.02(1). This low value for the RRR is comparable to other - structure compounds such as NbRe (1.3)ABK (), ReW (1.15) PKB (), ReHf (1.08) DSJ2 (), ReZr (1.09) DAM (), and ReTa (1.04) JAT1 (). These compounds are widely known to have strong electronic scattering, with a large residual resistivity due to occupational site disorder. This means that the resistivity of the materials with - structure is sensitive to disorder which is responsible for the poor conductivity. A superconducting transition can be seen clearly in the inset of Fig. 2(a) with a onset temperature = 2.7 0.2 K.

The normal-state resistivity for a non-magnetic metallic crystalline solid is analyzed using the Bloch-Grneisen (BG) model, which describes the resistivity arising due to electrons scattering from acoustic phonons. The temperature dependence of the resistivity, (T), is modeled as

(1) |

where is the residual resistivity due to the defect scattering and is essentially temperature independent whereas is the BG resistivity given by GG ()

(2) |

Here is the Debye temperature obtained from resistivity measurements, while C is a material dependent pre-factor and n takes values between 2-5 depending on the nature of the electron scattering ABA (). A fit to the data employing this model is shown in Fig. 2(a) and gave n = (4.2 0.2), Debye temperature = (376 2) K, C = (6.2 3.2) cm and residual resistivity = (130 2) cm. The value of the Debye temperature is close to that extracted from the specific heat data.

Figure 2(b) displays the temperature dependence of magnetic susceptibility measured in an applied field of 1 mT. Both zero-field cooled warming (ZFCW) and field-cooled cooling (FCC) regimes show a sharp diamagnetic transition around = 2.7 0.1 K, indicating the onset of bulk superconductivity. In the FCC data, the superconductor does not return to a full Meissner state, indicating strong flux pinning. The calculated Meissner volume fraction (4) is close to 100 , suggesting complete superconductivity within the sample.

To estimate the lower critical field , low-field magnetization data curves were measured at different fixed temperatures from 1.8 K to 2.5 K as shown in Fig. 3(a). The data follows a linear relation with the applied magnetic field below , whereas, above this point, magnetization starts to deviate from the straight line behavior due to flux penetration. The point of deviation was computed for different isotherm curves to obtain the temperature dependence of (T), as presented in the Fig. 3(a). The resulting graph is modeled using the Ginzburg-Landau relation

(3) |

When fitted to the experimental data, it yields (0) = 2.98 0.11 mT.

Figure 3(b) presents the high-field magnetization hysteresis loop collected in the superconducting state at T = 1.8 K, in the magnetic field range below 2 T. The magnetic behavior clearly corresponds to a type-II superconductor. The irreversibility field, , is defined as the field where the magnetic hysteresis collapses. At a given temperature for field , the magnetization is irreversible owing to pinned vortices, whereas for , it becomes reversible as the applied field is strong enough to depin the vortices. As shown in Fig. 3(b) at T = 1.8 K, = 0.38 0.02 T. The magnetization data were also collected in the normal state at T = 4 K and we conclude that there is no evidence for a magnetic impurity phase in our sample.

The magnetic hysteresis loop at different temperatures is displayed in Fig. 3(c). The hysteresis in magnetization M decreases with increasing temperature and magnetic field, characteristic of a conventional type-II superconductor. The gradient has a clear discontinuity at a field that is identified as the upper critical field [see inset Fig. 3(c)]. The resulting values of are determined using this discontinuity at different magnetic isotherms plotted in Fig. 3(d) (solid triangle markers). The temperature dependence of (T) is also determined from the ac susceptibility measurements (solid square markers Fig. 3(d)), where the shift in is evaluated under different applied magnetic fields up to 1.5 T [see inset of Fig. 3(d)]. The vs. curve can be described by the Werthamer-Helfand-Hohenberg (WHH) model EH (); NRW () by including the orbital breaking, the effect of Pauli spin paramagnetism () and spin-orbit scattering parameter (). The value of measures the relative strengths of the orbital and Pauli-limiting field, while is dominated by the atomic numbers of the elements. In the WHH model

(4) |

where t is the reduced temperature , is the diagamma function, , and is the dimensionless form of the upper critical field given by = (4.

The Maki parameter which measures the relative strengths of the orbital and Pauli-limiting field is calculated using the relation

(5) |

For a type-II superconductor in the dirty limit, the orbital limit of the upper critical field (0) is given by the WHH expression by

(6) |

For initial slope of T K near calculated from -T plot, it gives (0) = (3.36 0.07) T. The Pauli limiting field is given by (0) = 1.86 = (5.02 0.18) T. The Maki parameter is then calculated to be = 0.95.

Figure 3(d) show the temperature dependence of for two cases: = 0.95, = 0 and = 0, = 0, displayed by dotted black and blue lines. It is clear from the graph that the measured data is best described with = 0, whereas the calculation with = 0.95 fails to account for the experimental data. This certainly implies that is limited by orbital critical field and the Pauli limiting appears to have a limited effect if any at all.

The value for can therefore be calculated using the relation MAKI ()

(7) |

where it is considered = 0. For = 0, = (0) = (3.36 0.07) T.

The value for used to determine the Ginzburg Landau coherence length using the relation tin ()

(8) |

where (= 2.07 10 T m) is the magnetic flux quantum. For (0) = (3.36 0.07) T, we find = (99 1) Å.

Consequently, the values of and can be employed to calculate the Ginzburg-Landau penetration depth (0) from the expression tin ()

(9) |

Integrating the values of (0) = (2.98 0.11) mT and (0) = (99 1) Å in Eq. (9), yields (0) = (4608 88) Å. The Ginzburg-Landau parameter = / = 46.5 0.4 1/, therefore, NbOs classified as strong type-II superconductor.

The thermodynamic critical field (0) is given by

(10) |

and found to be (0) = (51 1) mT.

Recently unconventional vortex dynamics has been observed in some noncentrosymmetric superconductors where unusual vortex pinning mechanism and flux creep rates were proposed CFM1 (); CFM2 (). Therefore, it is necessary to measure the stability of vortex system against the various factors unsettling the vortex equilibrium. The measure of stability against the thermal fluctuations can be given by Ginzburg number . This is basically the ratio of thermal energy to the condensation energy associated with the coherence volume. Ginzburg number, , is given by GB ()

(11) |

where is anisotropy parameter which is 1 for the cubic NbOs. For (0) = (99 1) Å, (0) = (51 1) mT and = (2.7 0.1) K, we obtained = (1.1 0.1) 10. In conventional low superconductors, pinning is strong, whereas thermal fluctuations are weak, with 10. In high temperature superconductors, is high and hence the coherence volume is small, making it sensitive to thermal fluctuations, 10. In our compound value is more towards the low superconductors, suggesting that thermal fluctuations may not be playing any important role in vortex unpinning in this material.

Figure 4 shows the low temperature specific heat data measured in zero applied field. The plot vs is shown in the inset of Fig. 4, in the temperature range 3 K 100 K. A jump at around = 2.7 K, apparently confirms bulk superconductivity in NbOs. The normal state specific heat data above is contributed by both electronic and phononic parts given by: , where is the Sommerfeld coefficient and , are phononic contributions. The solid red line represents the best fit to the data with = 8.58 0.01 mJ mol K, = 0.123 0.002 mJ mol K, and = 0.43 0.04 J mol K [see inset Fig. 4]. The value for was used to determine the density of states at the Fermi level using the relation , where is the Fermi energy. For = 8.58 0.01 mJ mol K, it yields = 3.64 0.02 states eV f.u.. The Debye temperature = , where using R = 8.314 J mol K and N = 3, yields = 362 1 K. The value of = 362 K can be used to calculate the electron-phonon coupling constant using the McMillan formula WL (),

(12) |

where represents the repulsive screened Coulomb potential, usually given by = 0.13. With = 2.7 K and = 361.67 K, we obtained 0.52, which is similar to other weakly coupled NCS superconductors ABK (); SAS ().

The electronic contribution () to the specific heat can be calculated by subtracting the phononic contribution () from total specific heat data.

(13) |

Once the phononic contribution is subtracted, an equal entropy conservation line is drawn to estimate the normalized specific heat jump. The value for the specific heat jump 1.45, which is remarkebly similar to the value for a BCS isotropic gap superconductor (= 1.43) in the weak coupling limit, indicating that NbOs is a weakly coupled superconductor, consistent with the value obtained above.

The temperature dependence of the specific heat data in the superconducting state can best be described by the BCS expression for the normalized entropy S written as

(14) |

where f() = [exp(E()/)+1] is the Fermi function, E() = , where is the energy of normal electrons measured relative to the Fermi energy, y = , , and = tanh[1.82(1.018(()-1))] is the BCS approximation for the temperature dependence of the energy gap. The normalized electronic specific heat is calculated by

(15) |

Fitting the low temperature specific heat data using this model as shown by the solid red line in Fig. 4, yields = = 1.77 0.01. This is consistent with the value for a BCS superconductor = 1.764 in the weak coupling limit. Therefore, good agreement between the measured data (black symbols) and the BCS fit (solid red line), confirms an isotropic fully gapped BCS superconductivity in NbOs. It should be noted that in order to extract the true nature of the superconducting gap, the specific heat data up to the low temperature region need to be analyzed. It can give a complete understanding whether the superconducting gap is isotropic (exponential) or have nodes (power law). A summary of all the experimentally measured and estimated parameters is given in Table 2.

To further inspect the superconducting ground state of NbOs, we performed the muon spin relaxation and rotation measurements. Firstly, ZF-SR measurements were carried out to investigate the occurrence of TRS breaking in NbOs. ZF-SR being extremely sensitive to tiny changes in the internal magnetic fields, can unambiguously detect the presence of a TRS breaking signal. This technique was proved to be very crucial in establishing the TRS breaking in several superconductors such as LaIr JAT (), Re(Ti,Hf,Zr) RPS (); DSJ (); DSJ1 (), SrRuO SRO1 (); SRO2 () etc. Figure 5 shows the muon spin relaxation time spectra, collected below ( K < ) and above ( K > ) the superconducting transition temperature = 2.7 K. The absence of precessional signals suggests the absence of coherent internal fields which is generally associated with long range magnetic ordering. In the absence of atomic moments muon-spin relaxation in zero field can best be described by the Gaussian Kubo-Toyabe (KT) function RSH () together with a non-decaying constant background, ,

(16) | |||||

where is the initial asymmetry and accounts for the relaxation due to static, randomly oriented local fields associated with the nuclear moments at the muon site and is the electronic relaxation rate. Near identical relaxation signals, as seen in Figure 5, on the either side of the superconducting transition suggest the absence of any additional magnetic moments in the superconducting ground state, usually associated with exotic phenomena such as TRS breaking. This clearly suggests that the time-reversal symmetry is preserved in NbOs within the detection limit of SR.

Figure 6(a) shows the TF-SR precessional signals for NbOs collected above (3.5 K) and below (0.4 K) in an applied magnetic field of 20 mT. The measurements were done in the field-cooled mode where the field 20 mT was applied above the transition temperature. The sample was then subsequently cooled to the base temperature of 0.4 K.
In the normal state there is an almost homogeneous field distribution throughout the sample, where the weak depolarization is due to nuclear dipolar field. The depolarization rate in the superconducting state becomes more prominent, due to the formation of an inhomogeneous field distribution in the flux line lattice (FLL) state. The TF-SR spectra in 6(a) show a very small difference between above and below , indicating large magnetic penetration depth.

The TF-SR time spectra, for all temperatures, is best described by the sinusoidal oscillatory function damped with a Gaussian relaxation and an oscillatory background term arising from the muons implanted directly into the silver sample holder that do not depolarize,

(17) | |||||

Here and are the initial asymmetries of the sample and background signals, and are the internal and background magnetic fields, is the initial phase offset, = 135.53 MHz/T is the muon gyromagnetic ratio and is the depolarization rate. The depolarization rate comprised of the following terms: = , where is the depolarization arising due to the field variation across the flux line lattice and is the contribution due to nuclear dipolar moments. The superconducting contribution to the depolarization rate is calculated by using the above relation. Figure 6(b) shows the temperature dependence of where below , increases systematically as the temperature is lowered. The inset of Fig. 6(b) shows the temperature dependence of internal magnetic field, where the flux expulsion at is evident from the reduction of average feld inside the superconductor, and the corresponding background field is approximately constant over the temperature range. In an isotropic type-II superconductor with a hexagonal Abrikosov vortex lattice the magnetic penetration depth is related to by the equation EHB ():

(18) |

where is the reduced field. The above equation is valid for systems . For NbOs, = (46.5 0.4) and the temperature dependence of is shown in Fig. 3(d). Using the data of ) for NbOs, the temperature dependence of was extracted from Eq. (18), as displayed in Fig. 6(c). The temperature dependence of nearly constant below /3 0.9 K. This possibly suggests the absence of nodes in the superconducting energy gap at the Fermi surface. The solid line in Fig. 6(c) represent the temperature dependence of the London magnetic penetration depth within the local London approximation for a s-wave BCS superconductor in the dirty limit:

(19) |

where is the BCS approximation for the temperature dependence of the energy gap, where is the gap magnitude at zero temperature. Dirty-limit expression was used in accordance with the calculation done in the section below where we found > l. The fit yields energy gap = (0.43 0.02) meV and (0) = (6190 45) Å. The gap to ratio = 1.82 is close to the value of 1.764 expected from the BCS theory, suggesting that NbOs is a weakly-coupled superconductor.

Properties | unit | value |
---|---|---|

K | 2.7 0.2 | |

mT | 2.98 0.11 | |

T | 3.36 0.07 | |

mT | 51 1 | |

T | 5.02 0.18 | |

Å | 99 1 | |

Å | 4608 88 | |

46.5 0.4 | ||

mJmolK | 8.58 0.01 | |

mJmolK | 0.123 0.002 | |

K | 362 1 | |

0.52 0.02 | ||

D(E) | states/ev f.u | 3.64 0.02 |

1.45 0.03 | ||

1.77 0.01 |

## Iv Electronic properties

The Sommerfeld coefficient is related to the quasiparticle number density (n) per unit volume given by the relation ck ()

(20) |

where k is the Boltzmann constant, N is the Avogadro number, V is the volume of a formula unit and m is the effective mass of quasiparticles. The residual resistivity, , can be calculated, using the equation

(21) |

while the Fermi velocity can be written in terms of effective mass and the carrier density by

(22) |

For superconductors in the dirty-limit, where 1, the properties are affected due to the scattering of electrons. The dirty limit expression for the penetration depth (0) is then given by tin ()

(23) |

where is the BCS coherence length. The London penetration depth, , is given by

(24) |

The BCS coherence length and the Ginzburg-Landau coherence, (0), at T = 0 K in the dirty limit is related by the expression

(25) |

Using Eq.(20)-(25) form a system of equations which can be used to estimate the parameters m, n, l, and as done in Ref.DAM (); JAT1 (). The system of equations was solved simultaneously using the values = (8.58 0.01) mJ molK, = (99 1) Å and = (130 5) cm. The values were estimated for 4608 Å and 6190 Å, tabulated in Table 3. It is clear that > l, indicating that NbOs is in the dirty limit as previously asserted. This also accounts for the low RRR and residual resistivity that have been measured.

The bare-band effective mass m can be related to m, which contains enhancements from the many-body electron phonon interactions GG () m = m(1+), where is the electron-phonon coupling constant. Using this value of = 0.52, a value for m can be found, given in Table 3.

Properties | unit | ||
---|---|---|---|

Å | 4608 | 6190 | |

11.6 0.3 | 13.4 0.3 | ||

7.63 0.13 | 8.8 0.2 | ||

10m | 4.16 0.23 | 2.74 0.17 | |

Å | 64 5 | 90 6 | |

Å | 37 4 | 50 3 | |

Å | 2809 130 | 3720 164 | |

10m s | 4.95 0.23 | 3.75 0.16 | |

K | 1060 30 | 620 20 | |

0.0025 0.0001 | 0.0044 0.0001 |

## V Uemura Plot

For a 3D system the Fermi temperature T is given by the relation

(26) |

where n is the quasiparticle number density per unit volume. According to Uemura et al. YJU1 (); YJU2 (); YJU3 () superconductors can be conveniently classified according to their ratio. It was shown that for the unconventional superconductors such as heavy-fermion, high- , organic superconductors, and iron-based superconductors this ratio falls in the range 0.01 0.1. In Fig. 7, the orange region represents the band of unconventional superconductors. Uemura plot is presented in accordance with Ref. JAT1 (); YJU1 (); YJU2 (); YJU3 (), where it shows the superconductors with unconventional properties such as MoAlC, MgIrB LiPtB, LaNiC, LaIr and the family of Re-based - structure superconductors with broken TRS. Using the estimated value of n and m for NbOs in Eq. 26, it yields = 620 K, giving = 0.004. NbOs is located just outside the range of unconventional superconductors as shown by a solid red marker in Fig. 7, potentially borders an unconvnetional superconducting ground state. Interestingly, all the ReX RPS (); DSJ (); DSJ1 () noncentrosymmetric superconductors with TRS breaking are located in the same phase space in the Uemura plot. The close proximity of the entire family might also be pointing towards the common origin of TRS breaking, which seems to be Re SOC. However, other Re-based compounds for example ReX (X =W, Ta)PKB (); JAT1 () show preserved TRS. This may be due to the reduced Re composition in these compounds which may have a major effect on the underlying electronic characteristics of the spin-triplet channel. In addition, our ZF-SR measurements on the Re-free NbOs also advocate the above conclusion, since TRS is found to be preserved in the superconducting state of this material. Therefore, this makes the role of Re element even more interesting. However, to make a generic comment about the contribution of Re in the observed TRS breaking in this family, it is important to study more Re-rich - materials and more importantly Re itself.

## Vi Summary and Conclusions

Detailed investigation of the normal and superconducting phase properties of NbOs was carried out using XRD, magnetization, resistivity, specific heat and SR measurements. Sample was prepared by standard arc-melting technique, where the phase purity and noncentrosymmetric - structure for NbOs was confirmed by XRD analysis. Our results suggest type-II dirty limit superconductivity in NbOs with superconducting transition temperature = 2.7 0.1 K. The specific heat measurements confirm the superconducting gap is isotropic and are, within error, the same as the BCS predicted values. ZF-SR measurement confirmed that time-reversal symmetry is preserved, within the detection limit of SR. The TF-SR measurements also confirm fully gapped BCS superconductivity with no point or line nodes. In addition, in the Uemura plot NbOs is placed just outside the borderline of the band of unconventionality. Re-based - structure superconductors with broken TRS are situated very close to each other in the Uemura plot. The negligible effect of ASOC and persistence occurrence of TRS breaking in this family of compounds hints at its common origin or mechanism. It could be due to the Re SOC which facilitates the origin of TRS breaking in this family. The above conclusion is supported by the absence of TRS breaking phenomena in Re-free NbOs having the same structure. Further confirmation is provided by the recent observation of TRS breaking in pure Re metal (centrosymmetric). Therefore, it is imperative and timely to search for TRS breaking in other Re-rich superconductors and employ both experimental and theoretical approach to identify the origin of TRS breaking.

## Vii Acknowledgments

R. P. S. acknowledges Science and Engineering Research Board, Government of India for the Young Scientist Grant No. YSS/2015/001799 and Department of Science and Technology, India (SR/NM/Z-07/2015) for the financial support and Jawaharlal Nehru Centre for Advanced Scientific Research (JNCASR) for managing the project. We thank ISIS, STFC, UK for the beamtime to conduct the SR experiments

## References

- (1) E. Bauer, G. Hilscher, H. Michor, Ch. Paul, E. W. Scheidt, A. Gribanov, Yu. Seropegin, H. NoÂšel, M. Sigrist, and P. Rogl, Phys. Rev. Lett. 92, 027003 (2004).
- (2) E. Bauer and M. Sigrist, Non-centrosymmetric superconductors:Introduction and Overview, Vol. 847 (Springer Science Business Media, 2012).
- (3) L. P. Gor’kov and E. I. Rashba, Phys. Rev. Lett. 87, 037004 (2001).
- (4) P. A. Frigeri, D. F. Agterberg, A. Koga, and M. Sigrist, Phys. Rev. Lett. 92, 097001 (2004).
- (5) E. Bauer, G. Hilscher, H. Michor, C. Paul, E. W. Scheidt, A. Gribanov, Y. Seropegin, H. NoÃ«l, M. Sigrist, and P. Rogl, Phys. Rev. Lett. 92, 027003 (2004).
- (6) N. Kimura, K. Ito, H. Aoki, S. Uji, and T. Terashima, Phys. Rev. Lett. 98, 197001 (2007).
- (7) I. Sugitani, Y. Okuda, H. Shishido, T. Yamada, A. Thamizhavel, E. Yamamoto, T. D.Matsuda, Y. Haga, T. Takeuchi, R. Settai, and Y. Onuki, J. Phys. Soc. Jpn. 75, 043703 (2006)
- (8) I. Bonalde, W. Brmer-Escamilla, and E. Bauer Phys. Rev. Lett. 94, 207002 (2005).
- (9) E. Bauer, G. Rogl, X.-Q. Chen, R. T. Khan, H. Michor, G. Hilscher, E. Royanian, K. Kumagai, D. Z. Li, Y. Y. Li, R. Podloucky, and P. Rogl, Phys. Rev. B 82, 064511 (2010).
- (10) V. M. Edelstein, Phys. Rev. Lett. 75, 2004 (1995).
- (11) S. Fujimoto, Phys. Rev. B 72, 024515 (2005).
- (12) G. Bian, T.R. Chang, R. Sankar, S.Y. Xu, H. Zheng, T. Neupert, C.K. Chiu, S.M. Huang, G. Q. Chang, I. Belopolski, D. S. Sanchez, M. Neupane, N. Alidoust, C. Liu, B. K. Wang, C.C Lee, H.T. Jeng, C. L. Zhang, Z. J. Yuan, S. Jia, A. Bansil, F. C. Chou, H. Lin, and M. Z. Hasan, Nat. Commun. 7, 10556 (2016).
- (13) S. Y. Guan, P. J. Chen, M. W. Chu, R. Shankar, F. Chou, H. T. Jeng, C. S. Chang, and T. M. Chuang, Sci. Adv. 2, e1600894 (2016).
- (14) Z. Sun, M. Enayat, A. Maldonado, C. Lithgow, E. Yelland, D. C. Peets, A. Yaresko, A. P. Schnyder, and P. Wahl, Nat. Commun. 6, 6633 (2015).
- (15) M. Neupane, N. Alidoust, M. M. Hosen, J.-X. Zhu, K. Dimitri, S.-Y. Xu, N. Dhakal, R. Sankar, I. Belopolski, D. S. Sanchez, T.-R. Chang, H.-T. Jeng, K. Miyamoto, T. Okuda, H. Lin, A. Bansil, D. Kaczorowski, F. C. Chou, M. Z. Hasan, and T. Durakiewicz, Nat. Commun. 7, 13315 (2016).
- (16) A. D. Hillier, J. Quintanilla, and R. Cywinski, Phys. Rev. Lett. 102, 117007 (2009).
- (17) J. Quintanilla, A. D. Hillier, J. F. Annett, and R. Cywinski, Phys. Rev. B 82, 174511 (2010).
- (18) K. Togano, P. Badica, Y. Nakamori, S. Orimo, H. Takeya, and K. Hirata, Phys. Rev. Lett.93, 247004 (2004).
- (19) P. Badica, T. Kondo, T. Kudo, Y. Nakamori, S. Orimo, and K. Togano, Appl. Phys. Lett. 85, 4433 (2004).
- (20) P. Badica, T. Kondo, and K. Togano, J. Phys. Soc. Jpn.74, 1014 (2005).
- (21) M. Nishiyama, Y. Inada, and Guo-qing Zheng, Phys. Rev. B 71, 220505(R)(2005).
- (22) H.Takeya, K.Hirata, K.Yamaura, K.Togano,M.El Massalami, R. Rapp, F. A. Chaves, and B. Ouladdiaf, Phys. Rev. B 72, 104506 (2005).
- (23) H. Q. Yuan, D. F. Agterberg, N. Hayashi, P. Badica, D. Vandervelde, K. Togano, M. Sigrist, and M. B. Salamon, Phys. Rev. Lett. 97, 017006 (2006).
- (24) M. Nishiyama, Y. Inada, and Guo-qing Zheng, Phys. Rev. Lett. 98, 047002 (2007).
- (25) H. Takeya, M. ElMassalami, S. Kasahara, and K. Hirata, Phys. Rev. B 76, 104506 (2007).
- (26) L. Fang, H. Yang, X. Zhu, G. Mu, Z.S. Wang, L. Shan, C. Ren, and H.H. Wen, Phys. Rev. B79, 144509 (2009).
- (27) A. B. Karki, Y. M. Xiong, N. Haldolaarachchige, S. Stadler, I. Vekhter, P. W. Adams, D. P. Young, W. A. Phelan, and J. Y. Chan, Phys. Rev. B 83, 144525 (2011).
- (28) P. K. Biswas, M. R. Lees, A. D. Hillier, R. I. Smith, W. G.Marshall, and D. M. Paul, Phys. Rev. B 84, 184529 (2011).
- (29) J. Chen, M.B. Salamon, S. Akutagawa, J. Akimitsu, J. Singleton, J. L. Zhang, L. Jiao, and H. Q. Yuan, Phys. Rev. B 83, 144529 (2011).
- (30) E. Bauer, G. Rogl, X.-Q. Chen, R. T. Khan, H. Michor, G. Hilscher, E. Royanian, K. Kumagai, D. Z. Li, Y. Y. Li, R. Podloucky, and P. Rogl, Phys. Rev. B82, 064511 (2010).
- (31) A. B. Karki, Y. M. Xiong, I. Vekhter, D. Browne, P. W. Adams, D.P.Young,K.R.Thomas,J.Y.Chan, H.Kim, and R.Prozorov, Phys. Rev. B 82, 064512 (2010).
- (32) T. Mochiku, H. Takeya, T. Wuernisha, K. Mori, T. Ishigaki, T. Kamiyama, H. Fujii, and K. Hirata, Physica C 445, 57 (2006).
- (33) R. Khasanov, I. L. Landau, C. Baines, F. La Mattina, A. Maisuradze, K.Togano,and H.Keller, Phys. Rev. B 73,214528 (2006).
- (34) M. M. Doria, S. Salem-Sugui, Jr., P. Badica, and K. Togano, Phys. Rev. B 73, 184524 (2006).
- (35) H. Q. Yuan, D. Vandervelde, M. B. Salamon, P. Badica, and K. Togano, AIP Conf. Proc. 850, 633 (2006).
- (36) S. Tsuda, T. Yokoya, T. Kiss, T. Shimojima, K. Ishizaka, S. Shin, T. Togashi, S. Watanabe, C. Zhang, C. Chen, I. Hase, H. Takeya, K. Hirata, and K. Togano, J. Phys. Soc. Jpn. 78, 034711 (2009).
- (37) A. Shimamura, Y. Furukawa, K. Kumagai, H. Takeya, and K. Hirata, Physica C 460, 663 (2007).
- (38) D. C. Peets, G. Eguchi, M. Kriener, S. Harada, Sk. Md. Shamsuzzamen, Y. Inada, G.-Q. Zheng, and Y. Maeno, Phys. Rev. B 84, 054521 (2011).
- (39) G. Eguchi, D. C. Peets, M. Kriener, S. Yonezawa, G. Bao, S. Harada, Y. Inada, G.-q. Zheng, and Y. Maeno, Phys. Rev. B 87, 161203(R)(2013).
- (40) M. Smidman, A. D. Hillier, D. T. Adroja, M. R. Lees, V. K. Anand, R. P. Singh, R. I. Smith, D. M. Paul, and G. Balakrishnan Phys. Rev. B 89, 094509 (2014).
- (41) V. K. Anand, A. D. Hillier, D. T. Adroja, A. M. Strydom, H. Michor, K. A. McEwen, and B. D. Rainford Phys. Rev. B 83, 064522 (2011).
- (42) V. K. Anand, D. Britz, A. Bhattacharyya, D. T. Adroja, A. D. Hillier, A. M. Strydom, W. Kockelmann, B. D. Rainford, and K. A. McEwen Phys. Rev. B 90, 014513 (2014).
- (43) K. Matano, S. Maeda, H. Sawaoka, Y. Muro, T. Takabatake, B. Joshi, S. Ramakrishnan, K. Kawashima, J. Akimitsu, and G.-Q. Zheng,J. Phys. Soc. Jpn. 82, 084711 (2013).
- (44) X. B. Yan, Y. Xu, L. P. He, J. K. Dong, H. Cho, D. C. Peets, J.-G. Park, and S. Y. Li, Supercond. Sci. Technol. 29, 065001 (2016).
- (45) J. A. T. Barker, D. Singh, A. Thamizhavel, A. D. Hillier, M. R. Lees, G. Balakrishnan, D McK. Paul, and R. P. Singh, Phys. Rev. Lett. 115, 267001 (2015)
- (46) P. K. Biswas, H. Luetkens, T. Neupert, T. StÃŒrzer, C. Baines, G. Pascua, A. P. Schnyder, M. H. Fischer, J. Goryo, M. R. Lees, H. Maeter, F. BrÃŒckner, H. H. Klauss, M. Nicklas, P. J. Baker, A. D. Hillier, M. Sigrist, A. Amato, and D. Johrendt, Phys. Rev. B 87, 180503 (2013).
- (47) R. P. Singh, A. D. Hillier, B. Mazidian, J. Quintanilla, J. F. Annett, D. M. Paul, G. Balakrishnan, and M. R. Lees, Phys. Rev. Lett. 112, 107002 (2014).
- (48) D. Singh, J. A. T. Barker, A. Thamizhavel, D. McK. Paul, A. D. Hillier, and R. P. Singh, Phys. Rev. B 96, 180501(R) (2017).
- (49) D. Singh, Sajilesh K. P., J. A. T. Barker, D. McK. Paul, A. D. Hillier, and R. P. Singh Phys. Rev. B 97, 100505(R) (2018).
- (50) M.J. Winiarski, J Alloys Compd. 616, 1 (2014).
- (51) Mojammel A. K., A. B. Karki, T. Samanta, D. Browne, S. Stadler, I.Vekhter, Abhishek Pandey, P.W.Adams, D.P.Young, S.Teknowijoyo, K. Cho, R. Prozorov, and D. E. Graf, Phys. Rev. B 94, 144515 (2016).
- (52) T. Shang, M. Smidman, S. K. Ghosh, C. Baines, L. J. Chang, D. J. Gawryluk, J. A. T. Barker, R. P. Singh, D. Mck. Paul, G .Balakrishnan, E. Pomjakushina, M. Shi, M. Medarde, A. D. Hillier, H. Q. Yuan, J. Quintanilla, J. Mesot, T. Shiroka, Phys. Rev. Lett. 121, 257002 (2018).
- (53) B. T. Matthias, V. B. Compton, and E. Corenzwit, J. Phys. Chem. Solids 19, 180 (1961).
- (54) B. T. Matthias, T. H. Geballe, and V. B. Compton, Rev. Mod. Phys. 35 (1963).
- (55) S. L. Lee, S. H. Kilcoyne, and R. Cywinski, Muon Science: Muons in Physics, Chemistry and Materials (SUSSP Publications and IOP Publishing, Bristol, 1999).
- (56) R. M. Waterstrat and R. C. Manuszewski Journal of the Less Common Metals, 51, 55-67(1977).
- (57) D. Singh, A. D. Hillier, A. Thamizhavel, and R. P. Singh, Phys. Rev. B 94, 054515 (2016).
- (58) D. A. Mayoh, J. A. T. Barker, R. P. Singh, G. Balakrishnan, D. McK. Paul, and M.R.Lees, Phys. Rev. B 96 ,064521 (2017).
- (59) J. A. T. Barker, B. D. Breen, R. Hanson, A. D. Hillier, M. R. Lees, G. Balakrishnan, D. McK. Paul, and R. P. Singh Phys. Rev. B 98, 104506 (2018).
- (60) G. Grimvall, The Electron-Phonon Interaction in Metals(North Holland, Amsterdam, 1981).
- (61) A. Bid, A. Bora, and A. K. Raychaudhuri, Phys. Rev. B74, 035426 (2006).
- (62) E. Helfand and N. R. Werthamer, Phys. Rev. 147, 288 (1966).
- (63) N. R. Werthamer, E. Helfand, and P. C. Hohenberg, Phys. Rev. 147, 295 (1966).
- (64) K. Maki, Phys. Rev. 148, 362 (1966).
- (65) M. Tinkham, Introduction to Superconductivity, 2nd ed. (McGraw-Hill, New York, 1996).
- (66) C. F. Miclea, A. C. Mota, M. Nicklas, R. Cardoso, F. Steglich, M. Sigrist, A. Prokofiev, and E. Bauer, Phys. Rev. B 81, 014527 (2010).
- (67) C. F. Miclea, A. C. Mota, M. Sigrist, F. Steglich, T. A. Sayles, B. J. Taylor, C. A. McElroy, and M. B. Maple, Phys. Rev. B 80, 132502 (2009).
- (68) G. Blatter, M. V. Feigel’man, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur Rev. Mod. Phys. 66, 1125 (1994).
- (69) W. L. McMillan, Phys. Rev. 167, 331 (1968).
- (70) M. Isobe, A.Masao, and N. Shirakawa, Phys. Rev. B 93, 054519 (2016).
- (71) G. M. Luke, Y. Fudamoto, K. M. Kojima, M. I. Larkin, J. Merrin, B. Nachumi, Y. J. Uemura, Y. Maeno, Z. Q. Mao, Y. Mori, H. Nakamura, and M. Sigrist, Nature (London) 394, 558 (1998).
- (72) J. Xia, Y. Maeno, P. T. Beyersdorf, M. M. Fejer, and A. Kapitulnik, Phys. Rev. Lett.97, 167002 (2006).
- (73) R. S. Hayano, Y. J. Uemura, J. Imazato, N. Nishida, T.Yamazaki, and R. Kubo, Phys. Rev. B 20, 850 (1979).
- (74) E. H. Brandt, Phys. Rev. B 68, 054506 (2003).
- (75) C.Kittel, Introduction to Solid State Physics 8th edn. (Wiley, New York, 2005).
- (76) Y. J. Uemura et al., Phys. Rev. B. 38, 909 (1988).
- (77) Y. J. Uemura et al., Phys. Rev. Lett. 62, 2317 (1989).
- (78) Y. J. Uemura, L. P. Le, G. M. Luke, B. J. Sternlieb, W. D. Wu, J. H. Brewer, T. M. Riseman, C. L. Seaman, M. B. Maple, M. Ishikawa, D. G. Hinks, J. D. Jorgensen, G. Saito, and H. Yamochi, Phys. Rev. Lett.66, 2665 (1991).