Superconducting Gap Structure of the Noncentrosymmetric Topological Superconductor Candidate HfRuP

: We investigate the gap symmetry of the topological superconductor candidate HfRuP, which crystallizes in a noncentrosymmetric hexagonal crystal structure, using muon spin rotation/relaxation ( µ SR) measurements in transverse-ﬁeld (TF) geometry. The temperature and magnetic ﬁeld dependencies of the superconducting relaxation rate derived from the TF-µ SR spectra can be well described by an isotropic s -wave gap. The superconducting carrier density n s = 1.41(1) × 10 26 m − 3 and the magnetic penetration depth, λ (0) = 603(2) nm, were calculated from the TF-µ SR data. Inter-estingly, the ratio between the superconducting transition temperature and the superﬂuid density, T c / λ − 2 (0) ∼ 3.3, is very close to those of unconventional superconductors. Further, our zero-ﬁeld (ZF) µ SR results do not show any considerable change in the muon spin relaxation above and below the superconducting transition temperature, suggesting that time-reversal symmetry is preserved in the superconducting state of this superconductor.

HfRuP crystallizes in a hexagonal Fe 2 P-type structure without any inversion symmetry in the crystal structure [9,11,17].Recent ARPES studies on this compound show that HfRuP has 12 pairs of type-II Weyl points [17].Thus, it is very important to understand the superconducting gap structure of this interesting superconductor.Further, the noncentrosymmetric structure of this compound, containing heavy elements with strong spin-orbit interactions, sparks the possibility of mixed spin-singlet and spin-triplet pairing [18][19][20], which breaks the time-reversal symmetry.Tempted by these intriguing aspects, we carried out detailed muon spin rotation/relaxation (µSR) measurements on HfRuP.
In the case of a type-II superconductor, a vortex state creates an inhomogeneous spatial distribution of local magnetic fields.Such a field distribution influences the muon spin depolarization rate in the superconducting state.From the superconducting state depolarization rate, we can directly calculate the magnetic penetration depth λ.The temperature dependence of λ is related to the superconducting gap structure.Further, zero-field µSR is a well-established experimental technique to check if the time-reversal symmetry is broken in the superconducting state [21].Here, we demonstrate that HfRuP contains a single s-wave pairing gap and the time-reversal symmetry is preserved in the superconducting sate of this compound.

Experimental Details
A polycrystalline sample of HfRuP was prepared by pulverizing single crystals grown using the Cu-P flux method.To grow the single crystals, high-purity Hf, Ru, P, and Cu elements were combined in a 1:1:5.6:16.8molar ratio in an alumina crucible within a glove box.The crucible was then sealed inside a tantalum ampule under an argon atmosphere and subsequently sealed in an evacuated quartz tube.The quartz tube was heated to 1170 • C for 15 h and then left for 5 h.Then, the tube was slowly cooled down to 1100 • C and immediately quenched into ice water to form HfRuP with a hexagonal structure.Finally, needle-like single crystals were yielded after dissolving the Cu-As fluxes in nitric acid.Primary characterization data have been published elsewhere [17].Transverse-field (TF)and zero-field (ZF)-µSR experiments were carried out at the Paul Scherrer Institute (Villigen, Switzerland).The temperature-dependent TF and ZF experiments were performed on the DOLLY spectrometer and field-dependent TF experiments at 1.6 K were carried out on a GPS spectrometer.We first ground the sample of HfRuP to a powder, which was in the form of small single crystals, and pressed it into a 9 mm diameter pellet, using a hydraulic press, that was then mounted on a Cu holder using GE varnish.It was then placed in the appropriate spectrometer cryostat.DOLLY and GPS spectrometers are equipped with a standard veto setup [22] which allows a low-background µSR signal.For the TF experiments, a field was first applied above T c , followed by cooling the sample.The µSR time spectra were analyzed using the MUSRFIT software package [23].

Crystal Structure and Magnetization
The Rietveld refinement of the XRD pattern of HfRuP presented in Figure 1a reveals the single phase nature of the polycrystalline sample with an Fe 2 P-type hexagonal structure with space group P 62m (No. 189, Z = 3).The obtained lattice parameters and atomic position parameters are in good agreement with the values reported previously in the literature [9,24].Figure 1b,c show the layered hexagonal structure of HfRuP.Each layer in the hexagonal lattice is occupied by either Hf and P atoms or Ru and P atoms.The triangular clusters of three Ru atoms are formed in the ab-plane, while the corner-sharing triangular lattice is formed by Hf atoms.Furthermore, the Hf ions are arranged in a geometrically frustrated quasikagome network in the hexagonal basal plane.The Fe 2 P-type hexagonal structure does not have a center of inversion and hence the crystal structure of HfRuP is noncentrosymmetric.The structure has two kinds of mirror symmetries, m z and m x , which are vital to define the mirror Chern numbers.Figure 1d shows the temperature dependence of the zero-field-cooled (ZFC) and field-cooled (FC) dc magnetic susceptibility in an applied magnetic field of 20 Oe.A clear diamagnetic signal is observed in both the ZFC and FC curves, with an onset superconducting critical temperature T C = 9 K. Further, the values of ZFC magnetic susceptibility below T C reveal a 100% superconducting volume fraction in HfRuP.

TF-µSR Measurements
In Figure 2a, we present TF-µSR asymmetry spectra for HfRuP at temperatures above (20 K) and below (0.27 K) T c .Temperature-dependent TF-µSR measurements were performed in an applied magnetic field of 30 mT, which is higher than the lower critical field, H c1 (0) = 6.05 mT, and lower than the upper critical field, H c2 (0) = 13.4 T. The TF-µSR spectra above T c show a very small relaxation due to the presence of random local fields associated with the nuclear magnetic moments.For a type-II superconductor, in the mixed state, due to the formation of the flux line lattice (FLL), the relaxation rate of the µSR signal is enhanced.Thus, we observe a rapidly relaxing signal below T c for HfRuP.Assuming a Gaussian field distribution, the observed TF-µSR asymmetry spectra can be described using where A 0 refers to the initial asymmetry, γ µ /(2π) 135.5 MHz/T is the muon gyromagnetic ratio, ϕ is the initial phase of the muon spin ensemble, B int corresponds to the internal magnetic field at the muon site, and σ is the total relaxation rate.σ is related to the superconducting relaxation rate, σ SC , via the relation σ = σ 2 nm + σ 2 SC , where σ nm is the nuclear contribution which is assumed to be temperature independent.We considered the value of σ nm obtained from fitting TF spectra above T c and kept it fixed.The fits to the observed spectra with Equation (1) are shown in solid lines in Figure 2a. Figure 2b shows the Fourier transform amplitudes of the TF-µSR asymmetry spectra recorded at 20 K and 0.27 K (Figure 2a).We observed a sharp peak in the Fourier amplitude around 30 mT (external applied field) at 20 K, confirming homogeneous field distribution throughout the sample.In contrast, a fairly broad signal, with a peak position slightly shifted to a lower value (diamagnetic shift), was observed at 0.27 K, confirming that the sample is in the superconducting mixed state.The formation of FLL results in such broadening of the line shape.Considering a perfect triangular vortex lattice, the muon spin depolarization rate σ sc (T) can be related to the London magnetic penetration depth λ(T) by [25,26]: where Φ 0 = 2.068 × 10 −15 Wb is the magnetic flux quantum.Equation ( 2) is only valid when the separation between the vortices is smaller than λ, which is presumed to be field independent in this model [25].Figure 3a shows the temperature dependence of λ −2 (T) derived from σ sc (T) using Equation (2).The inset of Figure 3a shows the temperature evolution of the relative change in the internal field (i.e., change in field normalized to the applied field), . As expected for type-II superconductors, we observe that below T c the internal field values in the superconducting state are lower than the applied field due to a diamagnetic shift.
To comprehend the superconducting gap structure and, further, to extract quantitative information of different superconducting parameters of HfRuP, we analyzed the temperature dependence of the magnetic penetration depth, λ(T), which is directly associated with the superconducting gap.Within the London approximation (λ ξ), λ −2 (T) can be correlated with the superconducting gap of an s-wave superconductor using: [23,27,28] here, f = [1 + exp(E/k B T)] −1 is the Fermi function and ∆ i (T) = ∆ 0,i Γ(T/T c ). ∆ 0,i is the maximum gap value at T = 0 K.The temperature dependence of the gap is expressed as Γ(T/T c ) = tanh 1.82[1.018(Tc /T − 1)] 0.51 [29].As seen from Figure 3a, the experimentally obtained λ −2 (T) can be well described by a momentum-independent single-gap s-wave model (red solid line) with a gap value of ∆ 0 = 1.34(1) meV and T c = 9.22(6) K.We note that the T c of our sample is very similar to that reported for the single crystal of HfRuP (T c = 9 K) by Qian et al. [17], however, lower than the T c = 12 K reported by Meisner et al. [9,10].Interestingly, the crystal structure of the samples investigated remains the same.This discrepancy is most likely due to the difference in the sample preparation methods followed.Qian et al. [17] argued that the difference in annealing temperatures is accountable for the discrepancy in T c 's.Thus, we also anticipate that the lower T c in the HfRuP sample used in the present study is because of a different heat treatment procedure.Furthermore, the ratio of the superconducting gap to T c was estimated to be (2∆ 0 /k B T c )∼ 3.38(4), which is slightly smaller than the BCS value of 3.53.Our conclusion of the gap symmetry from the µSR study is in agreement with the temperature dependence of the magnetic penetration depth measurement using the tunnel-diode oscillator (TDO) method, which also found an s-wave superconductivity in HfRuP [30].We obtained an absolute value of the penetration depth λ(0) = 603(2) nm, which is slightly higher than that obtained from the TDO experiments [30].It is to be noted that the µSR method is more accurate than the TDO measurements in measuring the absolute value of the penetration depth.In Table 1, we have compared various superconducting parameters of different superconductors belonging to the TT X family.In London theory [31], the penetration depth, λ(0), is related to different superconducting parameters such as the effective mass, m * , and the superconducting carrier density, n s , via the relation λ 2 (0)= m * /µ 0 n s e 2 .Here, m * can be determined using the relation m * = (1 + λ e−ph )m e , where λ e−ph is the electron-phonon coupling constant and m e is the rest electron mass.
Using McMillan's relation, it is also possible to determine the electron-phonon coupling constant (λ e−ph ) [32]: where µ * is the repulsive screened Coulomb parameter, usually assigned as µ * = 0.13, and Θ D is the Debye temperature.As the value of Θ D for HfRuP is not available in the literature, we have used the Θ D = 345 K of ZrRuP [33,34], and using the scaling relation given in ref. [35] we have estimated the Θ D =294.6 K for HfRuP, and using Equation ( 4) we have estimated m * =1.799m e , where m e is the mass of the free electron for HfRuP.Using this value of m * , we estimated the value of n s = 1.41 × 10 26 carriers/m 3 .This value of n s is comparable to that observed in other TT X members such as ZrRuAs (n s = 2.1 × 10 26 m −3 ) [36], HfIrSi (n s = 6.6 × 10 26 m −3 ) [37], and ZrIrSi (n s = 6.9 × 10 26 m −3 ) [38].Notably, this value of n s is of the same order of magnitude as reported in some other TSCs, e.g., Nb 0.25 Bi 2 Se 3 (n s = 0.25 × 10 26 m −3 ) [39] and T d -MoTe 2 (n s = 1.67 × 10 26 m −3 ) [40].It is worth highlighting that the relatively high value of T c and low value of n s in HfRuP evinces possible unconventional superconductivity in this compound.The T c /λ −2 (0) ratio is often considered as a crucial parameter to qualitatively address any unconventional origin of the superconductivity.For HfRuP, the ratio T c [K]/λ −2 (0) µm −2 is estimated to be ∼3.3, which resides between the values observed for electron-doped (T c /λ −2 (0) ∼ 1) and holedoped (T c /λ −2 (0) ∼ 4) cuprates [41][42][43].This is indicative of a possible unconventional mechanism of superconductivity in HfRuP.
The presence of a single s-wave gap is further corroborated by field-dependent TF measurements at 1.6 K.During these measurements, for each field point we field cooled the sample from 15 K (above T c ) to 1.6 K.The field dependence of the relaxation rate σ sc (B) is presented in Figure 3b.The observed σ sc (B) can be described using the Brandt formula (for an s-wave superconductor) [26], Thus, we obtained an upper critical field of µ 0 H c2 (0) = 16(3) T, which is very similar to that obtained from the electrical resistivity measurements [17].

ZF-µSR Measurements
We also performed ZF-µSR experiments at temperatures above and below T c , to test whether the time-reversal symmetry is preserved or broken in HfRuP.As seen from Figure 4, the ZF-µSR spectra can be well described by a damped Gaussian Kubo-Toyabe depolarization function [44] A ZF (t) = A 0 G KT exp(−Λt), where A 0 is the initial asymmetry, G KT is the Gaussian Kubo-Toyabe (KT) function [44] representing the contribution from an isotropic Gaussian distribution of randomly oriented static (or quasistatic) local fields at the muon sites, and Λ is the exponential relaxation rate.We do not observe any noticeable difference in the spectra collected at 10 K and at 0.27 K, suggesting the absence of any spontaneous field in the superconducting state and, hence, the time-reversal symmetry is preserved in the superconducting state of HfRuP.

Summary
The temperature dependence of the magnetic penetration depth determined using TF-µSR experiments confirm a single s-wave gap structure in HfRuP.The gap to T c ratio, 2∆(0)/k B T c , is slightly smaller than that expected for a BCS superconductor.The results from field-dependent TF measurements also confirm single-gap s-wave superconductivity.The ZF-µSR results suggest that time-reversal symmetry is preserved in the superconducting state of HfRuP.Although the admixture of spin-singlet and spin-triplet pairing states is allowed in the noncentrosymmetry crystal structure of HfRuP, our TF-µSR study, and also TDO measurements [30], do not suggest the presence of a node or spin-triplet pairing.Similar conclusions on the s-wave gap symmetry and preservation of the TRS have been reported in the noncentrosymmetric topological superconductor ZrRuAs [36].The results of the present study will be important to understand the superconductivity in other topological superconductors.

Figure 1 .
Figure 1.(a) Room temperature Rietveld-fitted XRD pattern of HfRuP powder sample.(b,c) The hexagonal crystal structure of HfRuP.The red, purple, green, and cyan colored atoms are Hf, Ru, and the two types of P environment.(d) Temperature dependence of dc magnetic susceptibility, collected under ZFC and FCC conditions.

Figure 2 .
Figure 2. (a) Transverse-field (TF)-µSR time spectra obtained above (20 K-red) and below (0.27 K-blue) T c = 9.2 K for HfRuP (after field cooling the sample from above T c ) in 30 mT.(b) Fourier transform of the µSR time spectra from panel (a) at 20 K (red) and 0.27 K (blue).

Figure 3 .
Figure 3. (a) Temperature dependence of inverse squared London penetration depth (λ −2 ) for HfRuP, measured in an applied field µ 0 H = 30 mT.The solid line corresponds to the fitting of the experimental data using an isotropic single s-wave gap model.(b) Field dependence of the superconducting muon spin depolarization rate at 1.6 K.The solid red line represents fitting with an isotropic single s-wave gap model.

Table 1 .
Superconducting parameters of some superconductors belonging to the TT X family.