Spin-triplet superconductivity in Weyl nodal-line semimetals

Topological semimetals are three dimensional materials with symmetry-protected massless bulk excitations. As a special case, Weyl nodal-line semimetals are realized in materials having either no inversion or broken time-reversal symmetry and feature bulk nodal lines. The 111-family, including LaNiSi, LaPtSi and LaPtGe materials (all lacking inversion symmetry), belongs to this class. Here, by combining muon-spin rotation and relaxation with thermodynamic measurements, we find that these materials exhibit a fully-gapped superconducting ground state, while spontaneously breaking time-reversal symmetry at the superconducting transition. Since time-reversal symmetry is essential for protecting the normal-state topology, its breaking upon entering the superconducting state should remarkably result in a topological phase transition. By developing a minimal model for the normal-state band structure and assuming a purely spin-triplet pairing, we show that the superconducting properties across this family can be described accurately. Our results demonstrate that the 111 materials reported here provide an ideal test-bed for investigating the rich interplay between the exotic properties of Weyl nodal-line fermions and unconventional superconductivity.


INTRODUCTION
Topological materials are at the forefront of current condensed matter and material science research due to their great potential for applications. Among their defining characteristics is the symmetry-protected metallic surface state, arising from a nontrivial bulk topology. Recently, the experimental observation of many topological semimetals has shifted the research focus towards this subclass of topological materials 1,2 . Contrary to Diracor Weyl semimetals, which have point-type band crossings, in nodal-line semimetals band crossings occur in the form of lines or rings along special k-directions of the Brillouin zone. In this case, near the nodes, the low-energy excitations are nodal-line fermions with rather exotic properties 1,2 . Weyl nodal-line semimetals can be realized in systems lacking inversion symmetry or with broken time-reversal symmetry (TRS), provided the nodal lines that are protected by additional symmetries. Recently, the isostructural noncentrosymmetric 111-type materials LaNiSi, LaPtSi, and LaPtGe have been predicted to be Weyl nodal-line semimetals, protected by nonsymmorphic glide planes 3 . In addition, at low temperatures, all of them become superconductors [4][5][6] .
The breaking of additional symmetries in the superconducting state, besides the global gauge symmetry of the wave function, is a key characteristic of unconventional superconductors 7,8 . The combination of intriguing fundamental physics with far-reaching potential for applications has made unconventional superconductors one of the most investigated classes of materials. Broken time reversal symmetry in the superconducting state, one of the typical indications of unconventional superconductivity (SC), is manifested by the spontaneous appearance of magnetic fields below the superconducting transition temperature T c 9 . Recently, by using the muon-spin relaxation technique, several noncentrosymmetric superconductors (NCSCs) have been found to break TRS in their superconducting state. Otherwise they appear to exhibit the conventional properties of standard phonon-mediated superconductors [9][10][11][12][13][14] . In NCSCs, singlet-triplet admixed pairings can be induced by antisymmetric spin-orbit coupling (ASOC), however, ASOC itself cannot break TRS 15,16 . Noncentrosymmetric superconductors also provide a fertile ground for topological superconductivity, with potential applications to topological quantum computing [17][18][19] .
According to electronic band-structure calculations, the ASOC strength increases progressively from LaNiSi to LaPtSi to LaPtGe 3 . Hence, the 111-family of materials is a prime candidate for investigating the relationship between ASOC and unconventional SC with TRS-breaking, here made even more interesting by the interplay with the exotic nodal-line fermions. Recent muon-spin relaxation and rotation (μSR) studies on LaNiSi and LaPtSi reported an enhanced muon-spin relaxation at low temperatures, seemingly an indication of TRS breaking. However, their unusual temperature dependence (here resembling a Curie-Weiss behavior), the lack of any distinct features near T c 20 , and the absence of an additional muon-spin relaxation in LaPtGe (below its T c ) 21 , all seem to suggest that TRS is preserved in the superconducting state of these 111 materials. We recall that, in the past, inconsistent μSR results have been also reported in UPt 3 22,23 , whose TRS breaking could be independently proved by optical Kerr effect only a decade later 24 . Clearly, it is highly desirable to investigate the 111 materials, too, with other techniques such as the Kerr effect, in order to confirm their TRS breaking. Here, by combining extensive and thorough μSR measurements with detailed theoretical analysis, we show that, contrary to previous reports, all the above 111-type materials spontaneously break TRS at the superconducting transition and exhibit a fully-gapped pure spin-triplet pairing.

Bulk superconductivity
We synthesized three isostructural LaNiSi, LaPtSi, and LaPtGe samples and investigated systematically their physical properties via magnetic-susceptibility, specific-heat, electrical-resistivity, and μSR measurements. As shown in Fig. 1a, the 111-type materials crystallize in a noncentrosymmetric body-centered tetragonal structure. The corresponding I4 1 md space group (No. 109), confirmed by refinements of the powder x-ray diffraction (XRD) patterns (see e.g., in Fig. 1b), is nonsymmorphic and has a Bravais lattice with point group C 4v (4mm) (see details in Supplementary Note 1 and Table 1). Upon zero-field cooling, full diamagnetic screening (i.e., bulk SC) is found in the magnetic susceptibility measurements in an applied field of 1 mT (Fig. 1c). Consistent with previous studies 4-6 , we find T c = 1.28, 3.62, and 3.46 K for LaNiSi, LaPtSi, and LaPtGe, respectively. A prominent specific-heat jump at each superconducting transition (see below) confirms once more the bulk SC nature of these materials.

Lower and upper critical fields
For reliable transverse-field (TF) μSR measurements in a superconductor, the applied magnetic field should exceed the lower critical field H c1 and be much less than the upper critical field H c2 , so that the additional field-distribution broadening due to the fluxline lattice (FFL) can be quantified from the muon-spin relaxation rate, the latter is directly related to the magnetic penetration depth and thus, to the superfluid density. The H c1 values determined from field-dependent magnetization data are summarized in Fig. 2a-c, which provide lower critical fields μ 0 H c1 (0) = 3.9(5), 9.6(2), and 11.8 (1) mT for LaNiSi, LaPtSi, and LaPtGe, respectively. These H c1 (0) values are fully consistent with those determined from magnetic penetration depth (see below). We investigated also the upper critical fields H c2 of these 111 materials, here shown in Fig. 2d-f versus the reduced temperature T/T c (0) for LaNiSi, LaPtSi, and LaPtGe, respectively. Three different models, including Ginzburg-Landau (GL) 25 , Werthamer-Helfand-Hohenberg (WHH) 26 , and twoband model 27 were used to analyze the H c2 (T) data. In LaNiSi and LaPtGe, H c2 (T) is well described by the WHH model, yielding μ 0 H c2 (0) = 0.10(1) and 0.44(1) T, respectively. Conversely, in LaPtSi, both WHH and GL models fit H c2 (T) reasonably well only at low fields, i.e., μ 0 H c2 < 0.2 T. At higher magnetic fields, both models  Fig. 1. c Temperature dependence of the magnetic susceptibility. The zero field-cooled (ZFC) and field-cooled (FC) magnetic susceptibility were measured in a field of μ 0 H = 1 mT. The well distinct ZFC-and FC curves are consistent with type-II SC, as confirmed also by μSR measurements.
T. Shang et al. deviate significantly from the experimental data. Such a discrepancy most likely hints at multiple superconducting gaps in LaPtSi, as evidenced also by the positive curvature of H c2 (T), a typical feature of multigap superconductors. Indeed, here the two-band model shows a remarkable agreement with the experimental data and provides μ 0 H c2 (0) = 1.17(2) T. The presence of multiple superconducting gaps is also supported by the field-dependent electronic specific-heat coefficient (see Supplementary Fig. 11 and Note 7). However, note that, due to the similar gap sizes (or to small relative weights) of LaPtSi, the multigap features are not easily discernible in the superfluid density or in the zero-field electronic specific heat 28,29 .

ZF-μSR and evidence of TRS breaking
Zero-field (ZF)-μSR is a very sensitive method for detecting weak magnetic fields (down to~0.01 mT 30 ) due to the large muon gyromagnetic ratio (851.615 MHz T −1 ) and to the availability of nearly 100% spin-polarized muon beams. Therefore, the ZF-μSR technique has been successfully used to study different types of unconventional superconductors with broken TRS in their superconducting state 10,11,13,14,22,[31][32][33][34] . To search for the presence of TRS breaking in the superconducting state of 111 materials, ZF-μSR measurements were performed at various temperatures, covering both their normal-and superconducting states. Representative ZF-μSR spectra are shown in Fig. 3a-c for LaNiSi, LaPtSi, and LaPtGe, respectively. The ZF-μSR spectra exhibit small yet clear differences between 0.02 K and temperatures above T c (e.g., 1.9 K) for LaNiSi, which become more evident in the LaPtSi and LaPtGe case.
In general, in absence of external magnetic fields, the muonspin relaxation is mostly determined by the interaction of muon spins with the randomly oriented nuclear magnetic moments. Thus, the ZF-μSR asymmetry can be described by means of a phenomenological relaxation function, consisting of a combination of Gaussian-and Lorentzian Kubo-Toyabe relaxations [see Eq.
(2)] 35,36 . While σ ZF (T) is found to be nearly temperature independent (see Supplementary Fig. 7), as shown in Fig. 3d-f, all three compounds show a clear increase of the muon-spin relaxation in the Λ ZF channel below T c . Conversely, for T > T c , Λ ZF (T) is flat, thus excluding a possible origin related to magnetic impurities (the latter typically follow a Curie-Weiss behavior 28 ). Furthermore, longitudinal-field (LF) μSR measurements at base temperature (see Fig. 3a-c) indicate that a field of only 10 mT is sufficient to decouple the muon spins from the TRS breaking relaxation channel in all three compounds, indicating that the weak internal fields are static within the muon lifetime. Furthermore, the LF-μSR results rule out an extrinsic origin for the enhanced Λ ZF (T). Considered together, the ZF-μSR and LF-μSR results clearly evidence the increase in Λ ZF (T) below T c with the occurrence of spontaneous magnetic fields 10,11,13,14,22,[31][32][33][34] and, hence, the breaking of TRS in the superconducting state of LaNiSi, LaPtSi, and LaPtGe.

TF-μSR and nodeless superconductivity
To investigate the superconducting order parameters of LaNiSi, LaPtSi, and LaPtGe, the temperature dependence of their magnetic penetration depth was determined via TF-μSR measurements. The development of a flux-line lattice in the mixed state of a superconductor broadens the internal field distribution and leads to an enhanced muon-spin relaxation rate. Since the latter is determined by the magnetic penetration depth and, ultimately, by the superfluid density, the superconducting order parameter can be evaluated from the temperature-dependent TF-μSR measurements (see "Methods" section). Following a field-cooling protocol down to 0.02 K, the TF-μSR spectra were collected at various temperatures upon warming, covering both the superconducting and the normal states. As shown in Fig. 4a-c, below T c , the fast decay in the TF-μSR asymmetry caused by the FLL is clearly visible. By contrast, the weak decay in the normal state, is attributed to the nuclear magnetic moments, being similar to the ZF-μSR in Fig.  3a-c. The TF-μSR spectra were analyzed by means of Eq. (3). Above T c , the relaxation rate is small and temperature-independent, but below T c it starts to increase due to the formation of a FLL and the corresponding increase in superfluid density. At the same time, a diamagnetic field shift appears below T c (see Fig. 4d-f). The effective magnetic penetration depth and the superfluid density were calculated from the measured superconducting Gaussian relaxation rates (see "Methods" section). The normalized inversesquare of the effective magnetic penetration depth λ À2 eff ðTÞ (proportional to the superfluid density) vs. the reduced temperature T/T c for LaNiSi, LaPtSi, and LaPtGe is presented in Fig. 4g-i, respectively. Although these three NCSCs exhibit different T c values and ASOC strengths, below T c /3, their λ À2 eff values are practically independent of temperature. The low-T invariance of λ À2 eff ðTÞ and, consequently, of the superfluid density, clearly suggests the lack of low-energy excitations and, hence, a nodeless superconductivity in LaNiSi, LaPtSi, and LaPtGe, in good agreement with the low-T electronic specific-heat data (see below) and magnetic penetration depth measurements via the tunnel-diodeoscillator technique. The solid lines through the data in Fig. 4g-i are fits using a fully-gapped s-wave model with a single superconducting gap. These yield gap values Δ 0 = 1.95(5), 1.80 (5), and 2.10(5) k B T c , and λ 0 = 300(3), 228 (3), and 219(2) nm for LaNiSi, LaPtSi, and LaPtGe, respectively. When employing the dirty-limit model we find similar parameters (see Supplementary   Fig. 9 and Note 6). The gap values determined via TF-μSR are highly consistent with those derived from the specific-heat measurements (see Supplementary Table 5).
Minimal two-band model and electronic specific heat The three 111 materials, LaNiSi, LaPtSi, and LaPtGe share similar band structures and are inherently multiband systems, with several orbitals contributing to the density of states (DOS) at the Fermi level 3 . The lack of an inversion center implies that an antisymmetric spin-orbit coupling is naturally present in these materials. Here, the ASOC splits the bands near the Fermi level, with an increasingly larger strength from LaNiSi to LaPtSi to LaPtGe 3 . In their normal state, all of them are nonmagnetic and, thus, preserve TRS. These materials have been predicted to exhibit four Weyl nodal rings around the X point, at~0.5 eV below the Fermi level, a topological feature protected by nonsymmorphic glide mirror symmetry and TRS 3 . The size of the nodal rings increases with increasing ASOC strength and, due to their presence, the 111-type Weyl nodal-line semimetals are expected to show interesting magneto-transport properties. The antisymmetric spin-orbit coupling, however, does not change the topology of the Fermi surfaces qualitatively. In the 111 materials, in the absence of ASOC, three spin-degenerate bands cross the Fermi level giving rise to three Fermi surfaces with similar shapes. However, only two of them contribute significantly (~96%) to the DOS at the Fermi level (see Supplementary Fig. 12 and Note 9) 3 . The low-energy properties of 111 materials are thus dominated by these two Fermi surfaces. To capture qualitatively their topology, we construct a minimal two-band tight-binding model by suitably choosing the chemical Fig. 3 Probing TRS breaking superconductivity via ZF-μSR. a-c Zero-field μSR spectra collected above and below T c in LaNiSi (a), LaPtSi (b), and LaPtGe (c). In all cases, the lack of any oscillations implies a lack of long-range magnetic order. Solid lines through the data in a-c are fits to Eq. (2). The flat μSR datasets in a-c correspond to LF-μSR spectra, suggesting that even a small longitudinal field is sufficient to decouple muon spins from the local field. d-f Zero-field muon-spin relaxation rate Λ ZF versus the reduced temperature T/T c for LaNiSi (d), LaPtSi (e), and LaPtGe (f). Solid lines in d-f are guides to the eyes. A consistent increase of Λ ZF below T c reflects the onset of spontaneous magnetic fields, indicative of a breaking of TRS in the superconducting state, while the σ ZF is almost temperature independent (see details in Supplementary  Fig. 7 and Note 5). The error bars of Λ ZF are the SDs obtained from fits to Eq. (2) by the musrfit software package 52 .
To account for the ASOC effects in the minimal two-band model, we note that the form of ASOC for the corresponding C 4v point group is V ASOC ¼ α xy ðk y σ x À k x σ y Þ þ α z k x k y k z ðk 2 x À k 2 y Þσ z , where σ = (σ x ,  σ y , σ z ) is the vector of Pauli matrices in the spin space, while α xy and α z are the strengths of the two types of ASOC terms allowed by symmetry 15 . Note that, the second term is of fifth order in k and leads to spin splitting. On the other hand, the Rashba term, with strength α xy , is expected to be dominant because of the quasi-2D nature of the two Fermi surfaces (see Fig. 5a). Hence, only the Rashba ASOC term is phenomenologically relevant in the minimal model. In general, this term would have both interband and intraband contributions. However, to maintain the topology of the Fermi surfaces in presence of ASOC similar to that in its absence and to correctly describe the experimental observations in the 111 materials, we need to work in the limit where the interband contribution is large compared to the intraband one (see Supplementary Table 2  Due to the inherent multiband nature and to the presence of nonsymmorphic symmetries, the usual classification of the possible superconducting order parameters (based on pointgroup symmetries) in the effective single-band picture is insufficient for the 111 materials (see Supplementary Note 9 for details). Indeed, nonsymmorphic symmetries can lead to additional symmetry-imposed nodes along the high symmetry directions on the zone faces 37,38 . Even a loop supercurrent state 39 , which has a uniform onsite singlet pairing and proposed to be realized in some of the fully-gapped TRS breaking superconductors, is not allowed in the case of 111 materials, because there are only two symmetrically distinct sites within a unit cell. However, we note that the two Fermi surfaces under consideration have large sections in the Brillouin zone which are almost parallel and close to each other (see Fig. 5b). Hence, to consistently explain the phenomenon of TRS breaking at T c despite the presence of a full SC gap, we expect that an internally antisymmetric nonunitary triplet (INT) state 40,41 , which features a uniform pairing between same spins in the two different bands, to become the dominant instability. In this state, the pairing potential matrix isΔ ¼Δ S Δ B , whereΔ S andΔ B are the pairing potential matrices in the spinand band-space, respectively.Δ B ¼ iτ y gives the required fermionic antisymmetry.Δ S ¼ ðd:σÞiσ y , where d = Δ 0 η, with |η| 2 = 1, is the d-vector characterizing the triplet pairing state, which is nonunitary because q = i(η × η * ) ≠ 0. Δ 0 is an overall pairing amplitude. Similarly, nonunitary triplet superconducting ground state was also proposed to be realized in the recently discovered 3D Dirac semimetals LaNiGa 2 42 . We compute the Bogoliubov quasiparticle energies E n (k), n = 1, …4 for the effective model in the INT ground state using the Bogoliubov-de-Gennes (BdG) formalism (see details also in Supplementary Note 9, Figs. 13-15, Tables 3 and 4). The thermodynamic properties are computed by assuming that the temperature dependence comes only from the pairing amplitude in the form ΔðTÞ ¼ Δ 0 tanhf1:82½1:018ðT c =T À 1Þ 0:51 g and ignoring any weak temperature dependence of the q-vector. To reproduce the experimental specific-heat results for the three materials, three fitting parameters, namely, Δ 0 /(k B T c ), the direction of d-vector, and α had to be tuned to get the best fits in the weakcoupling limit (see Fig. 6). Note that, for all the three materials we can reproduce the specific-heat data rather well (especially at low temperatures) and the fitting process naturally preserves the trend of increasing ASOC strength in the 111 family. More importantly, the derived superconducting energy gaps are highly consistent with the values determined from TF-μSR measurements. The nonzero real vector q, found from the fits, points in different directions for the three materials and encodes the effective TRSbreaking field arising from spin-polarization caused by Cooperpair migration due to the nonunitary nature of pairing 41 .

DISCUSSION
According to ZF-μSR results in the 111 materials, the spontaneous magnetic fields or the magnetization in the superconducting state of LaPtSi or LaPtGe are much larger than in LaNiSi, here reflected in significantly larger variations of Λ ZF between zero-temperature and T c in the former two cases as compared to LaNiSi. Therefore, the TRS breaking effect is more prominent in the superconducting state of LaPtSi and LaPtGe than in LaNiSi (see Fig. 3). Previous ZF-μSR studies indicate that, although LaNiSi and LaPtSi exhibit an enhanced muon-spin relaxation rate at low temperatures, their Λ ZF (T) resembles a Curie-Weiss behavior [i.e., λ ZF (T) ∝ T −1 ]. This, and the lack of a distinct anomaly in Λ ZF (T) across T c 20 , are inconsistent with a TRS breaking effect. In general, an enhanced muon-spin relaxation with Curie-Weiss features might be related to either intrinsic or to extrinsic spin fluctuations. As for the intrinsic case, a typical example is that of the ThFeAsN iron-based superconductor. It exhibits strong magnetic fluctuations at low temperatures (confirmed also by nuclear magnetic resonance measurements), which are reflected in a steadily increasing Λ ZF (T) as the temperature is lowered 43 . As for the extrinsic case, a typical example is that of the ReBe 22 multigap superconductor. Here, Λ ZF (T) increases remarkably with decreasing temperature due to tiny amounts of magnetic impurities, whose contribution is enhanced near zero temperature 28 . Conversely, in case of a truly Fig. 6 Electronic specific heat. a-c Normalized electronic specific heat C e /γT (with γ the normal-state electronic specific-heat coefficient) as a function of the reduced temperature T/T c for LaNiSi (a), LaPtSi (b), and LaPtGe (c). C e /T was determined by subtracting the phonon contribution from the raw specific-heat data (see Supplementary Fig. 10 and Note 7). Solid lines through the data represent theoretical calculations corresponding to the INT state with the fit parameters listed in the figures. For LaPtSi, the reduced specific-heat jump at T c is mostly caused by the broadening of the superconducting transition.
broken TRS-for instance, in the 111 materials we report here-Λ ZF is almost independent of temperature for T > T c , strongly suggesting that the enhanced Λ ZF is induced by the spontaneous fields occurring in the superconducting state.
In LaPtGe, previous ZF-μSR data show similar features to our results (see Fig. 3f), i.e., a small yet clear difference in the ZF-μSR spectra between 0.3 and 4.5 K 21 , the latter dataset referring to the normal state. However, the authors claimed that, the temperaturedependent σ ZF and Λ ZF exhibit no visible differences and, thus, a preserved TRS was concluded 21 . By contrast, our systematic ZF-μSR measurements suggest the presence of spontaneous magnetic fields, hence, the broken TRS in the superconducting state of LaPtGe. Such discrepancies in ZF-μSR results might be related to the different sample quality, purity, or disorder. For example, the previous study reports a residual resistivity ρ 02 00 μΩ cm 21 , three times larger than that of current LaPtGe sample, ρ 0~6 6 μΩ cm (see Supplementary Fig. 2 and Note 2). Moreover, the residual resistivity ratio of the current LaPtGe sample is twice larger than that of the previous sample. Nevertheless, to independently confirm the TRS breaking in the superconducting state of 111 materials, the use of other techniques, as e.g., Josephson tunneling, SQUID, or optical Kerr effect, is highly desirable. In particular, the optical Kerr effect, another very sensitive probe of spontaneous fields in unconventional superconductors, is renown for confirming TRS breaking in Sr 2 RuO 4 and UPt 3 24,44 . In addition, to exclude disorder effects, search for possible non-s-wave behavior, and confirm the TRS breaking in the 111 materials, in the future, measurements on high-quality single crystals will clearly be helpful.
According to the Uemura plot 45 , clearly, the 111 materials studied here lie in the band where there are other families of superconductors that are found to break TRS in the superconducting state (see details in Supplementary Fig. 16 and Note 10). Apart from the La-based 111 materials studied here, also the isostructural ThTSi compounds (with T = Co, Ni, Ir, and Pt) are superconductors (with critical temperatures between 2 and 6.5 K) 46,47 . Similar to the La-based cases, the Th-based materials, too, exhibit a large ASOC upon replacing the 3d Ni and Co with 5d Pt and Ir 48 . Recently, superconductivity with T c = 5.07 K was reported in ThIrP, which also adopts a LaPtSi-type structure 49 . Therefore, it would be interesting to search for possible TRS breaking and, hence, unconventional superconductivity in these Th-based 111 materials. In addition, La-based 111 materials, in particular LaNi 1−x Pt x Si, represent ideal candidate systems for investigating the effect of ASOC on spontaneous magnetization and unconventional superconductivity.
Generally, in noncentrosymmetric superconductors, the ASOC can induce a mixing of singlet and triplet states. However, in the 111 materials under consideration, it does more than that, becoming crucial in stabilizing even a purely triplet state. Moreover, the necessity of a dominant interband contribution to the ASOC in achieving a fully gapped spectrum in the INT state, further justifies the interband pairing in the superconducting state. We also note that the triplet d-vectors, obtained from analyses of the specific-heat data, correspond to a partially spinpolarized (|q| < 1 ≠ 0) superconducting state. In this case, the spontaneous magnetization results from a migration of Cooper pairs from the majority to the minority spin species 41 .
The normal state of 111 materials has a non-trivial topology due to the Weyl nodal lines protected by the nonsymmorphic glide symmetry and TRS. Apart from the usual photoemission studies 1,2 , the corresponding drumhead surface states can also be investigated by inspecting the correlation effects on the surfaces 50 . Since TRS is spontaneously broken at T c , it is of interest to investigate the fate of the bulk nodal lines. Our results demonstrate that 111 materials represent a rare case of Weyl nodal-line semimetals, which break time-reversal symmetry in the superconducting state. As such, they epitomize the ideal system for investigating the rich interplay between the exotic properties of topological nodal-line fermions and unconventional superconductivity.

Sample preparation
Polycrystalline LaNiSi, LaPtSi, and LaPtGe samples were prepared by arc melting La (99.9%, Alfa Aesar), Ni (99.98%, Alfa Aesar), Pt (99.9%, ChemPUR), Si (99.9999%, Alfa Aesar) and Ge (99.999%, Alfa Aesar) in high-purity argon atmosphere. To improve homogeneity, the ingots were flipped and re-melted more than five times. The as-cast ingots were then annealed at 800°C for two weeks. The crystal structure and purity of the samples were checked using powder x-ray diffraction at room temperature using a Bruker D8 diffractometer with Cu K α radiation. All three compounds crystallize in a tetragonal noncentrosymmetric structure with a space group I4 1 md (No. 109). The estimated lattice parameters are listed in the Supplementary Table 1.

Sample characterization
The magnetization, heat-capacity, and electrical-resistivity measurements were performed on a Quantum Design magnetic property measurement system (MPMS) and a physical property measurement system (PPMS). The lower critical field H c1 was determined by field-dependent magnetization measurements at various temperatures up to T c , while the upper critical field H c2 was determined by measuring the temperature-dependent heat capacity under various magnetic fields, and by field-dependent magnetization at various temperatures.

μSR experiments
The μSR experiments were conducted at the general-purpose surfacemuon (GPS) and at the low-temperature facility (LTF) instruments of the Swiss muon source (SμS) at Paul Scherrer Institut (PSI) in Villigen, Switzerland. Once implanted in a material, at a typical depth of 0.3 mm, the spin-polarized positive muons (μ + ) act as microscopic probes of the local magnetic environment via the decay positrons, emitted preferentially along the muon-spin direction. The spatial anisotropy of the emitted positrons (i.e., the asymmetry signal) reveals the distribution of the local magnetic fields at the muon stopping sites 36,51 . For TF-μSR measurements, the applied magnetic field is perpendicular to the muonspin direction, while for LF-μSR measurements, the magnetic field is parallel to the muon-spin direction. In both the TF-and LF-μSR cases, the samples were cooled in an applied magnetic field down to the base temperature (1.5 K for GPS and 0.02 K for LTF). Field cooling reduces flux pinning and ensures an almost ideal flux-line lattice. The μSR spectra were then collected upon heating. For the ZF-μSR measurements, to exclude the possibility of stray magnetic fields, the magnets were quenched before the measurements, and an active field-nulling facility was used to compensate for stray fields down to 1 μT.

Analysis of the μSR spectra
All the μSR data were analyzed by means of the musrfit software package 52 . In absence of applied external fields, in the nonmagnetic LaNiSi, LaPtSi, and LaPtGe, the relaxation is mainly determined by the randomly oriented nuclear magnetic moments. Therefore, their ZF-μSR spectra can be modeled by means of a phenomenological relaxation function, consisting of a combination of Gaussian-and Lorentzian Kubo-Toyabe relaxations 35,36 : Here A s and A bg represent the initial muon-spin asymmetries for muons implanted in the sample and the sample holder, respectively. The σ ZF and Λ ZF represent the zero-field Gaussian-and Lorentzian relaxation rates, respectively. Since σ ZF shows an almost temperature-independent behavior, the Λ ZF values in Fig. 3a-c could be derived by fixing σ ZF to its average value, i.e., σ av ZF ¼ 0:094, 0.103, and 0.104 μs −1 for LaNiSi, LaPtSi, and LaPtGe, respectively (see details in Supplementary Fig. 7).
In the TF-μSR case, the time evolution of asymmetry was modeled by: T. Shang et al.
Here A s and A bg are the same as in ZF-μSR. B s and B bg are the local fields sensed by implanted muons in the sample and the sample holder (i.e., silver plate), γ μ is the muon gyromagnetic ratio, ϕ is the shared initial phase, and σ is a Gaussian relaxation rate reflecting the field distribution inside the sample. In the superconducting state, σ includes contributions from both the flux-line lattice (σ sc ) and a smaller, temperature-independent relaxation, due to the nuclear moments (σ n , similar to σ ZF ). The former can be extracted by subtracting the σ n in quadrature, i.e., σ sc = ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi σ 2 À σ 2 n p . Since σ sc is directly related to the effective magnetic penetration depth and, thus, to the superfluid density (σ sc / 1=λ 2 eff $ ρ sc ), the superconducting gap and its symmetry can be investigated by measuring the temperature-dependent σ sc .
The effective penetration depth λ eff could be calculated from σ sc by considering the overlap of vortex cores 53 :

Superconducting gap symmetry
Since the 111 materials exhibit an almost temperature-independent superfluid density below 1/3T c , to extract the superconducting gap, the superfluid density ρ sc (T) of LaNiSi, LaPtSi, and LaPtGe was analyzed by using a fully-gapped s-wave model, generally described by: with f ¼ ð1 þ e E=kBT Þ À1 the Fermi function 54,55 and λ 0 the effective magnetic penetration depth at zero temperature. The temperature evolution of the superconducting energy gap follows ΔðTÞ ¼ Δ 0 tanhf1:82½1:018ðT c =T À 1Þ 0:51 g, where Δ 0 is the gap value at zero temperature 56 .

Theoretical analysis
We use the Bogoliubov-de-Gennes formalism to compute the quasiparticle energy bands for the minimal two-band model in the INT state 8 . Here, we assume the temperature dependence comes only from the pairing amplitude. The Bogoliubov quasiparticle energy bands E n (k); n = 1…4 are used to compute the specific heat by means of the formula: C ¼ X n;k k B β 2 2 E n ðkÞ þ β ∂E n ðkÞ ∂β E n ðkÞ sech 2 βE n ðkÞ 2 where β ¼ 1 kBT and k B is the Boltzmann constant. Further details about the theoretical analysis can be found in the Supplementary Note 9.

DATA AVAILABILITY
All the data needed to evaluate the reported conclusions are presented in the paper and/or in the Supplementary Material. Additional data related to this paper may be requested from the authors. The μSR data were generated at the SμS (Paul Scherrer Institut, Switzerland). Derived data supporting the results of this study are available from the corresponding authors or beamline scientists. The musrfit software package is available online free of charge at http://lmu.web.psi.ch/musrfit/technical/ index.html.