# Quadrupole-Driven Non-Fermi Liquid and Magnetic-Field Induced Heavy Fermion States in a Non-Kramers Doublet System

###### Abstract

Orbital degrees of freedom in condensed matters could play important roles in forming a variety of exotic electronic states by interacting with conduction electrons. In 4-electron systems, because of strong intra-atomic spin-orbit coupling, an orbitally degenerate state inherently carries quadrupolar degrees of freedom. The present work has focussed on a purely quadrupole-active system PrIrZn showing superconductivity in the presence of an antiferroquadrupole order at = 0.11 K. We observed non-Fermi liquid (NFL) behaviors emerging in the electrical resistivity and the 4 contribution to the specific heat, , in the paramagnetic state at . Moreover, in magnetic fields 6 T, all data set of and are well scaled with characteristic temperatures ’s. This is the first observation of the NFL state in the nonmagnetic quadrupole-active system, whose origin is intrinsically different from that observed in the vicinity of the conventional quantum critical point. It implies possible formation of a quadrupole Kondo lattice resulting from hybridization between the quadrupoles and the conduction electrons with an energy scale of . At 0.13 K, and exhibit anomalies as approaches 5 T. This is the manifestation of a field-induced crossover toward a Fermi-liquid ground state in the quadrupole Kondo lattice.

## I Introduction

In metals and alloys, localized and/or electrons inherently possessing ‘spin’ and ‘orbital’ degrees of freedom are dominant sources of not only various magnetic phenomena but also unconventional superconductivity. Extensive studies on a lot of ‘spin’-active systems have successfully revealed the variety of the magnetic phenomena arising from competition of the inter-site Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction and the on-site magnetic Kondo effect, which was represented as the well-known ‘Doniach diagram’. On the other hand, there is less variety of researches on ‘orbital’-driven physics. Although the ‘orbital’ itself must have potential variety, there are few evidences of interplay between the orbital degrees of freedom and the conduction electron.

In the -electron systems, strong intra-atomic spin-orbit coupling forces the spin and orbital degrees of freedom to be described in terms of the total angular momentum J. When the magnetic moment due to J interacts strongly with itinerant conduction electrons, physically observable quantities follow the Fermi-liquid (FL) model of Landau, which is known as heavy fermion state. When some conditions were fulfilled unexpectedly, the FL state could become unstable, instead, an anomalous metallic state would emerge, so-called non-Fermi liquid (NFL) state.Lohneysen (); Umeo (); Si01 (); Custers03 () In the vicinity of the quantum critical point, there emerges unconventional superconductivity.Mathur98 ()

On the other hand, emergence of a different type of NFL state was predicted theoretically, in case electric quadrupoles of the localized electrons, which become active only in an orbitally degenerate state, interact with the conduction electrons, that is the quadrupole (two-channel) Kondo effect.Cox98 () It is quite different from the ordinary (single-channel) Kondo effect in terms of scattering process of the conduction electrons; the scattering source is not magnetic dipole but the time-reversal electric quadrupole. Thereby, the impurity quadrupoles were overcompensated by the conduction electrons, leading to the residual entropy of (0.5)ln2 and the NFL behavior: the magnetic susceptibility and the specific heat divided by temperature, , show ln dependence, and the electrical resistivity follows , where is residual resistivity and is a coefficient.Affleck91 (); Cox98 () NFL state and the residual entropy were also expected in a quadrupole Kondo ‘lattice’, in which quadrupole moments are periodically placed.Jarrell96 () Very recently, the peculiar temperature dependence of and have been theoretically predicted.Tsuruta99 (); Tsuruta15 () Moreover, in the lattice model, new types of electronic ordered states have been proposed.Hoshino11 (); Hoshino13 (); Hastatic () Nevertheless, there is no experimental evidence on the exotic states, probably because there are rare systems carrying purely active quadrupoles.

Over the past few decades, experimental efforts have been paid to address the issues that could arise from the impurity quadrupole Kondo effect.Stewart01 () NFL behaviors were observed in uranium- and praseodymium-based systems with 5 and 4 configurations, respectively, such as UYPdSeaman91 (); McEwen03 (), UBeMcElfresh93 (), UThBeAliev93 (), UThRuSiAmitsuka94 (), and PrLaPbKawae06 (). However, there is no firm experimental evidence for the impurity quadrupole Kondo effect, since atomic disorder or uncertainty in the crystalline electric field (CEF) ground states particularly in the U-based systems prevent the clarification.

Recently, a family of the praseodymium-based systems Pr (: transition metals, : Al, Zn, and Cd) have emerged as a prototype to study the quadrupole Kondo effect, since the CEF ground state is a non-Kramers doubly degenerated state, that is labeled as the doublet in the cubic point group, having no magnetic dipoles but electric quadrupoles.Nasch97 (); Onimaru11 (); Sakai11 (); Sato12 (); Iwasa13 (); Yazici15 () In PrIrZn, an AFQ order occurs at = 0.11 K.Onimaru11 () Although the entropy release of ln2 is expected from an order of the doubly degenerated CEF ground state, the entropy at estimated from the specific heat is only 20% of ln2. Therefore, there should be another mechanism, except the AFQ order, which consumes the rest of the entropy above . On the other hand, below , a superconducting transition sets in at = 0.05 K, suggesting a possible interplay between the quadrupole fluctuations and the superconducting Cooper-pair formation.Onimaru10 (); Onimaru11 (); Ishii11 (); Higemoto12 () The coexistence of superconductivity with quadrupole order also manifests in isostructural compounds PrRhZn, PrTiAl, and PrVAl.Onimaru12 (); Sakai12 (); Tsujimoto14 (); Matsubayashi12 () In PrTiAl and PrVAl, the strong hybridization between the 4 elections and the conduction electrons was revealed by the Al-NMR and Pr 3-4 resonant photoemission measurements.Tokunaga13 (); Matsunami11 () The large Seebeck coefficient divided by temperature observed in PrAl ( = Ti, V and Ta) also suggests the strong hybridization effect, although that for PrIrZn at is two or three orders smaller than those of PrAl.Machida15 () PrVAl exhibits NFL behavior above = 0.6 K; the 4 contributions in both and follow for 220 KSakai11 (), however, in the temperature range, the magnetic degrees of freedom in the first excited state must interfere in the NFL behavior. Quadrupolar quantum criticality induced by application of magnetic field has been proposed.Shimura15 ()

In the present work, we study the transport, thermodynamic and magnetic properties for the non-Kramers doublet system PrIrZn in magnetic fields applied along the [100] direction. For B [100], the AFQ order collapses at around 5 T, whose value is much lower than the critical fields for the AFQ order in the [110] and [111] field directions.Onimaru11 (); Ishii11 (); Hattori14 () Since high quality of the sample is necessary to reveal the inherent behavior arising from the quadrupolar degrees of freedom as described above, we have chosen single-cryslltaine samples with low residual resistivity of 0.2 cm (residual resistivity ratio, RRR100).Onimaru11 () The manifestation of the AFQ order guarantees the high-quality of the crystals, otherwise it is destroyed by the atomic disorder. Our measurements of the specific heat and the electrical resistivity in magnetic fields have revealed emergence of the NFL state possibly arising from formation of a quadrupole Kondo lattice. Furthermore, we have observed a magnetic-field induced crossover from the NFL state to an exotic heavy-fermion ground state in the quadrupole Kondo lattice.

## Ii Experimental

Single-crystalline samples of PrIrZn used in the present work were grown by using the melt-growth method described in the previous papers.Onimaru10 () Magnetization was measured using a commercial SQUID magnetometer (Quantum Design MPMS) between 1.9 and 350 K in magnetic fields up to 5 T. Magnetization measurements at low temperatures down to 0.045 K were performed by a capacitive Faraday method using a high-resolution capacitive force-sensing device installed in a He-He dilution refrigerator.Sakakibara94 () Electrical resistance was measured by a standard four-probe dc method in a laboratory-built system with a He-He dilution refrigerator. Thermopower was measured using a laboratory-made probe by applying a temperature difference of 0.04-0.3 K along a bar-shaped sample. Specific heat in magnetic fields was measured by a relaxation method at temperatures between 0.4 K and 300 K. The measurements at low temperatures down to 0.06 K were done under quasi-adiabatic conditions with a He-He dilution refrigerator equipped with a superconducting magnet of 12 T.

## Iii Results and discussion

### iii.1 Specific heat

Figure 1 (a) shows the temperature dependence of the 4 contribution to the specific heat, , of PrIrZn in magnetic fields 12 T applied along the [100] direction. To estimate the 4 contribution to the total specific heat , we subtracted the nuclear and phonon contributions and as described below. The hamiltonian of a nuclear spin of a Pr nucleus, 5/2, in magnetic field B can be represented as

(1) |

where 0.052 K is a coupling constant of hyper-fine interaction of a Pr ion Kondo61 (), and J, 1.72, and are the total angular momentum of 4 electrons, nuclear -factor, and nuclear magneton, respectively. The component of J along the magnetic field direction was regarded as the isothermal magnetization measured at 0.045 K. Thereby, the nuclear specific heat was estimated by taking the eigenvalues of the hamiltonian into consideration. The nuclear contribution is zero in 0. Applying magnetic field, it gradually increases on cooling, , in = 3 T, is increased up to about 3.5 J/K mol at 0.1 K. The contribution of the phonon was subtracted by using the specific heat of the La analog LaIrZn. The main panel of Fig. 1 (b) shows the entropy estimated by integrating . The value of reaches ln2 at 2 K, supporting that the physical properties at 2 K are governed by the doublet ground state.

The sharp peak at 0.11 K for = 0 in Fig. 1 (a) is the manifestation of the AFQ order.Onimaru11 (); Ishii11 (); Higemoto12 () As shown with the open triangles in the inset of Fig. 1 (b), upon applying magnetic field for 3 T, the peak shifts to lower temperatures and disappears at 4.5 T. The behavior of the specific heat peak at in the magnetic fields is much different from the previous report where the peak is split into two for 13 T.Onimaru11 () This inconsistency results from an extrinsic effect of another grain contained in the previous sample whose crystalline axis is directed away from the -axis of the dominant crystal.

On the other hand, as shown with the open triangles in Fig. 1 (a), in = 6 T shows a broad peak at 0.6 K, which shifts to 2 K with increasing up to 12 T. The height and width can be explained by taking account of the Zeeman splitting of the ground state doublet. Because the split singlets lose quadrupolar degrees of freedom, no phase transition occurs in 6 T.

We pay attention to another broad peak in ( = 0) at around 0.4 K in Fig. 1(a). When is plotted vs ln in the inset of Fig. 2(a), we find the ln dependence between 0.2 and 0.8 K. Since this ln dependence of emerges below 2 K, it arises from the degrees of freedom of the ground state. Here, we define a characteristic temperature, , as the temperature where reaches ln2. This definition follows the Cox’s definition of the Kondo temperature for the impurity quadrupole Kondo system.Cox98 () As shown in Fig. 1(b), a defined remains at around 0.4 K up to 3 T and increases for 3 T (see Fig. 5 (a)). We plot in Fig. 2 (a) the values of at various fields, which are well scaled with respect to in the range of 0.8. This scaling behavior will be discussed in detail by combining the NFL behavior of the electrical resistivity.

### iii.2 Electrical resistivity

The temperature dependence of the electrical resistivity in magnetic fields between 0 and 9 T applied along the [100] direction are shown in Fig. 3. It is noted that all the data are plotted without offset. In 4 T, shows upward convex curvature at 1.0 K, that is the NFL behavior. On cooling below 0.5 K, is likely to follow the variation as shown with the dashed lines. The absolute value of for 4 T is increased with increasing the magnetic field. On the other hand, above 6 T, shows downward curvature. The residual resistivity is increased with increasing the magnetic field.

We pay our attention again to the NFL behavior of the upward convex curve at 1.0 K. Here, the data of are replotted in the inset of Fig. 2 (b) as vs , where is the residual resistivity. In the region where stays at 1.0, follows . The arrows denote the characteristic temperature, , where starts deviating from the dependence on heating. We represent the data of vs in Fig. 2 (b), where follows the scaling well in the temperature range of 0.5 3. As shown in the inset of Fig. 2 (b), stays at around 0.35 K up to 4 T, and increases significantly once exceeds 4 T (see the (red) opened diamonds in Fig. 5 (a)). The field dependence of coincides with that of described in the previous subsection, suggesting the same origin of the both and .Stewart01 ()

### iii.3 Non-Fermi liquid behavior for

Let us discuss possible mechanisms for the observed NFL behaviors of and in the magnetic fields well scaled with respect to and , respectively. One is the contribution of a rattling phonon, which gives rise to .Dahm07 () There is an optical phonon excitation at around 7 meV (80 K), which is attributed to the low-energy vibration of the Zn atom.Hasegawa12 (); Wakiya15 (); Wakiya16 () However, as this phonon excitation probably ceases at 1 K, interaction of the optical phonon mode with the conduction electrons must be too weak to modify the low-temperature electronic transport.

The second is the impurity quadrupole Kondo effect predicted by Cox et al., which leads to and as mentioned in the introduction.Cox98 () The dependence of may be applicable to the present data. However, in the present case of PrIrZn, the quadrupole moments are not included as impurities but placed periodically, therefore, formation of the quadrupole Kondo lattice is a promising candidate bringing about the NFL state.Jarrell96 () The theoretical analyses on the temperature variations of and with the two-channel Anderson lattice model have led the following relations;

(2) | |||||

(3) |

where and (1 and 2) are parameters.Tsuruta15 () The calculations are shown by the (red) solid curves in Fig. 2 (a) and (b). The curves are in better agreement with our experimental data for wider temperature range than those calculated by the impurity quadrupole Kondo model as shown with the (blue) dotted line.

In the present case, the NFL behavior does not persist down to the low-temperature limit. This is certain because this NFL behavior does not arise from the conventional quantum critical point but from the quadrupole Kondo effect as described in the introduction. The quadrupoles are over-screened by the conduction electrons, leading to the unstable electronic state with the residual entropy of (0.5)ln2 at zero temperature in principle. In real systems, however, the residual entropy should go to zero at = 0 by a mechanism following the third law of thermodynamics. In fact, in the magnetic fields 4 T, the residual entropy is released by the intersite quadrupole interaction, leading to occurrence of the quadrupole order. The NFL behavior certainly manifests itself in the limited temperature range, where the magnetic entropy goes down from ln2 to (0.5)ln2. In zero magnetic field, where the NFL behaviors of the electrical resistivity and the specific heat appear in the temperature range from 1.5 K to 0.2 K.

Taking these experimental results and the analyses on them, this is the first observation of the NFL state in the nonmagnetic quadrupole-active system, whose origin is intrinsically different from that observed in the vicinity of the conventional quantum critical point. The scaling plots of and using the characteristic temperatures and , respectively, as shown in Fig. 2(a) and (b), imply the possible formation of the quadrupole Kondo lattice resulting from hybridization between the quadrupoles and the conduction electrons with an energy scale of .

### iii.4 Magnetic-field induced Fermi liquid state

Looking close to the data of at = 4 T as shown with the solid triangles in Fig. 4, we find another knee at = 0.12 K above = 0.08 K. In contrast to the drop of at 3 T, stays at around 0.13 K for 4 and 5 T. As shown in the inset of Fig. 1 (b), this coincides with the broad peak of at = 0.12 K in = 4.5 and 5 T. The field dependences of observed in the and measurements are plotted with the blue and red closed circles, respectively, in the phase diagram in Fig. 5 (a). The line can be a boundary of two states, meaning a phase transition or a crossover. In the present case, the crossover is likely to occur at , because the peak of becomes very broad despite the rather robust against the magnetic field. On cooling below , follows at low temperature for 46 T, as shown in the inset of Fig. 4. The coefficient as a function of is plotted in Fig. 5 (c). It is strongly enhanced at around 4.5 T, where at 0.07 K is also peaked. Moreover, there is a peak at 5.5 T in the Seebeck coefficient divided by temperature, , at 0.08 K.Ikeura14 () The coincidence of these peaks at around 5 T suggests a peculiar feature of the heavy fermion state on cooling through in a narrow range of , 3.56 T.

Fig. 6 (a) and (b) show the isothermal magnetization and the derivative of the magnetization with respect to the magnetic field, dd, at 0.045, 0.08, 0.17, 0.55, and 1.2 K in the magnetic field of 14 T applied along the [100] direction, respectively. The isothermal magnetization shows metamagnetic behavior at around 5 T, where dd shows a peak. The peak shifts slightly to higher fields and becomes broader with increasing temperatures, which is likely to be connected to the Schottky anomaly due to the splitting of the doublet in the magnetic field for 6 T.

The contour plot in Fig. 5 (a) indicates the derivative of the magnetization with respect to the magnetic field, dd. Note that the dd is largely enhanced only in the narrow range of , 3.55.5 T. Furthermore, the elastic modulus of the mode, which is a part of the -type mode, shows strong softening in the magnetic field at around 5 T.Ishii11 () All these observations suggest development of an exotic Fermi liquid state on cooling through in a range of , 3.55.5 T.

There are some possible mechanisms to form the exotic Fermi liquid state at . One is a crossover from the NFL state to a FL state expected in the two-channel Kondo model taking in magnetic field effect.Yotsuhashi02 () Perhaps, emergence of a FL state accompanied with a free magnetic spin in the vicinity of the NFL state has been pointed out by taking into account possible competition between the magnetic and quadrupole Kondo effects on the basis of an extended two-channel Kondo model.Kusunose14 () Another is a cross-over due to the composite electronic order Hoshino11 (); Hoshino13 () or the Hastatic order Hastatic () in the quadrupole Kondo lattice as were mentioned above. In any cases, it is noted that the hybridization between the quadrupoles and the conduction electrons gives rise to a new type of the electronic ground state. It is highly desirable to detect the order parameter by microscopic techniques.

Indeed, we expected that the quadrupole induced quantum critical behavior in would appear when the AFQ order is fully suppressed. However, as shown in Fig. 5 (b), does not diverse in the vicinity of 4.5 T on cooling to 0.07 K at all. It seems that the Doniach picture is not valid for this pure quadrupole system. Although it is still an open question why the quantum criticality was not detected in PrIrZn, a peculiar feature of the quadrupole must give rise to the unique electronic phenomena.

## Iv Summary

We report the low-temperature transport, thermodynamic and magnetic properties on a cubic system PrIrZn in which the non-Kramers doublet ground state has purely the quadrupolar degrees of freedom. In the moderately wide temperature range at , non-Fermi liquid behaviors were clearly observed in the electrical resistivity and the specific heat . Both and in magnetic fields 6 T applied along the [100] direction can be well scaled with a characteristic temperature , suggesting the formation of the quadrupole Kondo lattice due to the hybridization between the quadrupoles and the conduction electrons. It indicates that the NFL behavior observed in the present system is intrinsically different from that observed in the vicinity of the conventional quantum critical point. Furthermore, and exhibit anomalies at 0.13 K in the vicinity of 5 T, where the coefficient for , , and , have significant enhancement as a function of . The concomitant increase in dd indicates formation of a magnetic-field induced Fermi liquid ground state to remove the residual entropy in the quadrupole Kondo lattice. These observations imply that the Doniach picture relevant to the spin Kondo systems should not be valid for this purely quadrupole-active system, possibly leading us beyond the Doniach picture.

## Acknowledgements

The authors would like to thank A. Tsuruta, K. Miyake, H. Kusunose, S. Hoshino, J. Otsuki, Y. Kuramoto, K. Uenishi, Y. Yamane, I. Ishii, T. Suzuki and K. Iwasa for helpful discussions. We also thank Y. Shibata for the electron-probe microanalysis performed at N-BARD, Hiroshima University. The magnetization measurements with MPMS and specific heat measurements with PPMS were carried out at N-BARD, Hiroshima University. This work was financially supported by JSPS KAKENHI Grant Numbers 21102516, 23102718, 26707017, and 15H05884, 15H05886 (J-Physics), and by The Mazda Foundation Research Grant, Japan and Hiroshima University Fujii Research Fund.

## References

- (1) H. v. Löhneysen, A. Rosch, M. Vojta, and P. Wolfle, Rev. Mod. Phys. 79, 1015 (2007).
- (2) K. Umeo, H. Kadomatsu, and T. Takabatake, J. Phys.: Cond. Matter 8, 9743 (1996).
- (3) J. Custers, P. Gegenwart, H. Wilhelm, K. Neumaier, Y. Tokiwa, O. Trovarelli, C. Geibel, F. Steglich, C. Pépin, and P. Coleman, Nature 424, 524 (2003).
- (4) Q. Si, S. Rabello, K. Ingersent, J. L. Smith, Nature 413, 804 (2001).
- (5) N. D. Mathur, F. M. Grosche, S. R. Julian, I. R. Walker, D. M. Freye, R. K. W. Haselwimmer, and G. G. Lonzarich, Nature 394, 39 (1998).
- (6) D. L. Cox, Z. Zawadowski, Adv. Phys. 47, 599 (1998).
- (7) I. Affleck and A. W. W. Ludwig, Nucl. Phys. B 352, 849 (1991).
- (8) M. Jarrell, H. Pang, D. L. Cox, and K. H. Luk, Phys. Rev. Lett. 77, 1612 (1996).
- (9) A. Tsuruta, A. Kobayashi, Y. Ōno, T. Matsuura, and Y. Kuroda, J. Phys. Soc, Jpn. 68, 2491 (1999).
- (10) A. Tsuruta and K. Miyake, J. Phys. Soc. Jpn. 84, 114714 (2015)
- (11) S. Hoshino, J. Otsuki, and Y. Kuramoto, Phys. Rev. Lett. 107, 247202 (2011).
- (12) S. Hoshino, J. Otsuki, and Y. Kuramoto, J. Phys. Soc. Jpn. 82, 044707 (2013).
- (13) P. Chandra, P. Coleman, and R. Flint, Nature 493, 621 (2013)
- (14) G. R. Stewart, Rev. Mod. Phys. 73, 797 (2001).
- (15) C. L. Seaman, M. B. Maple, B. W. Lee, S. Ghamaty, M. S. Torikachvili, J. -S. Kang, L. Z. Liu, J. W. Allen, and D. L. Cox, Phys. Rev. Lett. 67, 2882 (1991).
- (16) K. A. McEwen, J. -G. Park, A. J. Gipson, and G. A. Gehring, J. Phys.: Condens. Matter 15, S1923 (2003).
- (17) M. McElfresh, M. B. Maple, J. O. Willis, D. Schiferl, J. L. Smith, Z. Fisk, and D. L. Cox, Phys. Rev. B 48, 10395 (1993).
- (18) FF. G. Aliev, S. Vieira, R. Villar, H. P. van der Meulen, K. Bakker, A. V. Andreev, JETP Lett. 58, 762 (1993).
- (19) H. Amitsuka, and T. Sakakibara, J. Phys. Soc. Jpn. 63, 736 (1994).
- (20) A. A. Yatskar, W. P. Beyermann, R. Movshovich, and P. C. Canfield, Phys. Rev. Lett. 77, 3637 (1996).
- (21) H. Tanida, H. S. Suzuki, S. Takagi, H. Onodera, and K. Tanigaki, J. Phys. Soc. Jpn. 75, 073705 (2006).
- (22) P. Morin, D. Schmitt, and du Tremolet de Lacheisserie, J. Magn. Magn. Mater. 30, 257 (1982).
- (23) T. Tayama, T. Sakakibara, K. Kitami, M. Yokoyama, K. Tenya, H. Amitsuka, D. Aoki, Y. Ōnuki, and Z. Kletowski, J. Phys. Soc. Jpn. 70, 248 (2001).
- (24) T. Onimaru, T. Sakakibara, N. Aso, H. Yoshizawa, H. S. Suzuki, and T. Takeuchi, Phys. Rev. Lett. 94, 197201 (2005).
- (25) T. Kawae, K. Kinoshita, Y. Nakaie, N. Tateiwa, K. Takeda, H. S. Suzuki, and T. Kitai, Phys. Rev. Lett. 96, 027210 (2006).
- (26) T. Nasch, W. Jeitschko, and U. C. Rodewald, Z. Naturforsch. B 52, 1023 (1997).
- (27) T. Onimaru, K. T. Matsumoto, Y. F. Inoue, K. Umeo, T. Sakakibara, Y. Karaki, M. Kubota, and T. Takabatake, Phys. Rev. Lett. 106, 177001 (2011).
- (28) A. Sakai and S. Nakatsuji, J. Phys. Soc. Jpn. 80, 063701 (2011).
- (29) T. J. Sato, S. Ibuka, Y. Nambu, T. Yamazaki, T. Hong, A. Sakai, and S. Nakatsuji, Phys. Rev. B 86, 184419 (2012)
- (30) K. Iwasa, H. Kobayashi, T. Onimaru, K. T. Matsumoto, N. Nagasawa, T. Takabatake, S. Ohira-Kawamura, T. Kikuchi, Y. Inamura, and K. Nakajima, J. Phys. Soc. Jpn. 82, 043707 (2013).
- (31) D. Yazici, T. Yanagisawa, B. D. White, and M. B. Maple, Phys. Rev. B 91, 115136 (2015).
- (32) I. Ishii, H. Muneshige, Y. Suetomi, T. K. Fujita, T. Onimaru, K. T. Matsumoto, T. Takabatake, K. Araki, M. Akatsu, Y. Nemoto, T. Goto, and T. Suzuki, J. Phys. Soc. Jpn. 80, 093601 (2011).
- (33) W. Higemoto, T. U. Ito, K. Ninomiya, T. Onimaru, K. T. Matsumoto, and T. Takabatake, Phys. Rev. B 85, 235152 (2012).
- (34) T. Onimaru, K. T. Matsumoto, Y. F. Inoue, K. Umeo, Y. Saiga, Y. Matsushita, R. Tamura, K. Nishimoto, I. Ishii, T. Suzuki, and T. Takabatake, J. Phys. Soc. Jpn. 79, 033704 (2010).
- (35) T. Onimaru, N. Nagasawa, K. T. Matsumoto, K. Wakiya, K. Umeo, S. Kittaka, T. Sakakibara, Y. Matsushita, and T. Takabatake, Phys. Rev. B 86, 184426 (2012).
- (36) A. Sakai, K. Kuga and S. Nakatsuji, J. Phys. Soc. Jpn. 81, 083702 (2012).
- (37) M. Tsujimoto, Y. Matsumoto, T. Tomita, A. Sakai, and S. Nakatsuji, arXiv:1407.0866.
- (38) K. Matsubayashi, T. Tanaka, A. Sakai, S. Nakatsuji, Y. Kubo, and Y. Uwatoko, Phys. Rev. Lett. 109 187004 (2012).
- (39) Y. Tokunaga, H. Sakai, S. Kambe, A. Sakai, S. Nakatsuji, and H. Harima, Phys. Rev. B 88, 085124 (2013).
- (40) M. Matsunami, M. Taguchi, A. Chainani, R. Eguchi, M. Oura, A. Sakai, S. Nakatsuji, and S. Shin, Phys. Rev. B 84, 193101 (2011).
- (41) Y. Machida, T. Yoshida, T. Ikeura, K. Izawa, A. Nakama, R. Higashinaka, Y. Aoki, H. Sato, A. Sakai, S. Nakatsuji, N. Nagasawa, K. Matsumoto, T. Onimaru, and T. Takabatake: J. Phys.: Conf. Ser. 592, 012025 (2015).
- (42) Y. Shimura, M. Tsujimoto, B. Zeng, L. Balicas, A. Sakai, and S. Nakatsuji, Phys. Rev. B 91, 241102(R) (2015).
- (43) K. Hattori and H. Tsunetsugu, J. Phys. Soc. Jpn. 83, 034709 (2014).
- (44) T. Sakakibara, H. Mitamura, T. Tayama, and H. Amitsuka, Jpn. J. Appl. Phys. 33, 5067 (1994).
- (45) Kondo, J. Internal magnetic field in rare earth metals. J. Phys. Soc. Jpn. 16, 1690-1691 (1961).
- (46) T. Dahm and K. Ueda, Phys. Rev. Lett. 99, 187003 (2007).
- (47) T. Hasegawa, N. Ogita, and M. Udagawa, J. Phys.: Conf. Ser. 391 012016 (2012).
- (48) K. Wakiya, T. Onimaru, S. Tsutsui, K. T. Matsumoto, N. Nagasawa, A. Q. R. Baron, T. Hasegawa, N. Ogita, M. Udagawa, and T. Takabatake, J. Phys.: Conf. Ser. 592, 012024 (2015).
- (49) K. Wakiya, T. Onimaru, S. Tsutsui, K. T. Matsumoto, N. Nagasawa, A. Q. R. Baron, T. Hasegawa, N. Ogita, M. Udagawa, and T. Takabatake, Phys. Rev. B 93, 064105 (2016).
- (50) T. Ikeura, T. Matsubara, Y. Machida, K. Izawa, N. Nagasawa, K. T. Matsumoto, T. Onimaru, and T. Takabatake, JPS Conf. Proc. 3, 011091 (2014).
- (51) S. Yotsuhashi and H. Maebashi, J. Phys. Soc. Jpn. 71, 1705 (2002).
- (52) H. Kusunose and T. Onimaru, J. Phys.: Conf. Ser. 592, 012099 (2015).