# Break up of heavy fermions at an antiferromagnetic instability

###### Abstract

We present results of high-resolution, low-temperature measurements of the Hall coefficient, thermopower, and specific heat on stoichiometric YbRhSi. They support earlier conclusions of an electronic (Kondo-breakdown) quantum critical point concurring with a field induced antiferromagnetic one. We also discuss the detachment of the two instabilities under chemical pressure. Volume compression/expansion (via substituting Rh by Co/Ir) results in a stabilization/weakening of magnetic order. Moderate Ir substitution leads to a non-Fermi-liquid phase, in which the magnetic moments are neither ordered nor screened by the Kondo effect. The so-derived zero-temperature global phase diagram promises future studies to explore the nature of the Kondo breakdown quantum critical point without any interfering magnetism.

## 1 Different types of quantum critical points

Quantum phase transitions in matter arise due to competing interactions, which result in competing ground-state properties. When a quantum phase transition is continuous it marks a quantum critical point (QCP). In the last decade, antiferromagnetic (AF) heavy-fermion metals turned out to be model systems to study quantum criticality. In these systems, QCPs are caused by the competition between the local Kondo and the non-local Ruderman-Kittel-Kasuya-Yoshida (RKKY) interaction. Studies on the interplay between these two phenomena have revealed different types of QCPs [1]. In the itinerant, spin-density-wave (SDW) scenario the heavy fermions keep their integrity at the QCP. In such a case the QCP can be treated as a continuous classical phase transition in an effective dimension , where is the spatial dimensionality and , the dynamic exponent defined via , describes the number of additional spatial dimensions that the time dimension corresponds to. and denote the correlation length and correlation time which both diverge at a QCP [2, 3, 4]. Several heavy-fermion compounds, e. g., CeCuSi [5], CeNiGe [6] and CeLaRuSi [7] were found to exhibit this type of itinerant AF QCP.

However, in a few heavy-fermion metals the AF instability appears to be accompanied by a breakdown of the Kondo effect [8, 9, 10], which is sometimes called a zero-temperature, -electron selective Mott transition [11]. CeCuAu [12, 13] and YbRhSi [14] are prototypical examples for such an unconventional quantum critical behavior. Here, we discuss some of the unique properties of YbRhSi in the vicinity of its AF QCP. In addition, we show how the AF and Kondo breakdown QCPs can be detached under chemical pressure.

## 2 Concurrence of AF and Kondo breakdown QCPs in YbRhSi

YbRhSi crystallizes in the body-centered tetragonal ThCrSi structure. The lowest-lying of the crystal-field Kramer’s doublets is well separated from the excited states, so that at 1 K, YbRhSi is an effective spin system. It exhibits very weak AF order below mK, which is continuously suppressed by a small critical magnetic field (60 mT for field applied perpendicular to the crystallographic axis, 660 mT for ) [15] accessing the QCP. On either side of the QCP, heavy Fermi-liquid (FL) behavior was found such as a quadratic form of the resistivity where denotes the residual resistivity[16]. Approaching leads to a suppression of the FL behavior and a divergence of indicating a critical slowing down of the heavy fermions at the QCP. Right at , pronounced non-Fermi-liquid (NFL) effects are observed down to the lowest accessible temperatures ( 10 mK), notably a linear dependence of the electrical resistivity [17]. We note that the quadratic form at reported in Ref. [18] might easily arise from either incorrect fine tuning of or heating effects. The NFL form of is accompanied by a logarithmic divergence of the Sommerfeld coefficient of the electronic specific heat, , between 0.3 and 10 K. Below 0.3 K, diverges stronger than logarithmically () [15].

The two FL phases adjacent to the QCP appear to posses different Fermi surfaces as inferred from a variety of electronic transport measurements closely related to the Fermi surface properties. Most important evidence stems from Hall effect measurements in the so called crossed-field geometry. Here, two perpendicular magnets are used to disentangle the two tasks of the magnetic field: One field, , generates the Hall response, and a second field, , tunes the ground state of the sample. The power of the crossed-field setup lies in the ability to extract the initial-slope Hall coefficient as a linear response to despite measuring at a finite tuning field .

Isotherms of the field dependent Hall coefficient depicted in Fig. 1(a) show a crossover as YbRhSi is tuned across its QCP [19, 20]. The crossover is situated on top of a background contribution. The two contributions are further illustrated in the inset of Fig. 1(a) plotting the derivative . Here, the crossover corresponds to the peak in proximity to the critical field whereas the linear background is seen as an offset. Despite strong sample dependences observed in the low- Hall coefficient, the crossover appears to be robust. Rather, the sample dependences are associated with the background contribution [20]. In fact, the considered single crystals which span the hole range of sample dependences [21] show that the crossover persists in its extrapolation to zero temperature [22]. The change from a positive to a negative value of observed in this nearly compensated metal at mK was shown to be consistent with results from (renormalized) band structure calculations [21]. As the temperature is lowered, the crossover shifts to lower fields as illustrated in the phase diagram in the inset of Fig. 1(b) plotting the crossover field (see Ref. [20] for the the fitting procedure used to extract ). In the extrapolation to zero temperature converges to the QCP thus showing that the quantum criticality is the origin of the crossover in .

The crossover sharpens as the temperature is lowered. The full width at half maximum (FWHM) of the peak in is displayed in Fig. 1(b). Importantly, the FWHM is proportional to temperature, i.e., extrapolates to zero for . This finding reflects the unconventional nature of the QCP in YbRhSi. On the one hand, the vanishing width of the Hall crossover at zero temperature implies a discontinuity of the Hall coefficient and, hence, a discontinuous evolution of the Fermi surface at the QCP [22]. Such a Fermi surface reconstruction is incompatible with the smooth crossover predicted for an SDW QCP. Rather, the Fermi surface reconstruction points towards a disintegration of the quasiparticles due to the breakdown of the Kondo effect. On the other hand, the FWHM being proportional to temperature is an indicator for energy over temperature, , scaling of the single electron excitations as the width of the Fermi surface crossover and thus also of the Hall crossover can be associated with the relaxation rate of the single electron excitations [20]. Such an scaling is fundamentally inconsistent with the scaling and the concomitant superlinear temperature dependence of and the FWHM expected for a SDW QCP in three dimensions (3D). Consequently, both the fundamental signatures of a Kondo breakdown QCP, the Fermi surface collapse and the unconventional scaling behavior are seen in YbRhSi.

The crossover in is accompanied by signatures in thermodynamic and transport properties which establish a new energy scale in the phase diagram [23]. This energy scale is associated with the Kondo breakdown and for the case of stoichiometric YbRhSi at ambient pressure, it converges together with the magnetic phase boundary and the crossover to the FL regime towards the QCP (cf. inset of Fig. 1(b)).

The change of the low- Hall coefficient from hole- to electron-dominated transport (upon increasing through ) is further confirmed by low-temperature thermopower results [24]. For a change from negative to positive occurs when cooling to below = 30 mK, whereas for the results indicate negative values down to . In the latter regime the thermopower is proportional to temperature, , which is typical for a FL. In the upper inset of Fig. 2(a), the saturation of to a constant value below reflects the FL character of the ground state on the high-field side of the QCP. At all fields , exhibits a pronounced maximum succeeded by a drastic drop at lower temperatures. This anomaly at mK as well as the sign change at are unrelated to the magnetic phase boundary as illustrated in Fig. 2(b). Rather, and are constant up to where they vanish discontinuously. This is to be contrasted with the continuous nature of the magnetic transition, which is smoothly suppressed at .

In addition, the magnetic field dependence of the thermopower exhibits a drastic change upon crossing the line. The maximum observed in the isothermal thermopower coefficient, , sharpens and shifts to lower fields upon cooling. The inflection point at the low-field side of the anomaly coincides well with the energy scale , cf. lower inset in Fig. 2(a) [24]. In the limit, the field-dependence of is expected to develop into a step-like function, i. e. a discontinuity to occur exactly at the QCP, indicating an abrupt change of the Fermi surface at the QCP.

The discontinuities in as well as in may straightforwardly be related to a jump from a large Fermi surface to a smaller one when decreasing . Remarkably, in both thermopower and Hall effect isotherms the values in the magnetically ordered ground state of YbRhSi (with localized 4 states in the presence of a small Fermi surface) approach the corresponding values of the non-magnetic reference compound LuRhSi [25, 26]. In the latter, no 4 states contribute to the (small) Fermi surface either.

At elevated temperatures, a pronounced logarithmic divergence is revealed in (cf. upper inset of Fig. 2(a)). This parallels observations on [14, 15]. However, in the FL regime (), the field dependence of the saturation value extrapolates to a finite value at , where diverges. These disparate field dependences are at strong variance to the prediction within theories addressing a SDW QCP [27].

The proximity of YbRhSi to the unconventional QCP even seems to affect the classical phase transition at as here, an unusual scaling exponent is observed. In general, universal static scaling dependences are valid for all classical second-order phase transitions and are only determined by the symmetry of the order parameter and the dimensionality of the critical fluctuations [28, 29]. The high-precision specific-heat measurements on a high-quality YbRhSi single crystal depicted in Fig. 3(a) allowed for a detailed analysis of the critical fluctuations [30]. The anomaly due to the onset of magnetic order is clearly visible. The critical exponent can be extracted utilizing the standard fit function [29],

(1) |

to describe the critical behavior with the reduced temperature , where () refers to (), respectively. The background contribution is approximated by a linear dependence () close to , while the power law (first term in Eq. 1) represents the leading contribution to the singularity in .

The excellent quality of the fit becomes evident in Fig. 3(b), where the curves are shown below (black lines) and above (red lines) , together with the experimental data. For YbRhSi, the fitting procedure yields a critical exponent of , which describes the data in the entire temperature range around the phase transition at mK. Also, no saturation (due to rounding effects of the phase transition) is observed at the lowest values. The exponent derived for stoichiometric YbRhSi strongly deviates from those obtained for any universality class. The largest positive value for is calculated for a 3D-Ising system, so that generally [29]. Our surprisingly large critical exponent is supported by thermal-expansion measurements in the vicinity of . Thermodynamic relations reveal the same critical exponent for the thermal expansion as for the specific heat [31]. Application of the same fitting procedure as for the specific heat gives a critical exponent of the thermal expansion of , in good agreement with the specific-heat result. It should be noted that minor variations of the measured length change around severely complicate the critical scaling analysis and impede a more accurate determination of [30]. Nevertheless, the discrepancy of the critical exponent found in YbRhSi with known universality classes is robust and likely to originate in the nearby unconventional QCP which appears to substantially influence the spatial fluctuations of the classical order parameter.

## 3 Detaching the AF instability and Kondo breakdown by chemical pressure

Having discussed before the concurrence of the AF instability and the Kondo breakdown in stoichiometric YbRhSi the question arises whether these two instabilities could be separated by changing an additional control parameter, i.e., volume. As was found earlier, pressurizing YbRhSi stabilizes magnetism [32]. In a recent study isoelectronic substitution was utilized to generate negative and positive chemical pressure [33]. The combination with results obtained under hydrostatic pressure confirmed the dominance of the volume over the disorder effect. For both positive physical and positive chemical pressure, magnetism is strengthened as expected: and increase with decreasing unit-cell volume. Ir substitution results in a unit-cell expansion and weakens magnetism. In fact, magnetic order is completely suppressed for Yb(RhIr)Si (Fig. 4(a)) and is slightly reduced for Yb(RhIr)Si (Fig. 4 (b)). Substituting Co leads to a unit-cell compression which strengthens magnetism as illustrated for Yb(RhCo)Si in Fig. 4(c) and for Yb(RhCo)Si in Fig. 4(d). A second transition discovered for YbRhSi under hydrostatic pressure[32] is also observed for Yb(RhCo)Si at K. The surprising finding is that the Kondo breakdown energy scale is nearly insensitive against volume changes. This is obvious from Figs. 4(a) to (d) showing the line remaining at almost the same position in the - phase diagram as for pure YbRhSi (cf. inset of Fig. 1(b)).

The essential result of the chemical pressure study is summarized in the () global phase diagram: Fig. 4(e) displays the critical fields and where and extrapolate to zero temperature. The decrease of for Ir substitution reflects the decline of magnetism under lattice expansion, while the increase for Co substitution reflects the rise of magnetism under lattice compression. Clearly, remains almost unchanged in the investigated range of substitution.

The different behaviors of and have strong consequences leading to two limiting cases. For Co substituted samples we observe a crossing of and at finite . This is anticipated for 3% Co content (Fig. 4(c) and Ref. [34]). For Co content the signatures of the -line are clearly observed within the magnetically ordered phase. This finding is confirmed by a study on pure YbRhSi under physical pressure [35]. Also no relation between the second transition at and the line is evident as only for Yb(RhCo)Si vanishes in proximity to line rendering this an accidental coincidence. With assigned to the Kondo breakdown, we infer from that the magnetism extends into the regime where the Kondo effect is operating. Hence, the magnetic instability is expected to obey the signatures of an SDW QCP. This is in accordance with the observed power-law dependence of , yielding an exponent in surprisingly good agreement with predicted for 3D SDW fluctuations [33].

In addition, the critical exponent of the specific heat recovers a conventional value for Yb(RhCo)Si. Fig. 3(a) includes specific heat data on two selected Co-concentrations ( and 0.27). Already the shape of the anomalies indicates a change of the critical fluctuations with increasing . Going from to , the phase transition anomaly at changes drastically, i.e., from a very sharp -type peak to a rather broad mean-field-like anomaly detected for concentrations (not shown) [36], which prevents the extraction of critical fluctuations. For one again observes a pronounced -type anomaly which allows to analyze the data in terms of critical fluctuations.

In Fig. 3(c), we present the analysis for sample. Here, K exceeds of pure YbRhSi by one order of magnitude. The best fit for reveals a negative critical exponent , which can already be inferred from the non-diverging -dependence in Fig. 3(c) (for comparison see Fig. 1 in Ref. [29]). To fit the data points at we have to allow for a Gaussian distribution of with mK to account for the rounding effects due to the high substitution level. This smearing of leads to a relatively small temperature range, , which determines the critical exponent, impeding a more accurate determination of . A similar measurement and analysis of the critical behavior was also carried out for , yielding again a negative critical exponent [37].

Comparing the critical exponents for the different materials, it is obvious that there is a drastic change on going from for to for . The latter value can be explained in terms of the classical universality classes in the Landau theory of phase transitions, for which generally holds true. However, the accuracy of the determined exponent for the samples with substitution is not sufficient to distinguish between the two applicable symmetry classes [29], namely the 3D-Heisenberg model or the 3D-XY model .

The evolution of the phase transition anomaly at in can be understood in the frame of Doniach’s phase diagram applied to Yb compounds [1]. At large Co concentrations (small average unit-cell volume) where the RKKY interaction predominates, classical ordering between local 4-derived magnetic moments occurs. Upon lowering the Kondo effect comes into play. This leads to mean-field type itinerant AF (SDW) order in the renormalized electron fluid. However, as the Kondo effect breaks down and the local moments reappear – fully unexpected within the Doniach phase diagram. Because of the small value of , the nearby quantum critical electronic fluctuations associated with this Kondo breakdown QCP add to the classical critical AF order parameter fluctuations and very likely contribute to the unusually large critical exponent found for pure YbRhSi.

Furthermore, our result of a conventional critical exponent for chemically pressurized YbRhSi is in contradiction to the quantum-tricritical-point scenario [38]. Here, Misawa et al. suggest a stabilized tricritical point at finite temperatures in pressurized YbRhSi. As the theoretical value at a classical tricritical point is , one would expect that the critical exponent becomes even larger for YbRhSi under chemical pressure, just opposite to our observation.

The findings on Co substituted YbRhSi are resembled by two other heavy-fermion systems: Both CeIn and CeRhCoIn display a Fermi surface reconstruction inside the antiferromagnetically ordered state [39, 40]. For the case of CeIn this suggests a transition from a large Fermi surface with the -electrons incorporated to a small Fermi surface with the electrons decoupled from the conduction sea. For CeRhCoIn the Fermi surface reconstruction appears to be accompanied by a change in the magnetic structure.

We now turn to the samples with Ir substitution, at negative chemical pressure where . Thus, for instance, for 6% Ir substituted YbRhSi mT is well below the Kondo-breakdown critical field mT. Consequently, there exists a finite field range at (indicated by question marks in Fig. 4(e)) in which the magnetic moments are neither ordered nor screened by the Kondo effect. The electrical resistivity exhibits pronounced NFL behavior in this regime which was tentatively ascribed to a metallic spin-liquid phase, hitherto unknown in a Kondo-lattice system [33]. Similar observations were made with other lattice expanded YbRhSi samples, substituted either with La on the Yb site or Ge on the Si site[41, 42].

The most prominent point in the global phase diagram (Fig. 4(e)) is taken by stoichiometric YbRhSi. Here, the AF instability and the Kondo breakdown concur, . Within experimental resolution this concurrence extends towards small Ir content as obvious from Fig. 4(b) [43]. For Yb(RhIr)Si the critical field of the Néel phase is slightly reduced compared to pure YbRhSi. Apparently, the critical field of the -line seems to be reduced by almost the same amount. Consequently, and concur within experimental accuracy for this slightly Ir substituted sample. Such a concurrence between the AF and electronic instabilities in a finite parameter range of the global phase diagram for a Kondo lattice has been predicted theoretically [44, 45]

## 4 Conclusion and Perspective

We have provided compelling evidence for a breakdown of the Kondo effect overlapping with the field induced AF QCP in pure YbRhSi. This is found to be accompanied by a reconstruction of the Fermi surface and highly unusual scaling behavior. Upon varying the average unit-cell volume of YbRhSi by partial isoelectronic substitution these very different instabilities can be separated from each other. Volume compression (Co substitution) furnishes an intersection of the Kondo-breakdown crossover line with the magnetic phase boundary . Here, the AF QCP () is of conventional 3D SDW type. Future studies on moderately Co substituted samples can establish the existence of a quantum phase transition between local and itinerant AF order.

When the unit cell of YbRhSi is moderately expanded, a novel low-temperature NFL phase, presumably of spin-liquid type, develops in a wide field range. Whether such a phase can in fact arise within the Kondo lattice model is an intriguing open theoretical question. Finally, YbRhSi with sufficiently expanded volume (e.g. containing more than 10% Ir) should be investigated to unravel the nature of the Kondo breakdown QCP without any interfering magnetism. At small Ir concentrations a finite range seems to exist where the magnetic and electronic instabilities concur. Specific heat measurements (e.g. on samples with % Ir) may reveal whether the unusual scaling behavior is linked to this concurrence.

###### Acknowledgements.

We greatly acknowledge valuable discussions with P. Coleman, G. Donath, P. Gegenwart, C. Klingner, M. Nicklas, N. Oeschler, S. Paschen and G. Zwicknagl. This works was supported by the DFG Research Group 960 “Quantum Phase Transitions” as well as by the NSF Grant No. DMR-1006985 and the Robert Welch Foundation Grant No. C-1411.## References

- [1] P. Gegenwart, Q. Si, and F. Steglich: Nature Phys. 4 (2008) 186.
- [2] J. A. Hertz: Phys. Rev. B 14 (1976) 1165.
- [3] T. Moriya: Spin Fluctuations in Itinerant Electron Magnetism (Springer, Berlin, 1985).
- [4] A. J. Millis: Phys. Rev. B 48 (1993) 7183.
- [5] P. Gegenwart, C. Langhammer, C. Geibel, R. Helfrich, M. Lang, G. Sparn, F. Steglich, R. Horn, L. Donnevert, A. Link, and W. Assmus: Phys. Rev. Lett. 81 (1998) 1501.
- [6] P. Gegenwart, F. Kromer, M. Lang, G. Sparn, C. Geibel, and F. Steglich: Phys. Rev. Lett. 82 (1999) 1293.
- [7] W. Knafo, S. Raymond, P. Lejay, and J. Flouquet: Nature Phys. 5 (2009) 753 .
- [8] Q. Si, M. S. Rabello, K. Ingersent, and J. L. Smith: Nature 413 (2001) 804.
- [9] P. Coleman, C. Pépin, Q. Si, and R. Ramazashvili: J. Phys.: Condens. Matter 13 (2001) R723.
- [10] T. Senthil, M. Vojta, and S. Sachdev: Phys. Rev. B 69 (2004) 035111.
- [11] I. Paul, C. Pepin, and M. R. Norman: Phys. Rev. Lett. 98 (2007) 026402.
- [12] H. v. Löhneysen, T. Pietrus, G. Portisch, H. G. Schlager, A. Schröder, M. Sieck, and T. Trappmann: Phys. Rev. Lett. 72 (1994) 3262.
- [13] A. Schröder, G. Aeppli, R. Coldea, M. Adams, O. Stockert, H. Lohneysen, E. Bucher, R. Ramazashvili, and P. Coleman: Nature 407 (2000) 351.
- [14] O. Trovarelli, C. Geibel, S. Mederle, C. Langhammer, F. M. Grosche, P. Gegenwart, M. Lang, G. Sparn, and F. Steglich: Phys. Rev. Lett. 85 (2000) 626.
- [15] J. Custers, P. Gegenwart, H. Wilhelm, K. Neumaier, Y. Tokiwa, O. Trovarelli, C. Geibel, F. Steglich, C. Pépin, and P. Coleman: Nature 424 (2003) 524.
- [16] P. Gegenwart, J. Custers, C. Geibel, K. Neumaier, T. Tayama, K. Tenya, O. Trovarelli, and F. Steglich: Phys. Rev. Lett. 89 (2002) 056402.
- [17] P. Gegenwart, T. Westerkamp, C. Krellner, M. Brando, Y. Tokiwa, C. Geibel, and F. Steglich: Physica B: Condensed Matter 403 (2008) 1184.
- [18] G. Knebel, R. Boursier, E. Hassinger, G. Lapertot, P. G. Niklowitz, A. Pourret, B. Salce, J. P. Sanchez, I. Sheikin, P. Bonville, H. Harima, and J. Flouquet: J. Phys. Soc. Jpn. 75 (2006) 114709.
- [19] S. Paschen, T. Lühmann, S. Wirth, P. Gegenwart, O. Trovarelli, C. Geibel, F. Steglich, P. Coleman, and Q. Si: Nature 432 (2004) 881.
- [20] S. Friedemann, N. Oeschler, S. Wirth, C. Krellner, C. Geibel, F. Steglich, S. Paschen, S. Kirchner, and Q. Si: Proc. Natl. Acad. Sci. USA 107 (2010) 14547.
- [21] S. Friedemann, S. Wirth, N. Oeschler, C. Krellner, C. Geibel, F. Steglich, S. MaQuilon, Z. Fisk, S. Paschen, and G. Zwicknagl: Phys. Rev. B 82 (2010) 035103.
- [22] S. Friedemann, N. Oeschler, S. Wirth, C. Krellner, C. Geibel, F. Steglich, S. Paschen, S. Kirchner, and Q. Si: to be published .
- [23] P. Gegenwart, T. Westerkamp, C. Krellner, Y. Tokiwa, S. Paschen, C. Geibel, F. Steglich, E. Abrahams, and Q. Si: Science 315 (2007) 969.
- [24] S. Hartmann, N. Oeschler, C. Krellner, C. Geibel, S. Paschen, and F. Steglich: Phys. Rev. Lett. 104 (2010) 096401.
- [25] U. Köhler, N. Oeschler, F. Steglich, S. Maquilon, and Z. Fisk: Phys. Rev. B 77 (2008) 104412.
- [26] S. Friedemann, N. Oeschler, C. Krellner, C. Geibel, S. Wirth, F. Steglich, S. Paschen, S. MaQuilon, and Z. Fisk: Physica B 403 (2008) 1251.
- [27] I. Paul and G. Kotliar: Phys. Rev. B 64 (2001) 184414.
- [28] L. D. Landau: Phys. Z. Sowjetunion 11 (1937) 26.
- [29] J. Wosnitza: J. Low Temp. Phys. 147 (2007) 249.
- [30] C. Krellner, S. Hartmann, A. Pikul, N. Oeschler, J. G. Donath, C. Geibel, F. Steglich, and J. Wosnitza: Phys. Rev. Lett. 102 (2009) 196402.
- [31] E. Scheer, J. Wosnitza, H. v. Löhneysen, R. Kürsch, M. Lang, and F. Steglich: J. Magn. Magn. Mat. 104-107 (1992) 175.
- [32] S. Mederle, R. Borth, C. Geibel, F. M. Grosche, G. Sparn, O. Trovarelli, and F. Steglich: J. Phys.: Condens. Matter 14 (2002) 10731.
- [33] S. Friedemann, T. Westerkamp, M. Brando, N. Oeschler, S. Wirth, P. Gegenwart, C. Krellner, C. Geibel, and F. Steglich: Nature Phys. 5 (2009) 465.
- [34] S. Friedemann, N. Oeschler, C. Krellner, C. Geibel, and F. Steglich: J. Phys.: Conf. Ser. 100 (2010) 012038.
- [35] Y. Tokiwa, P. Gegenwart, C. Geibel, and F. Steglich: J. Phys. Soc. Jpn. 78 (2009) 123708.
- [36] A. Steppke, private communication, 2010.
- [37] C. Krellner, C. Klingner, C. Geibel, and F. Steglich: Phys. Status Solidi B 247 (2010) 734.
- [38] T. Misawa, Y. Yamaji, and M. Imada: J. Phys. Soc. Jpn. 78 (2009) 084707.
- [39] N. Harrison, S. E. Sebastian, C. H. Mielke, A. Paris, M. J. Gordon, C. A. Swenson, D. G. Rickel, M. D. Pacheco, P. F. Ruminer, J. B. Schillig, J. R. Sims, A. H. Lacerda, M.-T. Suzuki, H. Harima, and T. Ebihara: Phys. Rev. Lett. 99 (2007) 056401.
- [40] S. K. Goh, J. Paglione, M. Sutherland, E. C. T. O’Farrell, C. Bergemann, T. A. Sayles, and M. B. Maple: Phys. Rev. Lett. 101 (2008) 056402.
- [41] F. Weickert, P. Gegenwart, J. Ferstl, C. Geibel, and F. Steglich: Physica B 378-380 (2006) 72.
- [42] J. Custers, P. Gegenwart, C. Geibel, F. Steglich, P. Coleman, and S. Paschen: Phys. Rev. Lett. 104 (2010) 186402.
- [43] S. Friedemann, T. Westerkamp, M. Brando, N. Oeschler, P. Gegenwart, C. Krellner, C. Geibel, S. MaQuilon, Z. Fisk, F. Steglich, and S. Wirth: J. Low Temp. Phys. 161 (2010) 67.
- [44] Q. Si: Phys. Status Solidi B 247 (2010) 476.
- [45] Q. Si: Physica B 378-380 (2006) 23.