Evidence for even parity unconventional superconductivity in Sr2RuO4 A. Chronister1†,∗ A. Pustogow1†,∗ N. Kikugawa2, D. A. Sokolov3, F. Jerzembeck3, C. W. Hicks3, A. P. Mackenzie3,4, E. D. Bauer5, and S. E. Brown1† 1Department of Physics & Astronomy, UCLA, Los Angeles, CA 90095, USA; 2National Institute for Materials Science, Tsukuba 305-0003, Japan; 3Max Planck Institute for Chemical Physics of Solids, Dresden 01187, Germany; 4Scottish Universities Physics Alliance, School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews KY16 9SS, UK; and 5Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA; (Dated: July 28, 2020) Unambiguous identification of the superconducting order parameter symmetry of Sr2RuO4 has remained elusive for more than a quarter century. While a chiral p-wave ground state analogue to superfluid 3He-A was ruled out only very recently, other proposed p-wave scenarios are still viable. Here, field-dependent 17O Knight shift measurements are compared to corresponding specific heat measurements, previously reported. We conclude that the shift results can be accounted for by the expected field-induced quasiparticle response only. An upper bound for the condensate magnetic response of < 10% of the normal state susceptibility is sufficient to exclude odd-parity candidates. Unraveling the secrets of the superconducting state in Knight shifts Ks(T <Tc, B) are a crucial probe of the Sr2RuO4 [1–3] has been a priority for unconventional su- order-parameter symmetry, in particular for Sr2RuO4. perconductivity research since its discovery in 1994 by In the normal state, Ks χn, with χn the susceptibility. Maeno and coworkers [4]. Particularly notable among In the superconducting∼ phase, a nonzero susceptibility several reasons for broad interest in Sr2RuO4 was the χsc associated with condensate polarization is expected suggestion of a p-wave triplet pairing state [5]. Among generally for triplet-paired, p-wave states. Hence, the states allowed by symmetry is the chiral state z(px ipy), observed reduction of the Knight shift for an applied in- which breaks time reversal symmetry and therefore± re- plane field excludes the chiral state [11], for which d c. quires two components. Soon after, the combination of Not eliminated are states characterized by d c. k + ⊥ results from NMR Knight shift [6] and µ SR [7] mea- At first sight, the most direct way to test for symmetry- surements lent support to the chiral p-wave description. allowed states with d c is to perform measurements Further evidence was inferred from the observed onset of with B c. However,⊥ the relevant upper critical field a non-zero Kerr rotation at Tc [8]. Unresolved issues re- k Bc2,[001] < 1 kG is very small [21] making such exper- mained, however. For example, thermal conductivity [9] iments particularly challenging because signal strength and specific heat [10] experiments were both interpreted and spectral resolution are reduced for very weak ap- as evidence for a nodal gap structure [3]. In a step toward 17 plied fields. Here, we take another approach, discussed clarification, recent O NMR measurements excluded previously in Refs. [12, 22]: the field orientation is fixed any candidate p-wave state with a d-vector aligned par- 17 in-plane, and the O shifts Ks are evaluated at low tem- allel to the c-axis [11, 12]. The NMR results published perature (25 mK) while varying B as much as experimen- so far leave open the possibility for states with an in- tally feasible. Quasiparticle creation is controlled by the plane d, as has since been explicitly discussed in several field strength, and also contributes to the magnetic re- theoretical works [13–15]. sponse. By way of comparing to previously reported spe- With these developments in mind, we recall other phys- cific heat results Ce(B)/T [10, 23], we estimate the up- ical properties of Sr2RuO4 deemed relevant to the order- per bound for the condensate portion to χsc/χn < 10%, parameter symmetry. Recognized early [4, 5] was the a value that contradicts the expectation for any of the evidence for ferromagnetic correlations, as originally in- pure p-wave order parameters relevant for Sr2RuO4. arXiv:2007.13730v1 [cond-mat.supr-con] 27 Jul 2020 ferred from the measured Wilson ratio exceeding unity. As in previous NMR studies on Sr2RuO4 [6], the la- Meanwhile, quantum oscillations [16], ARPES [17, 18], belled 17O (17I=5/2, 17γ=-5.772 MHz/T [24]) is intro- and optical conductivity experiments [19] all indicate duced by high-temperature annealing [6], here in 90% modest normal-state mass enhancement. At the same 17 ◦ O2 atmosphere at 1050 C. To facilitate access to time, an unconventional ground state was indicated by relatively low frequencies covering several octaves, we its very strong suppression with non-magnetic impuri- adopted a top tuning/matching configuration. Single- ties [20]. Thus, Sr2RuO4 constitutes a model system crystal dimensions were (3.5 mm x 1 mm x 0.2 mm), for the emergence of an unconventional superconducting with the shortest dimension corresponding to the out-of- state at low temperature, from a strongly correlated and plane [001]-direction, and the longest dimension parallel ultraclean Fermi-liquid normal state. to [100], see Fig. 1(a). The NMR coil containing the The temperature and field dependences of the NMR crystal under study, was mounted on a single-axis piezo- 2 rotator inside the mixing chamber of a bottom-loading [001] (a) O(2) (b) O(2) T = 25 mK dilution refrigerator. Sample alignment enabled in-plane [010] ◦ B = 1.50 T orientation to within 0.2 , based on RF susceptibility O(1 ) ± ⊥ O(1||) Ru O(1 ) measurements sensitive to Bc2, described in Ref. [11], and O(1||) ⊥ NMRIntensity units) (arb. discussed in the Supplemental Material [25]. 63Cu NMR [100] O(1 ) O(1 ) || relaxation rate measurements were used to determine the ⊥ Ep = 40 nJ equilibrium bath temperature T = 25 mK. As in our pre- B vious work [11], low-power RF experiments were carried O(2) 4 nJ out to make sure the results were not measurably altered 1.0 ) (d) (c) heating by B = 1.38 T by RF pulse heating effects. The applied field strength =0 s K 0.8 NMR pulses f normal B was determined to within uncertainties less than 10’s - normal 3 3 4 of µT from the NMR resonance of He in the He/ He 0.6 normal 130 nJ mixture of the dilution refrigerator. f )/( 0.4 0500 1000 SC SC =0 Addressed first are sample heating effects by the RF s B (T) K 0.8 f 40 nJ p - b⋅E 1.14 pulses, illustrated in Fig. 1, which turned out to be a ⊥ a+ 1 0.2 0.6 f 0.94 crucial issue [11, 12]. So as to enhance sensitivity to ( 0.4 0.73 0.4 nJ this potential artifact, we examined the transients with 0.0 1 10 100 1000 -40 -20 0 20 40 the field set to 1.38 T, a value just smaller than B . c2 E (nJ) f-f (kHz) Clear evidence for warming by the RF pulsing is in- p 0 ferred from a transient response corresponding to that FIG. 1. (a) Sr2RuO4 involves three distinct oxygen sites for of the normal-state (instead of the sought-after super- field direction B k [100]. (b) The three associated 17O NMR conducting state). Shown in Fig. 1(b,c) are 17O spectra central transitions (O(1k), O(2), O(1⊥) from left to right) are corresponding to central transitions for the three sites, independent of pulse energy Ep at 1.50 T > Bc2 ≃ 1.45 T. O(1k,2,1⊥), at applied magnetic fields slightly above and (c) Also at B = 1.38T . Bc2 the normal-state spectrum is −7 below Bc2. While at 1.5 T >Bc2 the line shape remains observed for Ep ≥ 10 J. Reducing to Ep = 40 nJ leads unaffected by changing the pulse energy, a normal state to doubled spectral features, most pronounced for O(1k,⊥), spectrum is produced also at 1.38 T < B when using which we assign to coexisting normal (dashed vertical lines) c2 and superconducting (solid) contributions around the first- a pulse energy Ep = 130 nJ. Decreasing Ep to 40 nJ order transition. Further reduction of Ep reveals the pure leads to a response where a new spectral line appears for superconducting-state spectrum. (d) O(1⊥) frequencies nor- each site, indicating the coexistence of normal and super- malized to normal-state (fnormal) and zero-shift (fKs=0; see conducting phases. This data set is particularly useful, Fig. 2) positions at B < Bc2 for variable Ep. Linear fits since the macroscopic phase segregation provides a quan- (solid lines, see inset) indicate that heating is less problem- titative measure of the magnetization jump ∆M at the atic at lower field due to larger Tc(B). Knight shifts Ks were discontinuous (first-order) transition [26, 27]. Note that determined using the frequency values leveling off at Ep → 0. these data are recorded following a single-pulse excita- tion. That is, the transient NMR response corresponds to a free induction decay (FID). All shift results of the are transition frequencies at each field, generated using present work were obtained from FID measurements car- the known normal state NMR parameters [25]. Then, the ried out with RF pulse energies low enough to avoid heat- frequency differences between closed and open symbols ing, as illustrated in Fig. 1(d). are proportional to the hyperfine fields, and constitute Having established a threshold for heating effects, the product of (normal-state) Knight shifts with applied 17 we now inspect the spectra recorded at variable field field, Ks,normal γB, for O(1k), O(2) and O(1⊥).