## Point Node Gap Structure of Spin-Triplet Superconductor UTe 2

Tristin Metz, 1 Seokjin Bae, 1 Sheng Ran, 1, 2 I-Lin Liu, 1, 2, 3 Yun Suk Eo, 1 Wesley T. Fuhrman, 1 Daniel F. Agterberg, 4 Steven Anlage, 1 Nicholas P. Butch, 1, 2 and Johnpierre Paglione 1, 5, ∗

1 Maryland Quantum Materials Center, Department of Physics,

University of Maryland, College Park, Maryland 20742, USA

2 3 Department of Materials Science and Engineering,

NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA

University of Maryland, College Park, MD 20742 USA

4 Department of Physics, University of Wisconsin, Milwaukee, Wisconsin 53201, USA

5 Canadian Institute for Advanced Research, Toronto, Ontario M5G 1Z8, Canada

(Dated: November 11, 2019)

Low-temperature electrical and thermal transport, magnetic penetration depth, and heat capacity measurements were performed on single crystals of the actinide superconductor UTe 2 to determine the structure of the superconducting energy gap. Heat transport measurements performed with currents directed along both crystallographic a - and b -axes reveal a vanishingly small residual fermionic component of the thermal conductivity. The magnetic field dependence of the residual term follows a rapid, quasi-linear increase consistent with the presence of nodal quasiparticles, rising toward the a -axis upper critical field where the Wiedemann-Franz law is recovered. Together with a quadratic temperature dependence of the magnetic penetration depth up to T/T c = 0 3, these measurements . provide evidence for an unconventional spin-triplet superconducting order parameter with point nodes. Millikelvin specific heat measurements performed on the same crystals used for thermal transport reveal an upturn below 300 mK that is well described by a divergent quantum-critical contribution to the density of states (DOS). Modeling this contribution with a T -1 / 3 power law allows restoration of the full entropy balance in the superconducting state and a resultant cubic power law for the electronic DOS below T c , consistent with the point-node gap structure determined by thermal conductivity and penetration depth measurements.

The uranium based superconductors URhGe, UCoGe, and UGe 2 have generated much interest due to the unusual coexistence of a ferromagnetic ground state with superconductivity present in all three systems [1]. Recently, the discovery of superconductivity in UTe 2 [2] has stimulated a flurry of experimental studies on this new actinide superconductor [2-7]. With a ground state close to ferromagnetic order as evidenced by quantum critical scaling of the magnetization [2] and magnetic fluctuations down to milliKelvin temperatures [7], the UTe 2 system adds the long-sought paramagnetic end member to the family of uranium-based ferromagnetic superconductors.

unitary pairing state is not formally possible given the orthorhombic crystal symmetry. More important, the entropy balance between normal and superconducting states that is required for a second-order transition ( i.e. , ∫ T c 0 C el s T dT = ∫ T c 0 γ n dT , where C el s is the electronic heat capacity in the superconducting state and γ n is the normal state electronic heat capacity) is not obeyed [2, 3].

glyph[negationslash]

The superconducting state of UTe 2 is quite exotic. Strong evidence for spin-triplet pairing is given by a lack of change in the NMR knight shift upon cooling into the superconducting state, and an upper critical field H c 2 that exceeds the paramagnetic limit for all field directions ( i.e. , by more than an order of magnitude for H ‖ b ) [2]. The NMR 1/T 1 relaxation time and the heat capacity were also observed to exhibit power law behavior consistent with point nodes in the gap structure [2]. Furthermore, a re-entrant superconducting phase was shown to develop in extremely high magnetic fields [4]. Intriguingly, heat capacity measurements in the superconducting state have consistently exhibited an apparent residual fermionic ( T -linear) term that appears to equal ∼ 50% of the normal state electronic DOS [2, 3], which has been interpreted as a signature of an unpaired fluid. However, this picture is problematic since a non-

Here we use thermal conductivity and penetration depth measurements to probe the superconducting gap structure directly by observing the properties of lowenergy quasiparticle excitations. Thermal conductivity κ has been established as a powerful directional tool for probing the gap structure of unconventional superconductors [8], and magnetic penetration depth is a well established tool for finding detailed information about the pairing state [9]. In our thermal transport measurements, we observe a residual quasiparticle contribution which is vanishingly small in the zero field limit, but is restored by magnetic field much more quickly than would be expected for either clean or dirty limit fully gapped superconductivity. Together with a quadratic temperature dependence of the measured magnetic penetration depth, our observations are consistent with a point node superconducting gap structure in UTe . We also present 2 measurements of the specific heat of the same samples used for thermal transport in order to probe the reported residual DOS in the superconducting state. We observe a low-temperature upturn in the T -linear term that is not consistent with a nuclear Schottky contribution, but

TABLE I. Candidate superconducting pairing states for the D 2 h point group symmetry.

| Representation   | Representative Function                 | Nodes   |
|------------------|-----------------------------------------|---------|
| A 1 g            | const+ c x kx 2 + c y k 2 y + c z k 2 z | none    |
| B 1 g            | ck z k y                                | lines   |
| B 2 g            | ck z k x                                | lines   |
| B 3 g            | ck x k y                                | lines   |
| A 1 u            | c x k x ˆ x + c y k y ˆ y + c z k z ˆ z | none    |
| B 1 u            | c 1 k y ˆ z + c 2 k x ˆ y               | points  |
| B 2 u            | c 1 k x ˆ z + c 2 k z ˆ x               | points  |
| B 3 u            | c 1 k x ˆ y + c 2 k y ˆ x               | points  |

rather a divergent term consistent with the presence of quantum critical degrees of freedom. Upon subtraction of the divergent term, we recover a complete entropy balance between superconducting and normal states in the electronic specific heat and the previously observed T 3 dependence [2, 3].

Single crystals of UTe 2 were synthesized using a vapor transport method as described in [2], using depleted uranium ( &gt; 99.9% 238 U) and natural abundance elemental 127 60 . Te (99.9999% purity). Crystals were shaped into rectangular plates with sample dimensions 1.5 × 0.5 × 0.1 mm 3 for the current j ‖ a sample (S1) and 1.8 × 1.0 × 0.3 mm 3 for the j ‖ b sample (S2), used for both thermal conductivity and specific heat measurements. Thermal conductivity κ was measured in a dilution refrigerator using a modular one-heater, two-thermometer probe [10], with contacts made by soldering gold wire leads with an In-Sn alloy to evaporated gold pads, resulting in ∼ mΩ contact resistances. Electrical resistivity ρ was measured in-situ using the same wires and contacts, with samples in the exact position used for thermal transport measurements. Both samples exhibit a residual resistance ratio ρ (300 K)/ ρ (0 K) of 22. As discussed below, the observation of the Wiedemann-Franz (WF) law in the field-induced normal state confirms the quality of measurements and rules out electron-phonon decoupling issues [11]. Measurements of the low temperature magnetic penetration depth ( λ ) were performed using a dielectric resonator technique [12] on two other crystals, labeled S3 and S4, and heat capacity C of samples S1 and S2 was measured using the relaxation time method.

The superconducting state in UTe 2 has been proposed [2-4] to be spin-triplet based on several experimental observations. Theoretically, because UTe 2 has an orthorhombic lattice that is symmorphic as confirmed at low temperatures by neutron scattering [13], Blount's theorem must hold [14] so that odd-parity ( i.e. , spintriplet) states will either have point nodes or be fully gapped. Conversely, even-parity (spin-singlet) states will have either line nodes or be fully gapped. More formally, assuming strong spin-orbit coupling and invoking the D 2 H point group yields the possible gap symmetries shown in Table 1. The possible symmetry-imposed nodal

FIG. 1. Thermal conductivity of UTe 2 samples S1 and S2. (a) Measured thermal conductivity in comparison to the estimated (elastic) charge contribution L /ρ 0 (solid and dashed lines) determined using the Wiedemann-Franz law and in-situ measured electrical resistivities of both samples. L /ρ T 0 ( ) is extrapolated below T c here assuming the resistivity follows the Fermi liquid temperature dependence ρ = ρ 0 + AT 2 . Panels (b) and (c) present signatures of the superconducting transition in thermal conductivity measurements for (b) j ‖ a and (c) j ‖ b orientations. The observed increase in κ/T below T c is thought to arise from a decrease in scattering due to the opening of the superconducting gap, and the suppression with increasing magnetic fields due to vortex scattering.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01069v2-figure-88cee506-b58a-4719-85f3-0512b374b584.jpg)

structures fall into three types: fully gapped ( A 1 g , A 1 u ), axial point nodes ( B 1 u , B 2 u , B 3 u with pairs of point nodes along a single high-symmetry axis), or polar line nodes ( B 1 g , B 2 g , B 3 g with lines node perpendicular to a single high symmetry axis).

Thermal conductivity is a directional probe that can be used to determine the momentum space position of nodes by studying low-energy quasiparticle excitations as a function of field angle, as shown in the cuprate superconductor YBCO [15, 16] and heavy fermion compound CeCoIn 5 [17], or as a function of the directional current as shown in UPt 3 [18, 19] and CeIrIn 5 [20]. Furthermore, the magnetic field evolution of the electronic thermal conductivity in the T = 0 limit ( κ 0 ( H )) is a reliable probe of both nodal gaps (e.g., UPt [18] and KFe As 3 2 2 [21]) and

multi-band or anisotropic gap superconductors [22].

Fig. 1(a) presents the total thermal conductivity of UTe 2 samples S1 and S2 measured up to 20 K and compared to the estimated (elastic) charge contribution determined using the in-situ measured resistivity ρ for each sample and the WF law ( i.e. , L /ρ 0 , where L 0 = π / 2 3( k B /e ) 2 ), showing that the majority of the thermal conduction in this range is not electronic. This makes distinguishing the quasiparticle contribution to thermal conductivity quite difficult, even within the superconducting regime. For both samples, there is an abrupt increase in κ/T at T c as highlighted in Figs. 1 (b) and (c), similar to that observed in CeCoIn 5 [23, 24] and YBCO [25, 26], where it is attributed to a suppression of the scattering of heat carriers (electrons and phonons, respectively) by the opening of the superconducting gap. In the case of UTe , the enhanced conductivity below 2 T c is likely due to a combination of electron and phonon carriers being affected by the gap opening, and, more important, is an indication of the high crystalline quality of the samples being measured. The rapid suppression of the enhancement in field is likely due to vortex scattering.

The nature of the superconducting gap structure can be determined by probing the value of the residual term of κ/T in the T = 0 limit, κ /T 0 , and its magnetic field dependence [27]. First, we confirm that the WF law is obeyed in the T =0 limit in the field-induced normal state at 7 T ( a -axis). To study the zero-field limit, we apply a small field of 50 mT to ensure soldered contacts are not superconducting and extrapolate to T = 0 by fitting the lowest temperature data to κ/T = κ /T 0 + AT α . As shown in Fig. 2 (a) and (b), we find that κ/T extrapolates to a vanishingly small residual value for both heat current orientations and rapidly rises with temperature. This is inconsistent with the behavior of a fully gapped superconductor [28], which would not conduct heat until at least ∼ 20% of T c , and is more consistent with a gap with nodal excitations [27].

To rule out the possibility of a strongly anisotropic but fully gapped scenario, we study the field dependence of κ /T 0 . For a full gap, quasiparticles remain localized to vortex cores and only contribute to the thermal conductivity for fields exponentially close to H c 2 , such as shown for the classic experiment on elemental Nb shown in Fig. 2(c) [29]. For a fully gapped superconductor in the dirty limit, κ /T 0 increases more quickly in field but still remains exponentially small in the low-field limit [30]). In multi-band superconductors with a disparity in gap values, a rapid rise can occur due to a very small field range H ∗ of activated behavior of a small gap ∆ 1 ( i.e. , where ∆ 1 / ∆ 2 ∼ √ H /H ∗ c 2 [22, 31]). For a nodal gap structure, low energy quasiparticle excitations are not bound and therefore follow a much stronger field dependence that depends on the specific gap structure (e.g. for line nodes, κ /T 0 ∝ √ H due to the Volovik effect).

Using the same extrapolation method as above, we plot in Fig. 2(c) the field dependence of κ 0 a /T and κ 0 b /T , which both exhibit a rapid quasi-linear rise with magnetic field, further confirmed by detailed 100 mK field-sweep measurements shown in the inset of Fig. 2(c). This behavior rules out a fully gapped scenario, including clean limit and dirty limit behaviors, and places a strong constraint on the presence of a very small gap (e.g. ∆ 1 glyph[lessmuch] ∆ ). Furthermore, we note that the temperature depen2 dence of κ/T in the presence of such a small gap should show a very rapid increase at low temperatures (or even an apparent residual term), making the extreme multiband scenario unlikely. Rather, the field dependence of κ /T 0 in UTe 2 is more reminiscent of that of nodal gap superconductors. Similar behavior was observed in line-node superconductors CeIrIn 5 [8], Sr 2 RuO 4 [32], and KFe As 2 2 [21], which are all accompanied by finite residual values in zero field. A line-node (e.g. d -wave) gap is expected to invoke a universal zero-field residual term given by κ /T 0 = 1 3 γ N F v 2 a glyph[planckover2pi1] 2 µ ∆ 0 . Using a conservative upper limit of κ /T 0 =0.1 mW/cm/K 2 from our experiment together with values 2∆ 0 = 3 5 . k B T c , γ N =100 mJ/mol/K 2 and µ = 2 for a standard d -wave gap slope at the nodes ( a is order unity) yields v F =6,000 m/s, a very small velocity similar to that of the extremely heavy p -wave superconductor UPt 3 [33], and likely not consistent with UTe . Interestingly, UPt 2 3 exhibits both vanishingly small κ 0 in zero field [33] and a linear increase with magnetic field [18], also shown for comparison in Fig. 2(c). Together with this striking similarity, the vanishingly small residual T -linear term and rapid quasi-linear rise in magnetic field observed for UTe 2 are most consistent with a point-node gap scenario.

Further information can be garnered from the anisotropy ratio ¯ κ ≡ κ /κ b a , shown in the inset of Fig. 2(d). Normalizing out the normal state anisotropy by using the ratio ¯(0 05 T) κ . /κ ¯(7 T), we find the superconducting state anisotropy approaches ∼ 0.8 in the lowT limit. This not as large as the value ( ∼ 0 5) observed . in UPt 3 [19], a more compelling indicator of point node direction. Rather, the near-unity anisotropy may arise due to nodal excitations directed in both basal plane directions (e.g. accidental), but will require detailed calculations based on candidate gap structures to draw firm conclusions.

Measurements of the low temperature magnetic penetration depth were performed using a cylindrical rutile dielectric resonator [34], which facilitates a microwave transmission resonance involving the ab -plane electrodynamic response at ≈ 11 GHz. The temperature dependence of the resonant frequency ∆ f 0 ( T ) and the quality factor Q T ( ) are measured and converted into surface impedance Z s through a geometrical factor estimated from the geometry of the resonator and the nature of the mode [35]. In the local electrodynamic limit, where the mean free path is smaller than the screening length, Z s can be converted to the complex conductivity

FIG. 2. Low temperature thermal conductivity of UTe 2 single crystal samples S1 (a) and S2 (b), measured with heat current applied along crystallographic a -and b -axes, respectively, and magnetic fields applied along the a -axis. The black × in (a) and (b) represent the converted normal state charge conductivity L /ρ 0 measured in-situ at 7 T and calculated using the Wiedemann-Franz law. (c) Magnetic field dependence of T → 0 K extrapolated values (, κ /T 0 ≡ κ T ( → 0) /T ) taken from the data in (a) and (b). Error bars correspond to uncertainty in the extrapolation from fitting to a power law over different temperature windows. For comparison we include data for elemental clean-limit fully gapped superconductor Nb [29] and cleanlimit nodal superconductors UPt 3 [18] and KFe As 2 2 [21]. The inset shows a comparison of the detailed low field dependence of ∆ κ H / ( ) ∆ κ N ≡ ( κ H ( ) -κ (0)) / κ ( (7 T) -κ (0)) measured in field-sweep mode at 100mK, with the T → 0 extrapolated values from the main panel. (d) Lorenz ratio L/L 0 calculated using the ratio of measured quantities κρ/T at 7 T from (a) and (b) to the Lorenz number L 0 . Linear extrapolations to T = 0 show agreement with L 0 to within 8% and 5% for j ‖ a and j ‖ b , respectively. The inset shows the zero-field anisotropy ratio ¯ = κ κ /κ b a normalized to the anisotropy ratio in the normal state ¯ κ N = ¯(7 T). κ

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01069v2-figure-29327336-0b57-402f-a1fb-850075b9694f.jpg)

˜ = σ σ 1 -iσ 2 of the samples through the local relation Z s = √ iµ ω/σ 0 ˜ [36]. The effective penetration depth is obtained from the imaginary part of the conductivity σ 2 = 1 /µ ωλ 0 2 eff ( T ), which yields λ eff (0) = 1126 nm for S3 and 947 nm for S4, both similar to the estimation ( λ (0) glyph[greaterorsimilar] 1000 nm) from a recent study [7].

As shown in Fig. 3, the change in screening length, ∆ λ eff ( T ), exhibits a rapid rise with temperature indicative of low-lying excitations [9]. Fitting ∆ λ eff ( T ) to a power law form aT n + b below 0 3 . T c yields a quadratic temperature exponent ( n = 2 00 . ± 0 01 for S3 and . n =

1 92 . ± 0 01 for S4) consistent with either a dirty limit line . node [37] or axial point node scenario [38, 39]. Considering the observed vanishingly small κ /T 0 values, and the deduced clean limit of the superconducting state ( meanfree-path of 95 nm estimated from the measured scattering rate Γ ≤ 0 005 . k B T c [34]) being much larger than the superconducting coherence length (4 ∼ 7 nm) deduced from H c 2 measurement [2]), we exclude the dirty limit scenario for this system and conclude that the quadratic dependence of ∆ λ eff ( T ) is more consistent with axial point nodes aligned in the direction of the vector poten-

FIG. 3. Effective penetration depth of UTe 2 samples S3 (red circles) and S4 (blue squares) plotted vs normalized temperature. The solid lines are power law fits of the form aT n + , b with extracted exponents n = 2 00 for S3 and . n = 1 92 for S4. . The vertical bars shown for each data set are the conservative uncertainty estimates. The inset presents the same data in terms of ( T/T c ) 2 to illustrate the quadratic temperature dependence, with linear dasned lines as guides.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01069v2-figure-18e20053-17bc-49d1-9d6b-5b58403d9867.jpg)

tial. In our setup, the applied microwave magnetic field generates circulating screening currents in the ab -plane [12], so that the direction of axial point nodes is assumed to lie in the ab -plane.

Concerning previous observations of an apparent residual DOS in specific heat [2, 3], our observations of a vanishingly small fermionic component of thermal conductivity in the zero-field, zero-temperature limit are inconsistent with a scenario that entails the presence of a sizeable residual non-paired fluid. To confirm this, Fig. 4 presents the measured specific heat of samples S1 and S2, which are the same crystals used for thermal transport. This measurement reproduces the large residual value in C/T observed in previous experiments [2, 3], but uncovers a sharp increase in C/T below about 300 mK. Although it appears similar to a nuclear Schottky anomaly, which arises from the splitting of nuclear spin states, such a contribution ( C N ) should be negligible since the only isotope with a nuclear spin is 125 Te, which is only 7% abundant in natural Te and does not have a nuclear quadrupole moment. Moreover, the observed upturn is much broader than the expected form C N /T ∝ T -3 . Instead, it is well fit by a weak power law ( C/T ∝ T -∼ 1 / 3 ) divergence as described below.

Motivated by this, we fit the total heat capacity over the entire temperature range ( i.e. , 50 mK to 5 K) using

$$\frac { C } { T } = A _ { p h } T ^ { 2 } + A _ { d i v } T ^ { - \alpha } + \begin{cases} A _ { s c } T ^ { 2 } & T < T _ { c } \\ \gamma _ { n } & T > T _ { c } \end{cases}. \quad ( 1 ) \quad \text{with} \\ \text{fermi} \\ \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Psi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi\Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phis \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi$$

The result is shown in Fig. 4(a), where the solid line rep-

FIG. 4. Analysis of the low temperature heat capacity of UTe 2 single-crystal samples S1 and S2 (same crystals used for thermal transport measurements). (a) The total measured heat capacity (open circles) of sample S1 decomposed into a weak power law divergence (dashed line), phonon term (dash-dotted line), and electronic heat capacity (filled circles) obtained by subtracting the diverging and phonon terms (see text for details). (b) The obtained electronic heat capacities of samples S1 and S2 plotted vs. T 2 . Solid lines are fits to the T 2 dependence in the superconducting state with T c determined by balancing entropy.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01069v2-figure-39129359-14be-41ab-a715-2424ea3d4cef.jpg)

resents the total fit, A ph T 2 is the phonon heat capacity (dash-dotted line), A div T -α is the power law divergence (dashed line) where α = 0.35 (0.33) for S1 (S2), and the last term is the electronic heat capacity which varies as T 2 in the superconducting regime ( T &lt; T c ) and is constant ( γ n ) above T c .

The extracted electronic component C /T el = C/T -A ph T 2 -A div T -α exhibits a nearly perfect entropy balance between normal and superconducting states, which is reflected by the fit result in Fig. 4(a) where the resultant T c (set by the entropy balance) matches the data very well. Fig. 4(b) shows that the specific heat for both S1 and S2 samples follows a cubic power law ( i.e. , C /T el ∼ T 2 ) through the superconducting temperature range, consistent with prior fits to raw data above 300 mK [2, 3]. While the divergent residual specific heat term requires further investigation, the lack of any sizeable contribution to zero-field thermal conduction means that this component arises from either a well-localized DOS, or itinerant fermionic carriers that are heavily scattered in the T → 0 limit.

In conclusion, detailed measurements of thermal conductivity, penetration depth and heat capacity in UTe 2 provide evidence for a superconducting gap structure with point nodes. The vanishingly small residual fermionic term in thermal conductivity, together with rapid increases with both temperature and magnetic field in the zero-field, zero-temperature limits, respectively,

confirm the presence of low-energy quasiparticle excitations and the proposed nodal gap structure, and rule out an anistropic or multi-band full gap scenario. A quadratic temperature dependence of the low temperature penetration depth is also consistent with an nodal gap structure, with point nodes in the ab -plane. Together with consideration of the D 2 H point group symmetry, these results provide evidence for spin-triplet Cooper pairing in UTe . 2 The origin of a localized or strongly scattered divergent quantum critical component of the specific heat, which underlies the superconducting state, will require further study to elucidate.

- ∗ paglione@umd.edu
- [1] D. Aoki et al. , J. Phys. Soc. Jpn. 88 , 022001 (2019).
- [2] S. Ran et al. , Science 365 , 684 (2019).
- [3] D. Aoki et al. , J. Phys. Soc. Jpn. 88 , 043702 (2019).
- [4] S. Ran et al. , Nat. Phys. (2019), doi:10.1038/s41567-0190670-x.
- [5] Q. Niu et al. , arXiv:1907.11118.
- [6] A. Miyake et al. , J. Phys. Soc. Jpn. 88 , 063706 (2019).
- [7] S. Sundar et al. , Phys. Rev. B 100 , 140502(R) (2019).
- [8] H. Shakeripour et al. , New J. Phys. 11 , 055065 (2009).
- [9] R. Prozorov and R.W. Giannetta, Supercond. Sci. Technol. 19 , R41 (2006).
- [10] M.A. Tanatar et al. , Rev. Sci. Instr. 89 , 013903 (2018).
- [11] M. Smith et al. , Phys. Rev. B 71 , 014506 (2005).
- [12] S. Bae et al. , Rev. Sci. Instrum. 90 , 043901 (2019).
- [13] V. Hutanu et al. , arXiv:1905.04377.
- [14] E.I. Blount, Phys. Rev. B 32 , 2935 (1985).
- [15] F. Yu et al. , Phys. Rev. Lett. 74 , 5136 (1995).
- [16] H. Aubin et al. , Phys. Rev. Lett. 78 , 2624 (1997).
- [17] K. Izawa et al. , Phys. Rev. Lett. 87 , 057002 (2001).
- [18] H. Suderow et al. , J. Low Temp. Phys. 108 , 11 (1997).
- [19] B. Lussier et al. , Phys. Rev. B 53 , 5145 (1996).
- [20] H. Shakeripour et al. , Phys. Rev. Lett. 99 , 187004 (2007).
- [21] J.-Ph. Reid et al. , Phys. Rev. Lett. 109 , 087001 (2012).
- [22] E. Boaknin et al. , Phys. Rev. Lett. 90 , 117003 (2003).
- [23] R. Mosovich et al. , Phys. Rev. Lett. 86 , 5152 (2001).
- [24] M. Tanatar et al. , Phys. Rev. Lett. 95 , 067002 (2005).
- [25] J.L. Cohn et al. , Phys. Rev. B 45 , 1314(R) (1992).
- [26] R.C. Yu et al. , Phys. Rev. Lett. 69 , 1431 (1992).
- [27] M. J. Graf et al ., Phys. Rev. B 53 , 15147 (1996).
- [28] J. Bardeen et al. , Phys. Rev. 113 , 982 (1959).
- [29] J. Lowell et al. , J. Low Temp. Phys. 3 , 65 (1970).
- [30] J.O. Willis and D.M. Ginsberg, Phys. Rev. B 14 , 1916 (1976).
- [31] M. Yamashita et al. , Phys. Rev. B 80 , 220509(R) (2009).
- [32] M. Tanatar et al. , Phys. Rev. B 63 , 064505 (2001).
- [33] R. Joynt and L. Taillefer, Rev. Mod. Phys. 74 , 235 (2002).
- [34] S. Bae et al. , arXiv:1909.09032.
- [35] M. A. Hein, High-Temperature Superconductor Thin Films at Microwave Frequencies (Springer, Heidelberg, 1999) pp. 45-46.
- [36] M. Dressel and G. Gruner, Electrodynamics of Solids (Cambridge University Press, 2002).
- [37] P. J. Hirschfeld and N. Goldenfeld, Phys. Rev. B 48 , 4219 (1993).
- [38] D. Einzel et al. , Phys. Rev. Lett. 56 , 2413 (1986).
- [39] F. Gross et al. , Z. Phy. B Cond. Matt. 64 , 175 (1986).
Tristin Metz Seokjin Bae Sheng Ran I-Lin Liu Yun Suk Eo Wesley T. Fuhrman Daniel F. Agterberg Steven Anlage Nicholas P. Butch Johnpierre Paglione
Point Node Gap Structure of Spin-Triplet Superconductor UTe2
Low-temperature electrical and thermal transport, magnetic penetration depth,
and heat capacity measurements were performed on single crystals of the
actinide superconductor UTe2 to determine the structure of the superconducting
energy gap. Heat transport measurements performed with currents directed along
both crystallographic a- and b-axes reveal a vanishingly small residual
fermionic component of the thermal conductivity. The magnetic field dependence
of the residual term follows a rapid, quasi-linear increase consistent with the
presence of nodal quasiparticles, rising toward the a-axis upper critical field
where the Wiedemann-Franz law is recovered. Together with a quadratic
temperature dependence of the magnetic penetration depth up to T/T_c=0.3, these
measurements provide evidence for an unconventional spin-triplet
superconducting order parameter with point nodes. Millikelvin specific heat
measurements performed on the same crystals used for thermal transport reveal
an upturn below 300 mK that is well described by a divergent quantum-critical
contribution to the density of states (DOS). Modeling this contribution with a
T^{-1/3} power law allows restoration of the full entropy balance in the
superconducting state and a resultant cubic power law for the electronic DOS
below T_c, consistent with the point-node gap structure determined by thermal
conductivity and penetration depth measurements.