J. Phys. Soc. Jpn. 91, 043702 (2022) [5 Pages]
LETTERS

Two-Dimensional XY-Type Magnetic Properties of Locally Noncentrosymmetric Superconductor CeRh2As2

+ Affiliations
1Department of Physics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan2Max Planck Institute for Chemical Physics of Solids, D-01187 Dresden, Germany

We performed 75As-NMR measurements to investigate the normal-state magnetic properties of CeRh2As2, a recently-discovered heavy-fermion superconductor. The magnitude and temperature dependence of the Knight shift at the As(2) site indicate easy-plane-type magnetic anisotropy in CeRh2As2. With regard to spin fluctuations, the temperature dependence of the nuclear spin–lattice relaxation rate 1/T1 arising from the 4f electrons decreases from high-temperature constant behavior on cooling at ∼40 K, which is a typical behavior of heavy-fermion systems. In addition, 1/T1 becomes constant at low temperatures, suggesting spatially two-dimensional antiferromagnetic fluctuations. Two-dimensional magnetic correlations in the real space are quite rare among heavy-fermion superconductors, and they may be a key factor in the unique superconducting multi phase in CeRh2As2.

©2022 The Author(s)
This article is published by the Physical Society of Japan under the terms of the Creative Commons Attribution 4.0 License. Any further distribution of this work must maintain attribution to the author(s) and the title of the article, journal citation, and DOI.

Unconventional superconductivity is mostly considered to be mediated by magnetic fluctuations.13) Thus, investigating magnetic properties in the normal state provides important information about the superconducting (SC) mechanism.4,5)

CeRh2As2 is a recently discovered heavy-fermion superconductor whose \(T_{\text{SC}}\) is approximately 0.3 K.6) The heavy-fermion superconductivity is characterized by a broad maximum in resistivity at \(T_{\text{coh}}\sim 40\) K and a large specific-heat jump at \(T_{\text{SC}}\). In many Ce-based heavy-fermion compounds, superconductivity appears in a narrow region near the antiferromagnetic (AFM) quantum critical point as a function of tuning parameters such as pressure, and it is considered to be mediated by AFM fluctuations.2) Even in CeRh2As2, AFM fluctuations are expected to play an important role. Further, CeRh2As2 exhibits additional features that are not seen in other systems.

One such feature is the lack of local inversion symmetry in Ce layers with global inversion symmetry. The crystal structure is of the tetragonal CaBe2Ge2-type with space group \(P4/nmm\) (No. 129, \(D_{4h}^{7}\)).7) CeRh2As2 has two crystallographically inequivalent As and Rh sites; As(1) [Rh(1)] is tetrahedrally coordinated by Rh(2) [As(2)], as shown in Fig. 1(a). The crystal structure looks similar to that of the typical heavy-fermion superconductor CeCu2Si2 (i.e., a ThCr2Si2-type structure).8) Both compounds have inversion symmetry in crystal structures. However, the stacking order of the block layers is different. In CeRh2As2, a Ce layer is located between two different block layers (Rh–As–Rh and As–Rh–As), and there is no inversion center at the Ce site. In contrast, in CeCu2Si2, a Ce layer is located between identical block layers (Si–Cu–Si). Further, in CeMSi3 (M = Rh or Ir), a Ce layer is sandwiched between Si–Si–M block layers and global inversion symmetry is broken,9) making CeRh2As2 a good comparison to investigate how a lack of local inversion symmetry affects unconventional SC and magnetic states. For superconductors with no inversion center at particular atomic sites with global inversion symmetry, theoretical studies have noted that odd-parity SC states such as the pair density wave state can be stabilized under magnetic fields even if a pairing interaction of the spin-singlet channel only exists.10) In fact, a first-order SC–SC phase transition at approximately 4 T was reported for \(H\parallel c\);6) this can be interpreted as a phase transition from a low-field even-parity state to a high-field odd-parity state inside the SC phase.1113)


Figure 1. (Color online) (a) Crystal structure of CeRh2As2. NMR spectrum arising from the central transition (\(-1/2\leftrightarrow 1/2\)) of CeRh2As2 at 4.2 K for (b) \(H\parallel c\) and (c) \(H\perp c\). The dashed lines indicate the position of \(K = 0\) when the contribution of the nuclear quadrupole interaction is taken into account.

Another feature of CeRh2As2 is the nonmagnetic and magnetic phase transition just above and below \(T_{\text{SC}}\), respectively.6,14) The specific heat shows a large anomaly at \(T_{\text{SC}}\) and a rather weak one at \(T_{0}\sim 0.4\) K. Magnetization and AC susceptibility measurements show no anomaly at \(T_{0}\), and the \(T_{0}\) anomaly increases with magnetic fields applied perpendicular to the c axis, suggesting that it is not a simple magnetic order. On the basis of recent renormalized band structure calculations, a quadrupole-density-wave transition at \(T_{0}\) has been proposed.15) In addition to the nonmagnetic transition at \(T_{0}\), our nuclear quadrupole resonance (NQR) measurements revealed an AFM order inside the SC phase.14) The NQR spectrum of two crystallographically inequivalent As sites showed site-dependent linewidth broadening below \(T_{\text{N}} = 0.25\) K. Because the internal magnetic field can be canceled out at the symmetric site due to the superposition of the transferred fields from different Ce neighbors, the site-dependent linewidth broadening is the evidence of AFM transition. Furthermore, the orbital order due to the existence of two nonequivalent positions of Ce atoms was theoretically predicted.16) Therefore, CeRh2As2 is a promising system to study how the absence of local inversion symmetry induces or influences unconventional nonmagnetic, AFM and SC states, as well as their interaction.

In this study, we performed 75As-NMR measurements to investigate the normal-state magnetic properties of CeRh2As2. The temperature dependence of the Knight shift and \(1/T_{1}\) indicates spatially two-dimensional magnetic fluctuations with easy-plane (XY)-type anisotropy in CeRh2As2. Two-dimensional magnetic correlations in the real space are quite different from those of other heavy-fermion superconductors, as a result of which it exhibits a unique SC state.

Single crystals of CeRh2As2 and a non-magnetic reference compound LaRh2As2 were grown by the Bi flux method.6) The bulk magnetic susceptibility was measured using a commercial magnetometer with a superconducting quantum interference device (Quantum Design, MPMS). For NMR measurements, we used a split SC magnet, that generated a horizontal field and combined it with a single-axis rotator to apply a magnetic field parallel or perpendicular to the c axis. The 75As-NMR spectra (nuclear spin \(I=3/2\), nuclear gyromagnetic ratio \(\gamma/2\pi=7.29\) MHz/T, and natural abundance 100%) were obtained as a function of the magnetic field at a fixed frequency (∼26 MHz). Figures 1(b) and 1(c) show typical NMR spectra for \(H\parallel c\) and \(H\perp c\), respectively. To ascribe the NMR peaks to the As sites and estimate the Knight shift \(K_{i}\) (\(i = c\) and ⊥), we computed the resonance magnetic field by diagonalizing the following Hamiltonian. \begin{align} \mathcal{H} &= - \gamma \hbar (1 + K_{i})I\cdot H\notag\\ &\quad + \frac{h\nu_{\text{Q}}}{6}\left[3I_{z}^{2}-I(I+1)+\frac{\eta}{2}(I_{+}^{2}+I_{-}^{2})\right], \end{align} (1) where h, \(\nu_{\text{Q}}=\frac{3heQV_{{zz}}}{2I(2I-1)}\), and \(\eta = |\frac{V_{yy}-V_{xx}}{V_{zz}}|\) are the Planck constant, NQR frequency, and asymmetry parameter, respectively. η is zero at each As site because of the 4-fold symmetry of the atomic position. Because \(\nu_{\text{Q}}\) of the As(1) site (∼30 MHz) is larger than that of the As(2) site (∼10 MHz),14) the site assignments can be made from the angular dependence of the NMR peaks. The temperature dependence of \(\nu_{\text{Q}}\) was determined from NQR measurements and the interval between NMR satellite signals. For \(H\parallel c\), because the contribution of the nuclear quadrupole interaction can be ignored at the center peak, the Knight shift is simply determined by \(K = (H_{0}-H)/H\). Here, \(H_{0}\) is the reference field. It was estimated using 63Cu (\(^{63}\gamma_{n}/2\pi = 11.285\) MHz/T) and 65Cu (\(^{65}\gamma_{n}/2\pi = 12.08\) MHz/T) NMR signals from the NMR coil. The nuclear spin–lattice relaxation rate \(1/T_{1}\) was determined by fitting the time variation of the nuclear magnetization probed with the spin–echo intensity after saturation to a theoretical function for \(I = 3/2\).17,18)

Figures 2(a) and 2(b) show the temperature dependence of the Knight shift in CeRh2As2 and LaRh2As2. \(K_{c}\) and \(K_{\bot}\) denote the Knight shift for \(H\parallel c\) and \(H\perp c\), respectively. Because \(\nu_{\text{Q}}\) at the As(1) site is large (∼30 MHz), \(K_{\bot}\) at the As(1) site could not be measured in the present measurement setup. For comparison, the temperature dependence of the magnetic susceptibility \(\chi_{i}\) (\(i = c\) and ⊥) is also shown. The Knight shift is proportional to the magnetic susceptibility as given by, \begin{equation} K_{i} = A_{\text{hf},i} \chi_{\text{spin},i} + K_{\text{orb},i}, \end{equation} (2) where \(A_{\text{hf},i}\), \(\chi_{\text{spin},i}\), and \(K_{\text{orb},i}\) are respectively the hyperfine coupling constant, spin susceptibility, and orbital part of the Knight shift in each direction (which is generally temperature independent). \(K_{i}\) in CeRh2As2 for both magnetic field directions increased on cooling. While \(K_{c}\) shows a broad maximum at 3 K [see the inset of Fig. 2(a)], \(K_{\bot}\) continued to increase down to 1.5 K. This behavior is consistent with the bulk magnetic susceptibility, which shows the XY-type magnetic anisotropy. Most heavy-fermion superconductors such as CeCu2Si219,20) and CeIrSi321,22) have this kind of XY-type magnetic anisotropy. Unlike in the case of CeRh2As2, \(K_{i}\) in LaRh2As2 slightly decreased on cooling. In conventional metals, \(K_{i}\) remains constant with changes in temperature. The weak temperature-dependent \(K_{i}\) in LaRh2As2 can be explained by the increase in the density of states due to thermally excited quasi-particles. This suggests the presence of a peak of the density of states near the Fermi energy \(E_{\text{F}}\), as in the case of iron-based superconductors.23,24) We note that \(K_{c}\) at the As(1) site in LaRh2As2 is larger than that in CeRh2As2. The reason for this is unclear, and this issue requires further study.


Figure 2. (Color online) Temperature dependence of Knight shift \(K_{i}\) (\(i = c\) and ⊥) for (a) \(H\parallel c\) and (b) \(H\perp c\) in CeRh2As2. For comparison, the temperature dependence of \(K_{i}\) in LaRh2As2 and magnetic susceptibility \(\chi_{i}\) are also depicted. Measurements were performed at ∼4 T. The \(K = 0\) line is indicated by a solid line. (Inset) Enlarged view in low-temperature range.

The hyperfine coupling constant \(A_{\text{hf},i}\) can be determined by the slope of the Kχ plot, as shown in Figs. 3(a) and 3(b). This slope changes between regions where \(\chi_{i}\) is small and large. This behavior is often observed in Ce-based heavy-fermion systems, and it is attributed to the development of the heavy-fermion state owing to cf hybridization.19,25) Indeed, the position of the slope change is roughly coincident with \(T_{\text{coh}}\), as indicated by the solid lines in Figs. 3(a) and 3(b). Table I summarizes the estimated hyperfine coupling constants \(A_{\text{hf}}\) and \(K_{\text{orb}}\). \(A_{\text{hf},c}\) at the As(2) site becomes larger on cooling, whereas \(A_{\text{hf},c}\) and \(A_{\text{hf},\bot}\) at the As(1) site become smaller on cooling. Assuming the classical dipole interaction, \(A_{\text{hf},i}\) was calculated to be 0.02–0.06 T/\(\mu_{\text{B}}\), which is of one order of magnitude smaller than the experimentally estimated \(A_{\text{hf},i}\). Therefore, the transferred hyperfine interaction is dominant in CeRh2As2.


Figure 3. (Color online) Kχ plot for (a) \(H\parallel c\) and (b) \(H\perp c\) in CeRh2As2. The broken lines indicate the linear fitting results. The solid lines indicate \(\chi_{i}\) value at \(T_{\text{coh}}\sim 40\) K.

Data table
Table I. Estimated hyperfine coupling constants \(A_{\text{hf}}\) and \(K_{\text{orb}}\). Temperatures in parentheses indicate the fitted temperature range.

Next, we discuss the spin dynamics in CeRh2As2. Figures 4(a) and 4(b) show the temperature dependence of \(1/T_{1}\). \(1/T_{1}\) of CeRh2As2 showed a complex temperature dependence; specifically, it decreased slowly at high temperatures, decreased significantly below ∼100 K, and finally became constant below 10 K. In contrast, \(1/T_{1}\) of LaRh2As2 is roughly proportional to temperature in both magnetic-field directions; this is the typical behavior of conventional metals. To extract the contribution of Ce \(4f\) moments \((1/T_{1})_{4f}\), \(1/T_{1}\) of LaRh2As2 was subtracted from that of CeRh2As2 as shown in the inset of Fig. 4. Although the data above 100 K is scattered owing to an experimental error, the constant \((1/T_{1})_{4f}\) at high temperatures reflects the localized \(4f\) moments. On cooling, \((1/T_{1})_{4f}\) decreased below 40 K, which is consistent with the broad maximum of the electrical resistivity at \(T_{\text{coh}}\sim 40\) K.6) Thus, this decrease in \((1/T_{1})_{4f}\) is related to the development of the heavy-fermion state. On further cooling, \((1/T_{1})_{4f}\) becomes constant again, unlike in the case of other heavy-fermion systems. Because the Ce \(4f\) moments have an itinerant character at low temperatures, the temperature dependence of \(1/T_{1}\) can be understood by the self-consistent renormalization (SCR) theory.4,26,27) According to the SCR theory, a constant \(1/T_{1}\) implies spatially two-dimensional AFM fluctuations, as is often observed in high-\(T_{\text{SC}}\) cuprates28,29) and iron-based superconductors,30,31) but not in heavy-fermion superconductors. In CeCu2Si2 and CeIrSi3, \(1/T_{1}\) is proportional to T20) and \(\sqrt{T}\)32) just above \(T_{\text{SC}}\), indicating a weakly correlated metallic state and the development of spatially three-dimensional AFM fluctuations, respectively. In addition, the \(1/T_{1}\) at low temperatures in the NQR measurement is almost the same as that in the NMR measurement, as shown in Fig. 4(a), and thus, the two-dimensional AFM fluctuations are related to the AFM order inside the SC phase rather than the magnetic field-induced magnetic fluctuations. Because some theories suggest that an SC multi phase appears only when in-plane interactions dominate over c-axis interactions,6,10) the two-dimensional magnetic correlation is considered to be the key feature in CeRh2As2. Indeed, some band calculations also indicate the quasi-two-dimensional electronic properties in this system.12,15,33)


Figure 4. (Color online) Information about magnetic fluctuations. Temperature dependence of \(1/T_{1}\) for (a) \(H\parallel c\) and (b) \(H\perp c\) in CeRh2As2 and LaRh2As2. The temperature dependence of \(1/T_{1}\) in NQR measurements is also depicted. The insets show the contribution from Ce \(4f\) moments \((1/T_{1})_{4f}\). The arrows indicate \(T_{\text{coh}}\). (c) Temperature dependence of \(R_{a}\) and \(R_{c}\) at As(2) site in CeRh2As2.

The XY-type magnetic anisotropy was also confirmed from the \(T_{1}\) measurements. In general, \(1/T_{1}\) divided by temperature \(1/T_{1}T\) can be described in terms of fluctuating hyperfine fields perpendicular to the applied magnetic field parallel to the z axis as follow: \begin{equation} \left(\frac{1}{T_{1}T}\right)_{z} {}\propto R_{x} + R_{y}, \end{equation} (3) where \(R_{i}\equiv |H_{\text{hf},i}(\omega_{\text{res}})|^{2}\) and \(|X(\omega)|^{2}\) denotes the power spectral density of a time-dependent random variable \(X(t)\). Because both As sites have 4-fold symmetry of the atomic position, \(R_{a} = R_{b}\); thus, the relations \((1/T_{1}T)_{4f,c} = 2R_{a}\) and \((1/T_{1}T)_{4f,\bot} = R_{a} + R_{c}\) can be applied. Figure 4(c) shows the temperature dependence of \(R_{a}\) and \(R_{c}\) at the As(2) site in CeRh2As2 as estimated using the above relations. \(R_{c}\) was small and showed a weak temperature dependence with a broad maximum at ∼3 K, and \(R_{a}\) continued to increase on cooling, indicating XY-type magnetic anisotropy. This behavior is similar to static magnetic susceptibility. From the temperature dependence and the anisotropy of \(1/T_{1}T\), it is concluded that CeRh2As2 shows two-dimensional fluctuations with the XY-type magnetic anisotropy, thereby satisfying the condition of a unique SC multi phase.

The XY-type magnetic anisotropy is often observed in heavy-fermion superconductors. Sakai et al. investigated the anisotropy of spin fluctuations in many heavy-fermion systems and found that \(T_{\text{SC}}\) increases with increasing \(\Gamma_{c}/\Gamma_{a}\).34) Here, \(\Gamma_{i}\equiv(k_{\text{B}}\gamma_{n}^{2}A_{\text{hf},i}^{2}/2\pi R_{i})^{1/2}\), which is the local magnetic fluctuation energy in the i direction. It is a quantity that becomes smaller near the quantum critical point. \(k_{\text{B}}\) and \(\gamma_{n}\) are the Boltzmann constant and nuclear gyromagnetic ratio, respectively. In CeRh2As2, \(\Gamma_{c}/\Gamma_{a}\) was calculated to be ∼6 at 1.5 K, indicating a large XY-type anisotropy of spin fluctuations. According to the relation between \(\Gamma_{c}/\Gamma_{a}\) and \(T_{\text{SC}}\) in Ref. 34, this value corresponds to \(T_{\text{SC}} = 10\) K; this is much higher than the actual \(T_{\text{SC}}\) of ∼0.3 K. The small \(T_{\text{SC}}\) in CeRh2As2 might be related to the uniqueness of its crystal structure.

Finally, we comment on the relationship between the magnetic anisotropy and a magnetic structure in the AFM state inside the SC phase. From our NQR measurements, the A-type AFM (in-plane ferromagnetic and inter-plane antiferromagnetic) order with magnetic moments parallel to the c axis or a helical order with in-plane moments are promising candidates for the magnetic structure.14) Considering the XY-type magnetic anisotropy observed in the present study, a magnetic structure with an in-plane ordered moment would be preferable, but it is not so simple in Ce-based compounds. Magnetic anisotropy is mainly affected by single-site crystal field effects, but the direction of ordered moments is determined by the interaction between the moments. As a result, the direction of ordered moments and magnetic easy axis can be different in some Ce-based compounds.35,36) In addition, in CeRh2As2, another phase transition, which is considered to be the quadrupole-density-wave transition,15) exists above the AFM transition temperature, and the order parameter of this state might be coupled with the AFM moments. Therefore, detailed measurements of low-temperature properties is important for the determination of the AFM structure.

In conclusion, we performed 75As-NMR measurements to investigate the normal-state magnetic properties of CeRh2As2. The magnitude and temperature dependence of the Knight shift at the As(2) site indicated XY-type magnetic anisotropy in CeRh2As2. For spin fluctuations, spatially two-dimensional XY-type antiferromagnetic fluctuations were confirmed from the temperature dependence of \(1/T_{1}\). These results should be crucial for clarifying the origin of the SC multi phase and the unique coexistence of superconductivity and magnetism in CeRh2As2.

Acknowledgments

The authors would like to thank A. Ikeda for experimental supports, and K. Nogaki, Y. Yanase, S. Ogata, Y. Maeno, and S. Yonezawa for valuable discussions. This work was partially supported by the Kyoto University LTM Center, and Grants-in-Aid for Scientific Research (KAKENHI) (Grants Nos. JP15H05745, JP17K14339, JP19K14657, JP19H04696, and JP20H00130). C.G. and E.H. acknowledge support from the DFG through grant GE 602/4-1 Fermi-NEST.


References

  • 1 P. Monthoux, D. Pines, and G. G. Lonzarich, Nature 450, 1177 (2007). 10.1038/nature06480 CrossrefGoogle Scholar
  • 2 C. Pfleiderer, Rev. Mod. Phys. 81, 1551 (2009). 10.1103/RevModPhys.81.1551 CrossrefGoogle Scholar
  • 3 K. Ishida, Y. Nakai, and H. Hosono, J. Phys. Soc. Jpn. 78, 062001 (2009). 10.1143/JPSJ.78.062001 LinkGoogle Scholar
  • 4 Y. Nakai, T. Iye, S. Kitagawa, K. Ishida, S. Kasahara, T. Shibauchi, Y. Matsuda, H. Ikeda, and T. Terashima, Phys. Rev. B 87, 174507 (2013). 10.1103/PhysRevB.87.174507 CrossrefGoogle Scholar
  • 5 K. Ishida, S. Matsuzaki, M. Manago, T. Hattori, S. Kitagawa, M. Hirata, T. Sasaki, and D. Aoki, Phys. Rev. B 104, 144505 (2021). 10.1103/PhysRevB.104.144505 CrossrefGoogle Scholar
  • 6 S. Khim, J. F. Landaeta, J. Banda, N. Bannor, M. Brando, P. M. R. Brydon, D. Hafner, R. Küchler, R. Cardoso-Gil, U. Stockert, A. P. Mackenzie, D. F. Agterberg, C. Geibel, and E. Hassinger, Science 373, 1012 (2021). 10.1126/science.abe7518 CrossrefGoogle Scholar
  • 7 R. Madar, P. Chadouet, J. P. Senateur, S. Zemni, and D. Tranqui, J. Less-Common Met. 133, 303 (1987). 10.1016/0022-5088(87)90241-4 CrossrefGoogle Scholar
  • 8 F. Steglich, J. Aarts, C. D. Bredl, W. Lieke, D. Meschede, W. Franz, and H. Schäfer, Phys. Rev. Lett. 43, 1892 (1979). 10.1103/PhysRevLett.43.1892 CrossrefGoogle Scholar
  • 9 N. Kimura, K. Ito, K. Saitoh, Y. Umeda, H. Aoki, and T. Terashima, Phys. Rev. Lett. 95, 247004 (2005). 10.1103/PhysRevLett.95.247004 CrossrefGoogle Scholar
  • 10 T. Yoshida, M. Sigrist, and Y. Yanase, Phys. Rev. B 86, 134514 (2012). 10.1103/PhysRevB.86.134514 CrossrefGoogle Scholar
  • 11 E. G. Schertenleib, M. H. Fischer, and M. Sigrist, Phys. Rev. Res. 3, 023179 (2021). 10.1103/PhysRevResearch.3.023179 CrossrefGoogle Scholar
  • 12 K. Nogaki, A. Daido, J. Ishizuka, and Y. Yanase, Phys. Rev. Res. 3, L032071 (2021). 10.1103/PhysRevResearch.3.L032071 CrossrefGoogle Scholar
  • 13 D. Möckli and A. Ramires, Phys. Rev. Res. 3, 023204 (2021). 10.1103/PhysRevResearch.3.023204 CrossrefGoogle Scholar
  • 14 M. Kibune, S. Kitagawa, K. Kinjo, S. Ogata, M. Manago, T. Taniguchi, K. Ishida, M. Brando, E. Hassinger, H. Rosner, C. Geibel, and S. Khim, Phys. Rev. Lett. 128, 057002 (2022). 10.1103/PhysRevLett.128.057002 CrossrefGoogle Scholar
  • 15 D. Hafner, P. Khanenko, E.-O. Eljaouhari, R. Küchler, J. Banda, N. Bannor, T. Lühmann, J. F. Landaeta, S. Mishra, I. Sheikin, E. Hassinger, S. Khim, C. Geibel, G. Zwicknagl, and M. Brando, Phys. Rev. X 12, 011023 (2022). 10.1103/PhysRevX.12.011023 CrossrefGoogle Scholar
  • 16 A. Ptok, K. J. Kapcia, P. T. Jochym, J. Łażewski, A. M. Oleś, and P. Piekarz, Phys. Rev. B 104, L041109 (2021). 10.1103/PhysRevB.104.L041109 CrossrefGoogle Scholar
  • 17 A. Narath, Phys. Rev. 162, 320 (1967). 10.1103/PhysRev.162.320 CrossrefGoogle Scholar
  • 18 D. E. MacLaughlin, J. D. Williamson, and J. Butterworth, Phys. Rev. B 4, 60 (1971). 10.1103/PhysRevB.4.60 CrossrefGoogle Scholar
  • 19 T. Ohama, H. Yasuoka, D. Mandrus, Z. Fisk, and J. L. Smith, J. Phys. Soc. Jpn. 64, 2628 (1995). 10.1143/JPSJ.64.2628 LinkGoogle Scholar
  • 20 S. Kitagawa, T. Higuchi, M. Manago, T. Yamanaka, K. Ishida, H. S. Jeevan, and C. Geibel, Phys. Rev. B 96, 134506 (2017). 10.1103/PhysRevB.96.134506 CrossrefGoogle Scholar
  • 21 Y. Okuda, Y. Miyauchi, Y. Ida, Y. Takeda, C. Tonohiro, Y. Oduchi, T. Yamada, N. Duc Dung, T. D. Matsuda, Y. Haga, T. Takeuchi, M. Hagiwara, K. Kindo, H. Harima, K. Sugiyama, R. Settai, and Y. Ōnuki, J. Phys. Soc. Jpn. 76, 044708 (2007). 10.1143/JPSJ.76.044708 LinkGoogle Scholar
  • 22 H. Mukuda, T. Ohara, M. Yashima, Y. Kitaoka, R. Settai, Y. Ōnuki, K. M. Itoh, and E. E. Haller, Phys. Rev. Lett. 104, 017002 (2010). 10.1103/PhysRevLett.104.017002 CrossrefGoogle Scholar
  • 23 F. L. Ning, K. Ahilan, T. Imai, A. S. Sefat, R. Jin, M. A. McGuire, B. C. Sales, and D. Mandrus, J. Phys. Soc. Jpn. 78, 013711 (2009). 10.1143/JPSJ.78.013711 LinkGoogle Scholar
  • 24 S. Kitagawa, Y. Nakai, T. Iye, K. Ishida, Y. Kamihara, M. Hirano, and H. Hosono, Phys. Rev. B 81, 212502 (2010). 10.1103/PhysRevB.81.212502 CrossrefGoogle Scholar
  • 25 N. J. Curro, B. Simovic, P. C. Hammel, P. G. Pagliuso, J. L. Sarrao, J. D. Thompson, and G. B. Martins, Phys. Rev. B 64, 180514 (2001). 10.1103/PhysRevB.64.180514 CrossrefGoogle Scholar
  • 26 T. Moriya and K. Ueda, J. Phys. Soc. Jpn. 63, 1871 (1994). 10.1143/JPSJ.63.1871 LinkGoogle Scholar
  • 27 J. L. Sarrao, E. D. Bauer, J. N. Mitchell, P. H. Tobash, and J. D. Thompson, Physica C 514, 184 (2015). 10.1016/j.physc.2015.02.031 CrossrefGoogle Scholar
  • 28 S. Ohsugi, Y. Kitaoka, K. Ishida, and K. Asayama, J. Phys. Soc. Jpn. 60, 2351 (1991). 10.1143/JPSJ.60.2351 LinkGoogle Scholar
  • 29 K. Ishida, Y. Kitaoka, N. Ogata, T. Kamino, K. Asayama, J. R. Cooper, and N. Athanassopoulou, J. Phys. Soc. Jpn. 62, 2803 (1993). 10.1143/JPSJ.62.2803 LinkGoogle Scholar
  • 30 Y. Nakai, T. Iye, S. Kitagawa, K. Ishida, H. Ikeda, S. Kasahara, H. Shishido, T. Shibauchi, Y. Matsuda, and T. Terashima, Phys. Rev. Lett. 105, 107003 (2010). 10.1103/PhysRevLett.105.107003 CrossrefGoogle Scholar
  • 31 S. Kitagawa, T. Kawamura, K. Ishida, Y. Mizukami, S. Kasahara, T. Shibauchi, T. Terashima, and Y. Matsuda, Phys. Rev. B 100, 060503(R) (2019). 10.1103/PhysRevB.100.060503 CrossrefGoogle Scholar
  • 32 H. Mukuda, T. Fujii, T. Ohara, A. Harada, M. Yashima, Y. Kitaoka, Y. Okuda, R. Settai, and Y. Ōnuki, Phys. Rev. Lett. 100, 107003 (2008). 10.1103/PhysRevLett.100.107003 CrossrefGoogle Scholar
  • 33 D. C. Cavanagh, T. Shishidou, M. Weinert, P. Brydon, and D. Agterberg, Phys. Rev. B 105, L020505 (2022). 10.1103/PhysRevB.105.L020505 CrossrefGoogle Scholar
  • 34 H. Sakai, H. Chudo, Y. Tokunaga, S. Kambe, Y. Haga, F. Ronning, E. D. Bauer, J. D. Thompson, Y. Homma, D. Aoki, Y. Nakano, F. Honda, R. Settai, and Y. Ōnuki, J. Phys. Soc. Jpn. 81 [Suppl. B], SB003 (2012). 10.1143/JPSJS.81SB.SB003 LinkGoogle Scholar
  • 35 S. Araki, N. Metoki, A. Galatanu, E. Yamamoto, A. Thamizhavel, and Y. Ōnuki, Phys. Rev. B 68, 024408 (2003). 10.1103/PhysRevB.68.024408 CrossrefGoogle Scholar
  • 36 C. Krellner and C. Geibel, J. Cryst. Growth 310, 1875 (2008). 10.1016/j.jcrysgro.2007.10.045 CrossrefGoogle Scholar