Evidence of s-wave superconductivity in the noncentrosymmetric La7Ir3

Superconductivity in noncentrosymmetric compounds has attracted sustained interest in the last decades. Here we present a detailed study on the transport, thermodynamic properties and the band structure of the noncentrosymmetric superconductor La 7 Ir 3 (Tc ~ 2.3 K) that was recently proposed to break the time-reversal symmetry. It is found that La7Ir3 displays a moderately large electronic heat capacity (Sommerfeld coefficient γn ~ 53.1 mJ/mol K2) and a significantly enhanced Kadowaki-Woods ratio (KWR ~32 μΩ cm mol2 K2 J−2) that is greater than the typical value (~10 μΩ cm mol2 K2 J−2) for strongly correlated electron systems. The upper critical field Hc2 was seen to be nicely described by the single-band Werthamer-Helfand-Hohenberg model down to very low temperatures. The hydrostatic pressure effects on the superconductivity were also investigated. The heat capacity below Tc reveals a dominant s-wave gap with the magnitude close to the BCS value. The first-principles calculations yield the electron-phonon coupling constant λ = 0.81 and the logarithmically averaged frequency ωln = 78.5 K, resulting in a theoretical Tc = 2.5 K, close to the experimental value. Our calculations suggest that the enhanced electronic heat capacity is more likely due to electron-phonon coupling, rather than the electron-electron correlation effects. Collectively, these results place severe constraints on any theory of exotic superconductivity in this system.

electron-phonon coupling. However, other properties of this putative TRS-breaking state have never been reported thus far. In this context, it is necessary to study the transport and thermodynamic properties of this intriguing state and search for the evidence of TRS breaking state. For example, the triplet pairing with broken TRS often leads to an upper critical field far above the Pauli paramagnetic limit and has nodal or anisotropic gap functions.
In this study, we explored the noncentrosymmetric La 7 Ir 3 superconductor via ultra-low temperature transport and thermodynamic measurements, complementary with the first principles calculations. Remarkably, the upper critical field in this system is lower than the Pauli limit and can overall be described by the one-band Werthamer-Helfand-Hohenberg (WHH) model. A moderately large electronic specific heat γ n ~ 53.1 mJ/mol K 2 was seen in the normal state and quasiparticle excitations in the superconducting state can be well characterized by a single s-wave gap with the magnitude of 2Δ g /k B T c = 3.56, very close to the BCS value of 3.5. Besides, the Hall effect, high pressure effect and the band structure were also investigated, which consistently support a phonon-mediated BCS pairing scenario for La 7 Ir 3 .

Results and Discussion
The schematic view of the crystal structure of La 7 Ir 3 is shown in Fig. 1. La 7 Ir 3 crystallizes in a hexagonal Th 7 Fe 3 structure with the space group P6 3 mc (No. 186). Its structure consists of alternate stacking of iridium monolayer and lanthanum bilayer along the c-axis. We optimized the coordinates of atoms with the experimental lattice parameters a = b = 10.2376Å, c = 6.4692Å 19 . Fig. 2(a) shows the temperature dependent resistivity of La 7 Ir 3 . A sharp superconducting transition is observed around 2.3 K, which is consistent with the diamagnetization measurement shown in the inset of Fig. 2(a). The normal state resistivity of the sample displays a typical metallic behavior. For T below ~25 K, the resistivity is well fitted by the Fermi-liquid (FL) expression, ρ(T) = ρ 0 + AT 2 . The resultant ρ 0 and A are 2.5 μΩ m and 9 × 10 −2 μΩ cm K −2 , respectively. It is worth noting, this A value is almost one order of magnitude larger than iron-based superconductor BaFe 2 (As 1−x P x ) 2 (~9 × 10 −3 μΩ cm K −2 for .  x 0 49) 21 , three orders of magnitude larger than recently reported Dirac nodal line semimetal PtSn 4 (~2 × 10 −4 μΩ cm K −2 ) 22 and PdSn 4 (~7 × 10 −4 μΩ cm K −2 ) 23 , yet one order of magnitude smaller than heavy fermion superconductor CePt 3 Si (~2.35 μΩ cm K −2 ) 2 .
The specific heat plotted as C/T vs T 2 is depicted in Fig. 2(b). A jump in C/T around 2.3 K also signifies a sharp transition to the superconducting state. It is noted that C/T has a small residual value as T → 0, indicating the non-superconducting fraction of our sample. The non-superconducting counterpart accounts for ~10% of the sample in volume. Through fitting the C(T) data above T c to the formula C = γ n T + β n T 3 + α n T 5 , the Sommerfeld coefficient γ n representing the electron contribution is extracted as 53.1 mJ/mol K 2 . This γ n value is much larger than that in BaFe 2 (As 1−x P x ) 2 (~7 mJ/mol K 2 ) and in PtSn 4 (~4 mJ/mol K 2 ), but smaller than that in CePt 3 Si (~390 mJ/mol K 2 ). This large γ n is consistent with the enhanced A coefficient, suggesting considerable density of states at the Fermi level in the normal state. The corresponding Kadowaki-Woods ratio (KWR) γ A/ n 2 is estimated to be 32 μΩ cm mol 2 K 2 J −2 . Note that this value is even larger than those found in many strongly correlated metals 24,25 , typically ~1 0 μΩ cm mol 2 K 2 J −2 . The transport and thermodynamic parameters for the above mentioned compounds are summarized in Table 1 for comparison.
The Hall effects are also studied in Fig. 3. The Hall resistivity ρ xy is found to grow linearly with field up to 9 T (Fig. 3 inset). The Hall coefficient R H is calculated as the slope of ρ xy (B) curves. As seen, the value of R H is negative in all temperature range studied, indicating the dominant electron carriers. R H increases with the decreasing temperature before saturating at T < 20 K. The fitting obtained residual resistivity ρ 0 is 2.5 μΩ m, and transport coefficient A is 9 × 10 −2 μΩ cm K −2 . Inset shows the temperature dependence of magnetization for both zero-field cooling (ZFC) and field cooling (FC) processes. (b) Temperature dependence of specific heat plotted as C/T versus T 2 . The red solid line represents the fit to the normal state C based on the equation C/T = γ n + β n T 2 + α n T 4 . The obtained parameters are γ n = 53.1 mJ/mol K 2 , β n = 3.96 mJ/mol K 4 and α n = 0.072 mJ/mol K 6 , respectively.

Comp.
A μΩ cm K −2 γ n mJ/mol K 2 KWR μΩ cm mol 2 K 2 J −2  Figure 4 shows the pressure (P) effects on La 7 Ir 3 . With increasing P, the normal state resistivity decreases slightly, and the superconducting transition temperature is progressively suppressed. As shown in Fig. 4(b), T c is linearly suppressed with P. The pressure suppression rate for T c is ~−0.15 K/GPa. Temperature sweeps and field sweeps for the superconducting transitions are shown in Fig. 5(a-c). S1 and S2 labeled the two major samples studied in the transport measurements. The temperature is measured down to 160 mK in the dilution refrigerator. As shown in Fig. 5(d), the upper critical field (μ 0 H c2 ) is extracted from both B-and Tsweeps, using both 50% and 90% of ρ n criteria (ρ n is the normal state resistivity). Note that the abscissa of Fig. 5(d) is the reduced temperature T/T c . It is found that μ 0 H c2 (T) curves from different criteria and different samples nearly collapse onto a single curve. These curves can be nicely fitted by the one-band Werthamer-Helfand-Hohenberg (WHH) model in the dirty limit. The reason of an upward curvature seen in 90% criterion of S1 is not clear, which  may be an artefact due to the impurity phase or due to an inhomogeneous current distribution arising from contact-sample interface issues. The fit gives zero-temperature μ 0 H c2 (0) close to 1.6 T, which is apparently smaller than the Pauli paramagnetic limit = .
 H T 1 84 p c 4 T by a factor of 2-3. We now turn to the detailed analysis of specific heat below T c , which provides unambiguous evidence for the superconducting gap symmetry. In Fig. 6(a), ΔC = C(0T) − C(2T) is shown. At μ 0 H = 2 T, the jump of C due to superconducting transition is totally suppressed (data overlap with the red fitting line in Fig. 2(b)). As the heat capacity from phonon is independent of magnetic field, ΔC excludes the phonon contribution as well as that from the nonsuperconducting part. We fit ΔC with different gap functions as was done in refs [26][27][28] . In the BCS theory, the zero field electronic specific heat in the superconducting state is derived from the entropy S es which is written as . For simplicity, we assume p-wave gap function to be Δ = Δ 0 cosφ. The electronic specific heat is calculated by C es = T(∂S/∂T). As seen in Fig. 6(a), at first sight, all three models fit the experimental data well. However, a blow-up view can clearly distinguish the differences among different models (see inset in Fig. 6(a)). At low temperatures, both p-wave and d-wave models deviate significantly from the experimental data, while the s-wave model reproduces the data very well. The entropy-conserving construction at T c gives ΔC/γ n T c = 1.22, close to the weak coupling BCS value of 1.43. The s-wave fit gives γ n ~ 49 mJ/mol K 2 , slightly smaller than the value (γ n ~ 53.1 mJ/mol K 2 ) obtained in Fig. 2(b), which means small amounts of the sample still nonsuperconducting. A rough estimate based on the ratio between these two values lead to a superconducting volume fraction around 92%. From the s-wave fit, the resultant α = 1 indicates the weak coupling BCS gap. We delve further into the temperature dependence of ΔC/T in the low-T range (T < 1/3T c ). In this low temperature limit, one would expect C es ~ T 2 , i.e., ΔC/T ~ aT−γ n for a clean d-wave superconductor, and for the nodeless s-wave model. Here, Δ g is the energy gap at zero temperature. As shown in Fig. 6(b), the experimental data of ΔC/T versus T clearly deviates from a linear relation, i.e., d-wave model is an inadequate description for La 7 Ir 3 . The s-wave model can again fit the experimental data better. The resultant Δ g gives Δ g /k B T c = 1.78, which is also close to the weak coupling BCS value (~1.76).
We further calculated the electronic structure and the phonon dynamics of the La 7 Ir 3 . As both La and Ir are relatively heavy elements, we have taken into account the spin-orbital coupling (SOC) in the calculations. For band calculations without SOC, there are four bands crossing the Fermi level (E F ) along the high symmetry path in the first Brillouin zone (BZ) as shown in Fig. 7(a). The corresponding total and partial densities of states are shown in the right panel. The bands crossing E F are constructed by the hybridization of La 5d and Ir 5d orbitals, and La has the dominant contribute around the Fermi level. The calculated band structure with SOC is shown in Fig. 7(b), where the number of bands doubles due to the lifting of the degeneracy. The calculated total density of states at E F is ~21 eV −1 per cell (2 formula units) with SOC. The corresponding bare specific heat coefficient γ 0 is 25.7 mJ/mol K 2 . The resultant Fermi surfaces (FSs) are plotted in 6(c) and (d), without and with the inclusion of SOC respectively, showing complex three dimensional characteristics. Figure 8 summarizes the lattice dynamics properties of La 7 Ir 3 . In the left panel of Fig. 8, we show the calculated phonon dispersions. The acoustic modes are not strongly anisotropic. There are 60 phonon bands in total extending up to ~140 cm −1 and the point group at Γ point is C 6v . Γ modes can be decomposed as 10  10  7  3  3  7   1  2  1  2  1 2 , with E 1 and E 2 modes doubly degenerate. The frequency of each mode at Γ is listed in Table 2.
In the middle panel of Fig. 8, we present Eliashberg spectral function α 2 F(ω), and electron-phonon coupling λ(ω). In the right panel of the same figure, we show the atom-projected phonon DOS. Similarly to the electronic bands, the phonon branches have fewer dispersions in the z direction. Analyzing the the phonon eigenvectors reveals that there is no clear separation between in and out-of-plane vibrations, as often happens in layered compounds. The phonon DOS distributes continuously in the frequency range up to 120 cm −1 . A small phonon gap of 10 cm −1 occurs between 120 cm −1 and 130 cm −1 . The projected phonon DOS shows that vibrations of La and Ir occupy the same frequency range, and the eigenvectors have a strongly mixed character. The vibration of Ir dominates the frequency range below 50 cm −1 , due to its relatively lager atom mass. In the density-functional perturbation theoretical (DFPT) 29 calculations, the Eliashberg spectral function depends directly on the electron-phonon matrix element: ) Here, N k is the number of k points used in the summation, N(E F ) is the density of states at the Fermi level, and ω qν are the phonon frequencies. The electron-phonon matrix element | | ν + g n m k k q , ( ) 2 is defined by the variation in the self-consistent crystal potential. From the λ(ω) curves in Fig. 8 c ln with the Coulomb parameter μ * = 0.15, we finally obtain T c = 2.5 K, which is close to the experimental value T c = 2.3 K. It then follows that La 7 Ir 3 is very likely a conventional electron-phonon superconductor, which seems consistent with our experimental results. From the above calculated electronic specific heat coefficient γ 0 25.7 mJ/mol K 2 and the

Conclusion
The transport and thermodynamic properties of La 7 Ir 3 presented in our study show clear evidence for the dominant spin singlet pairing, in agreement with the μ SR measurements. TRS-breaking spin triplet state usually has gap nodes (or anisotropy) in the order parameter. In broken TRS compounds LaNiC 2 and LaNiGa 2 32,33 , recent study proposed a novel triplet pairing state, namely, pairing between electrons with the same spins but on different orbitals. This unusual pairing inevitably leads to two even-parity (s-wave) gaps 34 . In our La 7 Ir 3 compound, however, no evidence for the two s-wave gaps is observed from H c2 and the heat capacity.
In summary, we presented the detailed physical properties of a recently reported nonnoncentrosymmetric superconductor La 7 Ir 3 in which time reversal symmetry was proposed to be broken via the detection of quasi-static internal magnetic field below T c . The La 7 Ir 3 superconductor is found to show a notably large electronic heat capacity (Sommerfeld coefficient γ n ~ 53 mJ/mol K 2 ) and a strongly enhanced Kadowaki-Woods ratio (KWR ~ 32 μΩ cm mol 2 K 2 J −2 ). The analysis of temperature dependent upper critical field and the electronic specific heat suggests a more likely conventional s-wave superconductivity in La 7 Ir 3 . The first-principles calculations on the electronic structure and the electron-phonon properties confirmed the experimental results.

Methods
Polycrystalline La 7 Ir 3 samples were prepared by the arc melting method as previously reported 19 . The sample phase has been verified through powder x-ray diffraction measurements (Fig. 1). Electrical transport measurements were performed by a standard four-probe method on the PPMS-9 T system (Quantum Design). To obtain the upper critical fields, measurements on both 3 He cooling system and a dilution refrigerator have been performed to measure the resistive superconducting transitions under field well below T c . High-pressure electrical measurements were carried out on PPMS using a piston cylinder type pressure cell (type: HPC-33) with Daphne 7373 as the pressure transmission medium. The specific heat data were obtained using a relaxation method via  the PPMS-9 T system equipped with a 3 He cooling option. Thermometers and the addenda were well calibrated prior to the measurements of the sample. The first-principles calculations of the band structure and the lattice dynamic properties were carried out based on experimental crystal structure. The electronic structure calculations with high accuracy were performed using the full-potential linearized augmented plane wave (FP-LAPW) method implemented in the WIEN2K code 35 . The generalized gradient approximation (GGA) 36 was applied to the exchange-correlation potential calculation. The muffin tin radii were chosen to be 2.5 a.u. for both La and Ir. The plane-wave cutoff was defined by RK max = 7.0, where R is the minimum LAPW sphere radius and K max is the plane-wave vector cutoff. Lattice dynamic properties including phonon dispersion, phonon density of states and electron-phonon coupling were performed using the Quantum-ESPRESSO 37 code with ultrasoft pseudopotential method and the plane wave basis. The cutoffs were chosen as 50 Ry for the wave functions and 500 Ry for the charge density. The generalized-gradient approximation of Perdew-Burke-Ernzerhof (PBE) 36 was used for the exchange-correlation energy function. The electronic integration was performed over a 8 × 8 × 8 k-point mesh. Dynamical matrices and the electron-phonon interaction coefficients were calculated on a 4 × 4 × 4 q-point grid. A dense 24 × 24 × 24 grid was used for evaluating the accurate electron-phonon interaction matrices.