Evidence for chiral d-wave superconductivity in URu2Si2 from the field-angle variation of its specific heat\abst
Low-energy quasiparticle (QP) excitations in the heavy-fermion superconductor URuSi were investigated by specific-heat measurements of a high-quality single crystal. The occurrence of QP excitations due to the Doppler-shift effect was detected regardless of the field direction in of the present clean sample, which is in sharp contrast to a previous report. Furthermore, the polar-angle-dependent measured under a rotating magnetic field within the plane exhibits a shoulder-like anomaly at and a sharp dip at () in the moderate-field region. These features are supported by theoretical analyses based on microscopic calculations assuming the gap symmetry of , whose gap structure is characterized by a combination of a horizontal line node at the equator and point nodes at the poles. The present results have settled the previous controversy over the gap structure of URuSi and have authenticated its chiral -wave superconductivity.
The pairing mechanism of unconventional superconductivity in the heavy-fermion compound URuSi ( K) has been one of the most tantalizing questions in condensed-matter physics. It coexists with the “hidden-order” (HO) phase occurring below K, a mysterious phase that exhibits neither an ordinary magnetic order nor a structure change . The true nature of the HO phase is a long-standing mystery and has been intensively studied by a number of experimental and theoretical groups. High-pressure experiments have revealed that superconductivity disappears and the HO phase transforms to the antiferromagnetic phase with a large moment at a critical pressure of GPa . The Shubnikov-de Haas experiments performed under pressure  indicate that the Fermi surface is similar between the HO and antiferromagnetic phases. On the basis of a first-principles theoretical approach, emergent rank-5 nematic order has been proposed as a possible origin of the HO phase , which can naturally explain the rotational symmetry breaking suggested from recent experiments [5, 6, 7, 8, 9].
Towards the identification of the pairing mechanism of superconductivity, the determination of the gap symmetry is essential because it is closely related to the pairing glue. The spin part of the order parameter for URuSi has been considered to be singlet because the Pauli-paramagnetic effect is observed in any field direction , although a change in the NMR Knight shift below has not yet been detected . Recent polar Kerr effect measurements  and other experiments [13, 14] have provided evidence for the spontaneous breaking of time-reversal symmetry in the superconducting state, i.e., a chiral state in which two order parameters are degenerated. Indeed, a giant Nernst effect has recently been observed above , and has been interpreted by asymmetric scatterings due to chiral superconducting fluctuations . It has been widely accepted that the superconducting gap possesses nodes on the basis of the power-law temperature dependence of thermodynamic quantities [17, 18, 19]. The invariance of the specific heat  and the thermal conductivity [19, 21] with changing in-plane field direction demonstrates the absence of vertical line nodes. From these experimental findings, chiral -wave superconductivity, the pairing state, has been inferred. For this state, the superconducting gap is expected to have a horizontal line node at and point nodes at .
Nevertheless, the detection of the detailed nodal structure of URuSi remains controversial. Previous measurements of the specific heat performed at 0.34 K  indicated the presence of point nodes at the north and south poles but did not detect the horizontal line node in the heavy-mass bands. This conclusion was derived from the observation that the specific heat showed -linear (nearly -linear) dependence in the low-field region for (). On the basis of the Doppler-shift effect , , the results suggest that Doppler-shift-active quasiparticles (QPs) are absent for , i.e., all the nodal QPs have their parallel to the direction. Here, is the Fermi velocity at a node and is the supercurrent velocity circulating around each vortex (). Meanwhile, thermal-conductivity studies [19, 21] concluded a horizontal line node in the light-mass bands and point nodes in the heavy-mass bands. Thus, experimental evidence for a horizontal line node in the heavy-mass bands has been missing.
Regarding the ordered-state Fermi surface of URuSi, recent experimental and theoretical studies indicate the presence of heavy-mass sheets around the M point [4, 23, 5]. Therefore, the absence of a horizontal line node on the heavy-mass band is inconsistent with the anticipated chiral -wave state. In order to settle this controversy, in this study we measured the field-angle-dependent specific heat down to a low temperature of 0.2 K. We used a high-quality single-crystalline sample (10.6 mg weight) grown by the Czochralski pulling method in a tetra-arc furnace . The residual resistivity ratio of another crystal from the same batch is 255. The specific heat was measured by the standard quasi-adiabatic heat-pulse and relaxation methods in a dilution refrigerator. The magnetic-field orientation was controlled with high accuracy using a vector magnet.
Figures 1(a) and 1(b) show the temperature dependence of measured under a magnetic field parallel to the and axes, respectively. The zero-field shows a sharp jump at of 1.41 K and takes relatively small values at low temperatures below 0.2 K. These results ensure that the present sample is of higher quality than the sample used in the previous study .
In the intermediate- region between 0.3 and 1 K, the zero-field exhibits a nearly -linear dependence as previously observed. The extrapolation of this -linear behavior to , which yielded a positive intersect in previous reports [17, 24, 20], becomes negative for the present clean sample. Apparently, a simple dependence does not hold in the whole temperature range below 1 K and the slope of becomes intrinsically smaller below 0.3 K. These features, along with a small jump of the specific heat at , , can be reproduced approximately from a single-band calculation by assuming the gap symmetry of with a gap size of , as represented by the dashed line in Fig. 1(b). Incorporating the multiband effect would improve the agreement between the calculated result and the experimental data.
The slight upturn in on cooling below 0.2 K is probably due to some impurity effect, because it becomes larger in lower-quality samples . With increasing the magnetic field, this upturn significantly develops in any field direction. The field-enhanced upturn in is too large to be explained by the nuclear Schottky contribution. To exclude the effect of this unusual upturn, we hereafter concentrate on the data measured above 0.2 K.
For , the sample is always in the normal state at 3 T. As shown in Fig. 1(b), the normal-state decreases slightly on cooling to 0.2 K. Similar behavior has also been observed in previous results .
Figures 1(c) and 1(d) represent the field dependence of at 0.2, 0.34, 0.6, and 1.6 K for and , respectively. The data measured at 1.6 K show a slight change with , reflecting its normal-state property. Because of the limited field range of the present experiment, the upper critical field ( T) was not reached for .
Let us first focus on the data at 0.34 K, taken under the same conditions as those described in Ref. References. In sharp contrast to the previous report, in which no behavior was observed for , for the present sample shows -like behavior not only for but also for in the low-field region. This change probably occurred owing to the improvement of the sample quality; in general, the low-energy QP excitations are obscured easily by the impurity-scattering effect  and thermal QP excitations. Indeed, by increasing the temperature to 0.6 K, the behavior becomes indiscernible, particularly for .
To more clearly see the low-energy QP excitations in URuSi, we measured at a lower temperature of 0.2 K. The rapid increase in at low fields is more clearly detected for both field directions. These low-field data can be fitted satisfactorily by using the function as demonstrated by dashed lines in Figs. 1(c) and 1(d). Furthermore, we have confirmed that the present and data measured for in T obey the scaling law  in a wide region, as shown in the inset of Fig. 1(a). Here, and is a scaling function that is constant (proportional to ) for small (large) . From these results, it is concluded that the heavy-mass bands have node(s) in regions where . In other words, even the existence of point nodes at poles is not justified by the data alone.
Then, we need to reexamine the gap structure of URuSi. Figure 2 shows the azimuthal-angle dependence of the specific heat measured by rotating a magnetic field within the plane at various field strengths and temperatures. Here, the field angle is measured from the axis. Recall that, according to Doppler-shift analyses [27, 28], the total QP density of states shows a local minimum when the magnetic field is along the nodal or gap-minimum direction. At K, the specific heat is invariant with at any field, consistent with the previous report . By contrast, at 0.2 and 0.6 K, we observed a clear fourfold oscillation in at 1 T (). Moreover, the sign of the oscillation changes at around 0.34 K, reminiscent of the behavior in CeCoIn . However, this fourfold oscillation is diminished by decreasing the magnetic field to 0.5 T (). This behavior, which is similar to the case of the anisotropic -wave superconductor CeRu , can be explained by the presence of gap minima and/or anisotropy of the Fermi velocity [31, 32]. Thus, we conclude that vertical line nodes for the rotation are absent.
Next, to study the horizontal line nodes and polar point nodes, we investigated the polar-angle dependence of the specific heat by rotating a magnetic field within the plane. Here, is the angle between the magnetic field and the axis. The data measured at 0.2 K are shown in Fig. 3. Note that the sample is in the normal state near at 3 T because T [see Fig. 1(d)]. In the high-field region, a large twofold oscillation, reflecting the anisotropy of that is predominantly due to the anisotropic Pauli-paramagnetic effect, is observed.
In the low-field region below 0.15 T, no significant anomaly was detected in within the sensitivity limits of our measurements: the data can be fitted by a simple twofold function, . By contrast, a striking feature is found in at 0.2 T; a shoulder-like anomaly appears at around . With increasing the magnetic field, this anomaly slightly moves to the higher side but is visible up to at least 1.5 T. In addition, the dip at becomes sharp and pronounced in the intermediate-field region between 0.3 and 1.5 T.
In order to examine whether these features of can be explained by the anticipated chiral -wave state, we calculate the polar-angle dependence of the zero-energy density of states, , on the basis of the microscopic Eilenberger theory assuming the gap symmetry of [inset of Fig. 4(a)] and a single-band spherical Fermi surface . For simplicity, the Pauli-paramagnetic effect and the anisotropy of the Fermi velocity are not taken into account. Figure 4(a) shows the field dependence of at () and (). Both data show behavior, as observed in the present experiment. At low fields, is slightly larger for than for . By increasing the magnetic field, this anisotropy is reversed at because of the slight anisotropy of () originating from the gap structure.
In Fig. 4(b), is plotted as a function of the polar angle at several fields. The calculated result at a low field of can reproduce the features observed in the present experiment, i.e., the shoulder-like anomaly and the dip structure in observed at moderate fields. With increasing the magnetic field, a sharp dip develops at and then the anisotropy in is reversed above . This reflects the anisotropy, as mentioned above. The dip at and the reversal of the anisotropy in the high-field region, demonstrated in Fig. 4(b), were not detected in the experiment because of the actual anisotropy owing to the strong Pauli-paramagnetic effect.
To determine the reason why the shoulder-like anomaly appears, the -resolved density of states, , obtained from the present calculations is mapped on the spherical Fermi surface in Fig. 4(c). Here, the magnetic field of is applied within the plane. Note that the integration of over the Fermi surface yields . When the magnetic field is applied along the axis (), prominent QP excitations occur anywhere near the horizontal line node while those at the point nodes are inactive. By tilting the magnetic field away from the axis, QPs are also excited around the point nodes.
An important change in can be seen between and ; for , the contribution from the part of the horizontal line node near is strongly reduced because the Doppler-shift effect becomes small there. As a result, a substantial suppression of the total QP density of states occurs for and leads to the shoulder-like anomaly around at . This dependence is less sensitive to the change in the Fermi-surface shape; qualitatively the same results can be obtained from calculations by using cigar-shaped and spheroidal Fermi surfaces.  Hence, the observed shoulder-like anomaly in can be attributed to the existence of a horizontal line node at . Note that, in the even-parity basis functions for tetragonal-lattice symmetry, the presence of a horizontal line node at and the absence of vertical line nodes only match with the chiral symmetry . Therefore, the present results strongly indicate that the gap symmetry of URuSi is of the type.
In summary, the field-orientation-dependent specific heat of URuSi was measured by using a high-quality single crystal. behavior of the specific heat was observed in any field direction for the present clean sample, whereas it was not observed for in the previous study . A fourfold oscillation of the specific heat was found in a rotating magnetic field within the plane, but it disappeared by lowering the magnetic field. This fact ensures the absence of vertical line nodes in the superconducting gap. In addition, we found a shoulder-like anomaly at and a sharp dip at in the polar-angle dependence of the specific heat. From theoretical analyses, these features can be explained by assuming the presence of a horizontal line node at the equator and point nodes at the poles in the gap. Thus, the gap symmetry of URuSi has been identified as a chiral -wave type described by .
Acknowledgements.We thank Y. Yanase for useful discussions. A part of the numerical calculations was performed by using the HOKUSAI GreatWave supercomputer system in RIKEN. This work was supported by a Grant-in-Aid for Scientific Research on Innovative Areas “J-Physics” (15H05883) from MEXT, and KAKENHI (15K05158, 15H03682, 25103716, 26400360, 15H05745, 15K17715, 15J05698) from JSPS.
- J. A. Mydosh and P. M. Oppeneer, Rev. Mod. Phys. 83, 1301 (2011).
- H. Amitsuka, K. Matsuda, I. Kawasaki, K. Tenya, M. Yokoyama, C. Sekine, N. Tateiwa, T. C. Kobayashi, S. Kawarazaki, and H. Yoshizawa, J. Magn. Magn. Mater. 310, 214 (2007).
- E. Hassinger, G. Knebel, T. D. Matsuda, D. Aoki, V. Taufour, and J. Flouquet, Phys. Rev. Lett. 105, 216409 (2010).
- H. Ikeda, M.-T. Suzuki, R. Arita, T. Takimoto, T. Shibauchi, and Y. Matsuda, Nat. Phys. 8, 528 (2012).
- T. Shibauchi, H. Ikeda, and Y. Matsuda, Philos. Mag. 94, 3747 (2014).
- R. Okazaki, T. Shibauchi, H. J. Shi, Y. Haga, T. D. Matsuda, E. Yamamoto, Y. Onuki, H. Ikeda, and Y. Matsuda, Science 331, 439 (2011).
- S. Tonegawa, K. Hashimoto, K. Ikada, Y.-H. Lin, H. Shishido, Y. Haga, T. D. Matsuda, E. Yamamoto, Y. Onuki, H. Ikeda, Y. Matsuda, and T. Shibauchi, Phys. Rev. Lett. 109, 036401 (2012).
- S. Tonegawa, S. Kasahara, T. Fukuda, K. Sugimoto, N. Yasuda, Y. Tsuruhara, D. Watanabe, Y. Mizukami, Y. Haga, T. D. Matsuda, E. Yamamoto, Y. Onuki, H. Ikeda, Y. Matsuda, and T. Shibauchi, Nat. Commun. 5, 4188 (2014).
- S. Kambe, Y. Tokunaga, H. Sakai, T. D. Matsuda, Y. Haga, Z. Fisk, and R. E. Walstedt, Phys. Rev. Lett. 110, 246406 (2013).
- J. P. Brison, N. Keller, A. Vernire, P. Lejay, L. Schmidt, A. Buzdin, J. Flouquet, S. R. Julian, and G. G. Lonzarich, Physica C 250, 128 (1995).
- Y. Kohori, K. Matsuda, and T. Kohara, J. Phys. Soc. Jpn. 65, 1083 (1996).
- E. R. Schemm, R. E. Baumbach, P. H. Tobash, F. Ronning, E. D. Bauer, and A. Kapitulnik, Phys. Rev. B 91, 140506(R) (2015).
- G. Li, Q. Zhang, D. Rhodes, B. Zeng, P. Goswami, R. E. Baumbach, P. H. Tobash, F. Ronning, J. D. Thompson, E. D. Bauer, and L. Balicas, Phys. Rev. B 88, 134517 (2013).
- I. Kawasaki, I. Watanabe, A. Hillier, and D. Aoki, J. Phys. Soc. Jpn. 83, 094720 (2014).
- T. Yamashita, Y. Shimoyama, Y. Haga, T. D. Matsuda, E. Yamamoto, Y. Onuki, H. Sumiyoshi, S. Fujimoto, A. Levchenko, T. Shibauchi, and Y. Matsuda, Nat. Phys. 11, 17 (2015).
- H. Sumiyoshi and S. Fujimoto, Phys. Rev. B 90, 184518 (2014).
- R. A. Fisher, S. Kim, Y. Wu, N. E. Phillips, M. Mcelfresh, M. S. Torikachili, and M. B. Maple, Physica B 163, 419 (1990).
- K. Matsuda, Y. Kohori, and T. Kohara, J. Phys. Soc. Jpn. 65, 679 (1996).
- Y. Kasahara, T. Iwasawa, H. Shishido, T. Shibauchi, K. Behnia, Y. Haga, T. D. Matsuda, Y. Onuki, M. Sigrist, and Y. Matsuda, Phys. Rev. Lett. 99, 116402 (2007).
- K. Yano, T. Sakakibara, T. Tayama, M. Yokoyama, H. Amitsuka, Y. Homma, P. Miranović, M. Ichioka, Y. Tsutsumi, and K. Machida, Phys. Rev. Lett. 100, 017004 (2008).
- Y. Kasahara, H. Shishido, T. Shibauchi, Y. Haga, T. D. Matsuda, Y. Onuki, and Y. Matsuda, New J. Phys. 11, 055061 (2009).
- G. E. Volovik, JETP Lett. 58, 469 (1993).
- S. Tonegawa, K. Hashimoto, K. Ikada, Y. Tsuruhara, Y.-H. Lin, H. Shishido, Y. Haga, T. D. Matsuda, E. Yamamoto, Y. Onuki, H. Ikeda, Y. Matsuda, and T. Shibauchi, Phys. Rev. B 88, 245131 (2013).
- T. D. Matsuda, E. Hassinger, D. Aoki, V. Taufour, G. Knebel, N. Tateiwa, E. Yamamoto, Y. Haga, Y. nuki, Z. Fisk, and J. Flouquet, J. Phys. Soc. Jpn. 80, 114710 (2011).
- C. Kübert and P. J. Hirschfeld, Solid State Commun. 105, 459 (1998).
- G. E. Volovik and N. B. Kopnin, Phys. Rev. Lett. 78, 5028 (1997).
- I. Vekhter, P. J. Hirschfeld, J. P. Carbotte, and E. J. Nicol, Phys. Rev. B 59, R9023 (1999).
- T. Sakakibara, A. Yamada, J. Custers, K. Yano, T. Tayama, H. Aoki, and K. Machida, J. Phys. Soc. Jpn. 76, 051004 (2007).
- K. An, T. Sakakibara, R. Settai, Y. nuki, M. Hiragi, M. Ichioka, and K. Machida, Phys. Rev. Lett. 104, 037002 (2010).
- S. Kittaka, T. Sakakibara, M. Hedo, Y. nuki, and K. Machida, J. Phys. Soc. Jpn. 82, 123706 (2013).
- P. Miranovi, N. Nakai, M. Ichioka, and K. Machida, Phys. Rev. B 68, 052501 (2003).
- P. Miranovi, M. Ichioka, K. Machida, and N. Nakai, J. Phys.: Condens. Matter 17, 7971 (2005).
- Y. Tsutsumi, T. Nomoto, H. Ikeda, and K. Machida, unpublished.
- M. Sigrist and K. Ueda, Rev. Mod. Phys. 63, 239 (1991).
- In the odd-parity category, with the chiral symmetry is possible, as classified in Ref. References. However, this state seems unlikely because the quadratic line node giving at zero field contradicts the observed at low temperatures, in addition to the Pauli-paramagnetic effect in any direction mentioned above.
- K. Machida, M. Ozaki, and T. Ohmi, J. Phys. Soc. Jpn. 65, 3720 (1996).