Microsoft Word - PHYS081124-f edited_checked.doc 71-83 SQU Journal For Science, 14 (2009) © 2009 Sultan Qaboos University 71 Anisotropy in RFe2 Intermetallics: An NMR Study C. Carboni Department of Physics, College of Science, Sultan Qaboos University, P.O. Box 36, Postal code 123, AlKhod, Muscat, Sultanate of Oman, Email: carlo@squ.edu.om. NMR دراسة : RFe2تباين خواص المركبات المعدنية كاربوني. ك باسـتخدام اعتمـاد طيـف الـرنين . 3oH Tb3++3mT+. النادرةاألرضلقد إستقصينا تباين خواص أيونات :خالصـة 2TmFe 2TbF 2HoFe الموجودة في المركبات mT169 , bT159, Ho165المغناطيسي النووي على شدة المجال لنوى تفاعل الدقيق لهذه األيونات على طيف الرنين المغناطيسي النووي و الذي يعطى شـدة يطغى ال . و المركبات ذات الصلة .كما قمنا بمقارنة القياسات بالحسابات النظرية. جيجاهيرتز7 الى 2االستقطاب في حيز التردد من ABSTRACT: The anisotropy of the Ho3+, Tb3+ and Tm3+ rare-earth ions in RFe2 compounds is investigated using the field dependence of the NMR spectra of 169Tm, 159Tb and 165Ho in polycrystalline TmFe2, TbFe2, HoFe2 and related pseudo-binary compounds. The very large hyperfine interaction in these ions dominates the NMR spectrum; the NMR frequencies in the 2 to 7 GHz range provide therefore a measure of the localized 4 f moment and its orientation in the crystal field. We compare the NMR measurements to single ion computations that include magnetostriction. KEYWORDS: Cubic laves phases, Rare-earth intermetallic, Hyperfine splitting, Magnetic anisotropy, NMR. 1. Introduction he cubic Laves phases series RFe2 where R is a rare-earth element is the simplest of the several inter- metallic compounds in which rare-earth and iron magnetism coexist. For more than 30 years the RFe2 compounds and their pseudo-binary (RR') Fe2 and R(FeM)2 alloys where M is a transition metal have been the subject of intense investigations. These compounds have provided and continue to provide (Sawickki et al, 2000) a fertile ground to test theoretical models. More recently epitaxially grown films (Lee et al, 2000, de la Fuente et al, 2001, Gordeev et al, 2001, Oster et al, 2005), X-ray magnetic circular dichroism investigations (Fujiwara et al, 2005, Dumesnil et al, 2002, Parlebas et al, 2006) and the hydrated RFe2 materials (Paul-Boncour, 2004) have attracted the attention of several researchers. However, there are some basic questions in the RFe2 that are not yet completely resolved. There is still uncertainty about the anisotropy of the iron sub-lattice, about the T C. CARBONI 72 contribution of magneto-elastic effects to the measured anisotropy (Mougin et al, 2000), about the crystal-field interaction at the lanthanide site and about the role of the ordered super-lattice of vacancies in the C15 Laves phases (Gratz et al, 1996). Also, the distribution of the total magnetization between the iron and lanthanide sub- lattices and the relative contribution of localized and itinerant moments is still an open question (Fujiwara et al 2005). In the present paper the anisotropy of the rare-earth ion is investigated using the field dependence of the NMR spectra of 169Tm(I = 1/2), 159Tb(I = 3/2) and 165Ho(I = 7/2) in polycrystalline TmFe2, TbFe2, HoFe2 and some pseudo-binary compounds. In the Ho3+, Tb3+ and Tm3+ ions the very large hyperfine interaction dominates the NMR spectrum; the NMR frequencies in the 2 to 7 GHz range provide therefore a measure of the localized 4 f moment and its orientation in the crystal field. We compare our NMR measurements to single-ion computations that include the distortion of the cubic environment of the lanthanide ion due to the magnetostriction. The results from the computation suggest that the magnetostriction does not have a significant effect on the NMR spectrum. 2. Theory In this section we describe the model used for the computation of the magnitude and the orientation of the moments in RFe2 compounds subjected to a magnetic field applied in any arbitrary direction but before doing so, we need to clarify a possible confusion between two senses of the term “anisotropy”. There are two distinct aspects to the anisotropy of a magnetic material, the energy anisotropy and the magnetization anisotropy. The former, conventionally specified by classical “anisotropy constants” denotes the angular variation of the free energy, whereas the latter denotes the angular variation of the magnitude of the ionic moment. For the RFe2 compounds in which the exchange interaction exH may, at least to a first approximation, be considered isotropic both types of anisotropy are manifestations of the crystal-field interaction cfH at the lanthanide site. The relationship between energy and magnetization anisotropy is not straightforward. For example, it is not necessarily the case that the hardest direction of magnetization (i.e. the direction for which the free energy is greatest) is the direction in which the ionic moment is most quenched. The crystal-field interaction cfH in the 2RFe is a small perturbation on the exchange interaction exH the energy anisotropy may thus be regarded as a quasi-classical effect of first order in cfH . The magnetization anisotropy, on the other hand, is an essentially quantum mechanical effect associated with those terms in cfH that do not commute with exH . The quenching of the ionic moment is a second order effect varying as the ratio of square of the strength of cfH to that of exH . A fuller discussion of these points can be found in section 4.4 of McCausland and Makenzie (1979). 2.1 The effective electronic Hamiltonian The orientation of the lanthanide and the iron moments in a crystal of RFe2 subjected to a magnetic field applied at an arbitrary angle is a compromise between several competing interactions. There is the ferrimagnetic coupling between the lanthanide and the iron, the crystal-field interaction, the interaction with the applied field and indirectly the magnetostriction which affects the crystal field. If one assumes that the iron sub-lattice is isotropic in both senses, the magnitude of the iron moment is constant and its direction always collinear to the total field FetotB at the iron site, the equilibrium orientation of the two moments in the presence of a field applied at an arbitrary direction is as shown on Figure 1. The total field at the lanthanide site is mBB λ−= a ln tot where m is the iron moment and λ the lanthanide-iron exchange constant. ANISOTROPY IN RFe2 INTERMETALLICS 73 The lanthanide ion is described by the effective electronic Hamiltonian Figure 1. The orientation with respect to the direction of the applied field Ba of the iron moment m, the lanthanide moment µ, the total field lntotB at the lanhanide site and the total field Fe totB at the iron site. Fe exB and Fe exB are respectively the exchange field seen by the iron and the exchange field seen by the lanthanide ion. µBVOBOB lntot26644el ⋅−++=H (1) where the n nB O denote the cubic crystal-field terms in the notation of Lea, Leask and Wolf (1962). The operators 4O and 6O are respectively combinations of fourth order and sixth order Stevens operators equivalent (Stevens 1952); the explicit form of these operators is given in the appendix. The term 2V is the second order crystal-field term arising from the distortion of the cubic environment of the ion caused by magnetostriction (Buschow 1977). Ba Ba m µ ln exB ln totB Fe exB Fe totB D ir ec ti on o f t he ap pl ie d fi el d C. CARBONI 74 [ ])(kk)(kk)(kk )1J2(J K2 311313233232122121 me 2 JJJJJJJJJJJJV +++++− = (2) where K is a unit vector parallel to the total angular momentum J of the 4 f electrons of the ion and meK is the magneto elastic constant. The indices 1, 2 and 3 denote respectively the directions of the three four-fold axis of the cubic unit cell. The crystal-field parameters are phenomenological parameters determined experimentally. In transparent compounds these parameters can be accurately obtained from the optical absorption spectrum. In metals the determination of the crystal field parameters is not straightforward; they are obtained by combinations of different experimental techniques which can produce conflicting results (Mcmorrow et al, 1989). For the RFe2 series there is to our knowledge no detailed comprehensive analysis for the crystal field parameters. To date the most widely accepted sets of crystal field parameters for the RFe2 are the ones obtained by Germano and Butera (1981). The parameters for the Tb, Ho and Tm compounds are given in Table 1 together with the exchange and the magnetoelastic constant. Germano and Butera obtained the crystal field parameters from specific heat measurements. In their analysis the magenostriction is not taken into account, the site symmetry is considered exactly cubic; this plus the fact that the contribution from the crystal field to the specific heat is relatively small may raise questions concerning the accuracy of the parameters. Nevertheless most of the predictions based on these parameters are in fairly good agreement with experiments; these parameters are used faute de mieu for the present work. The computed direction of spontaneous magnetization in HoFe2, TmFe2 and TbFe2 at zero temperatures is respectively the <100>, <110> and <111> direction. Table 1. The exchange fields, the crystal-field parameters and the magnetoelestic constants obtained by Germano and Butera (1981) for HoFe2, TbFe2 and TmFe2. λm (T) 104 B4 (K) 10 6B6 (K) Kme (K/ion) TbFe2 162 101 20 -13 HoFe2 108 -9.26 5.1 -1.7 TmFe2 70 55 35 -9.0 In TmFe2 the computation reveals a low-lying first excited state that strongly prefers the <111> direction. The existence of this excited state is responsible for the reorientation of the magnetization at about 55 K reported by Germano and Butera. Above 55 K the <111> direction is the preferred orientation for the magnetization. The experiments reported in the present work were carried out at liquid-helium temperatures where only the ground state is significantly populated therefore the <110> direction is strongly preferred in TmFe2. 2.2 Computation In zero applied field and in the absence of magnetostriction the ground state and ground state energy of the lanthanide ion can be readily obtained form the Hamiltonian 1. The expectation value of zJ for the ground state is then obtained to calculate the intra-ionic contribution to the hyperfine splitting and the rare-earth’s moment. In the presence of an applied field at an arbitrary orientation one must first determine the orientation of ln totB . To determine the equilibrium orientation of lntotB in an applied field we have minimized with respect to the polar angles of the iron moment, Feθ and Feφ the free energy per formula unit at absolute zero including the magnetostriction. At zero temperature the free energy is simply ANISOTROPY IN RFe2 INTERMETALLICS 75 , (3)el FeF U TS U H H= − = = + where ( )Fe aH 2 λ= −B m.µ is the contribution from the iron sub-lattice and the expectation values are calculated using the ground state of the lanthanide ion. In writing equation (3) it is implicitly assumed that the iron sub-lattice is isotropic in both senses. The minimization of the free energy (3) is not straightforward. The term 2V depends on the orientation of the lanthanide moment that is itself the object of the computation. The computation was done by iterations. An educated guess of the orientation of the lanthanide moment was used for the first iteration. The computation was repeated until self-consistency and minimum energy was achieved. It was observed that the term 2V does not affect significantly the outcome of the computation. In order to investigate the behavior of polycrystalline samples in applied fields we have computed the magnitude and orientation of the lanthanide moment in HoFe2 , TbFe2 and TmFe2 at various fields applied in 30 different orientations equally spaced over the portion of the unit sphere defined by the <100>, <111> and <110> directions. The results for the 8 T field are summarized in Table 2 and Figure 2. Figure 2 is a planar projection of the portion of the unit sphere under consideration. The orientation of the moment is represented on the Figure by a point corresponding to the intersection of the direction of the moment and the surface of the unit sphere. The shaded areas around the <100>, <110> and <111> axis represent the loci of the direction of the holmium, the thulium and the terbium moments respectively, when the direction of the applied field is scanned across the entire area of the unit sphere under consideration. One can see that an 8 T field applied in any direction hardly moves the thulium moment away from its <110> preferred orientation, 5 degrees at most when the field is along the <111> direction. Holmium, on the other hand, can be more readily moved away from the <100> preferred orientation, 25 degrees when the field is along the <110> axis. Table 2. The computed orientation of the iron’s magnetic moment m, the lanthanide magnetic moment µ and the total field lntotB at the lanthanide site for an 8 tesla field applied along each of the principal crystallographic directions. The column before last gives the computed intra-ionic contribution to the dipolar hyperfine splitting and the last column gives the total field-dependent contribution (intra-ionic plus the contribution from applied field). ψ is the angle between the applied field and the lanthanide moment. Ba m lntotB µ Ion θ φ θ φ θ φ θ φ a’ = az (MHz) a’+ γΒcosψ (MHz) 0 0 54.0 45 52 45 52 45 3171 3219 45 0 54.2 44 53 43 54 43 3171 3236 Tb3+ 54.7 45 54.7 45 54.7 45 54.7 45 3171 3250 0 0 0 0 0 0 0 0 6488 6559 45 0 25 0 27 0 26 0 6479 6546 Ho3+ 54.7 45 13 45 15 45 13.5 45 6484 6537 0 0 48 0 44 0 44 0 2308 2328 45 0 45 0 45 0 45 0 2309 2337 Tm3+ 54.7 45 45 0 45 5 45 4 2308 2331 C. CARBONI 76 30 10 20 30 40 40 2010 <110> <111> <100> TbFe2 HoFe2 TmFe2 Figure 2. The planar projection of the portion of the unit sphere under consideration. The angles along both axes are in degrees. The shaded areas around the <100>, <110> and <111> axis represent the loci of the direction of the holmium, the thulium and the terbium moments respectively, when the direction of an applied 8 T field is scanned across the entire area of the unit sphere under consideration. Table 2 gives the computed polar angles of the iron moment m , the lanthanide moment µ and the total field at the lanthanide site lntotB when the 8 tesla field aB is applied along each of the principal crystallographic orientations. The last two columns give computed field dependent contributions to the dipolar hyperfine splitting which will be discussed in the next section. 2.3 The hyperfine splitting The theory of the hyperfine splitting in rare-earth compounds is described in detail in McCausland and Makenzie (1979), Bunbury et al (1989) and Li et al (1996) whose notation we shall follow. The hyperfine splitting of the electronic states is described to first order by the effective nuclear spin Hamiltonian { })II(IH 22zz 31tt Pah −+= (4) where I is the nuclear spin, ta at and tP are respectively the dipolar and the quadrupolar hyperfine parameters. The z-axis is chosen to be parallel to J . Preliminary computations have shown that in the exchange dominated RFe2 compounds second order contributions to the hyperfine parameters are at most 0.05%. Therefore, the second order terms (Bunbury et al 1989) have been neglected in the present work. Because of the strong hyperfine coupling in Ho3+, Tb3+ and Tm3+ the total dipolar hyperfine constant is dominated by the dipolar intra- ionic term (5) ' za a J= ANISOTROPY IN RFe2 INTERMETALLICS 77 where a, the dipolar hyperfine coupling constant, is 6497 MHz for Ho3+, 3180 MHz for Tb3+ and 2333 MHz for Tm3+, the intra-ionic contribution a' is therefore very sensitive to small variations of the 4 f moment and provides an accurate measure of the moment. The extra-ionic contribution to the dipolar splitting is zB n nga h µ⎛ ⎞′′ ′′= ⎜ ⎟ ⎝ ⎠ (6) where zB′′ is the component along J of B′′ the extra-ionic field at the lanthanide nucleus. ′′B is the vector sum of various terms a dip p n′′ = + + +B B B B B (7) where aB is the applied field, dipB is the dipolar field, pB is the field due to conduction electrons polarized by the spin of the parent ion and nB is the transferred hyperfine field from neighboring rare-earth and iron ions. To a good approximation ( )dip p n+ +B B B can be regarded as independent of the applied field. Since we are interested in the field dependence of the hyperfine splitting, we will consider the field independent part of the extra-ionic contribution to the dipolar hyperfine splitting ( ) ( )thf n n dip p na g / h B B B= µ + + (8) as a free parameter to be determined from the zero-field NMR spectrum; thfa is equal to the difference between the computed intraionic contribution and the measured total dipolar hyperfine parameter in zero fields. The total dipolar hyperfine splitting is therefore t thfa a a B cosγ′= + + ψ (9) where ψ is the angle between the applied field and the lanthanide moment and γ is the gyromagnetic factor. The field dependent contribution ( )a B cosγ′ + ψ computed in a field of 8 T applied along each of the three principal crystallographic axis is given in Table 2. 3. Experiment. The compounds were prepared by melting the constituents in an arc furnace under an atmosphere of purified argon. For the thulium compounds, the starting materials contained a 3% excess of Tm to compensate for the anticipated loss by evaporation. The resulting ingots were re-melted three times, annealed for two weeks at 900 0C and finally crushed to form a fine powder sample. X-ray examination revealed less than 5% of non- Laves phases. Approximately 5 mg of the powder, embedded in epoxy resin, was located close to the current antinode of the central conductor of a tunable coaxial resonator of the type described by McCausland and Makenzie (1979). The NMR measurements were carried out at 1.3 K and 4.2 K in magnetic fields up to 8 T using the microwave spin-echo spectrometer described by Carboni et al (1989). The spectra were obtained by varying the frequency at fixed values of the applied field. C. CARBONI 78 Figure 3. The 169Tm NMR line of thulium in TmFe2 in zero fields and in an applied field of 7.5 T. For TmFe2 and TbFe2 the measured relaxation times are of the order of several micro seconds which is an order of magnitude larger than the pulse separation (200 ns) used to take the spectra so the NMR lines are negligibly distorted by relaxation. In pure Ho Fe2 the transverse relaxation time is of the order of 200 ns, too short to obtain a reliable spectrum. The seven lines of the 165Ho (I = 7/2) spectrum could not be unequivocally identified. The 165Ho NMR data presented here were obtained mainly from specimens with a few percent of holmium substituted in TmFe2, LuFe2 and YFe2 where the relaxation times are less prohibitive and the spectra in low fields are resolved. The spectrum of 159Tb(I = 3/2) consists of three lines; because of the strong quadrupolar hyperfine coupling the lines are very well resolved; the separation between the lines is of the order of 350 MHz. The spectrum of 169Tm(I = 1/2) consists of a single line. A summary of the hyperfine parameters measured in the various compounds investigated at liquid helium temperatures and in zero fields is given in Table 3. For comparison previously published data from other R': RFe2 materials have also been included in the table. 4. Discussion. The NMR line of 169Tm:TmFe2, in zero fields and in an applied field of 7.5 tesla are shown in Figure 3. The field dependence of the measured total hyperfine parameter ta of 169Tm:TmFe2, is shown in Figure 4-a; the variation is linear with a slope of (3.42 + 0.05) MHz/T. The slope is less than the value of 3.52 MHz/T expected for the nuclear spin parallel to the applied field. We also observe that the applied field causes only a small broadening of the NMR line (see Figure 5). ANISOTROPY IN RFe2 INTERMETALLICS 79 A B TbFe2 T = 4.2 K Figure 4. The measured total dipolar hyperfine parameter of 169Tm:TmFe2 and 159Tb:TbFe2 as a function of the applied field. These two observations confirm the prediction from our computations that the moment of the Tm3+ ion in all the crystallites of the specimen has remained close to the <110> orientation even at high fields. The computation C. CARBONI 80 shows that there is a significant quenching of if the moment is rotated away from the <110> orientation with a concomitant reduction of the hyperfine splitting. The reduction is 28% if the thulium moment is forced along the <100> direction corresponding to a shift of more than 350 MHz in the NMR frequency. Therefore, if the rare-earth moments rotate to become close to the direction of the applied field one expects a line broadening of the order of 200 to 300 MHz. This is not observed. When the thulium ion is placed in a non magnetic and isotropic host like 1% 169Tm (Tm:LuFe2) we observe that at decreases with applied field with a slope of (–13 + 2) MHz/T, much larger than the -3.52 MHz/T expected for the nuclear spin anti-parallel to the applied field. The line broadens rapidly with the applied field; the line width is 10 MHz in a field of 0.7 T and the line is too broad for reliable measurements in a field of 1 T. Assuming that the lutetium sub-lattice does not contribute significantly to the anisotropy these observations can be interpreted as follows: the applied fields causes a significant rotation of the iron moments towards a direction parallel to the applied field. Because of the strong ferimagnetic coupling between the thulium ion and the iron sublattice, the thulium moments are gradually rotated away from the strongly preferred <110> orientation hence the sharp reduction in the intraionic contribution to at and the line broadening. Finally we note from table 3 that when the thulium ion is placed in the HoFe2 host there is an increase of about 10 MHz in the dipolar hyperfine parameter and not the expected reduction of about 350 MHz if the thulium moment was forced along the <100> orientation preferred by the HoFe2 host. One can conclude that thulium moment has remained close to the <110> direction and that therefore the iron together with the holmium sublattice are not along the <100> direction. The increase in the dipolar hyperfine parameter is due to the larger transferred hyperfine field in the holmium host due the fact that Ho3+ has a larger projected spin (σHo = 2.0) than Tm 3+ (σTm = 1.0) (McCausland and Makenzie, 1979). A refined analysis of the data presented in this section is still in progress and will be presented in a separate publication (Carboni et al 2009, to be published). Table 3. The measured hyperfine parameters of Tm3+, Ho3+ and Ho3+ in various RFe2 hosts. (a) Li et al (1996), (b) De Azevedo (1986). (c) Al Assadi et al (1984). Nucleus and host at (MHz) Pt (MHz) 169Tm:TmFe2 2 +2456 - 1% 169Tm (Tm:Ho Fe2) 2 +2467 - 1% 169Tm (Tm:Lu Fe2) 2 +2445 - 169Tm (Tm:Y Fe2) (a) 2 +2422 - 169Tm (Tm:Gd Fe2) (a) 2 +2466 - 1% 165Ho (Ho:TmFe2) 5 +6920 1 +60 3% 165Ho (Ho:LuFe2) 10 +6850 2 +57 165Ho (Ho:YFe2) (c) 10 +6860 2 +59 165Ho (Ho:GdFe2) (c) 10 +6986 159Tb:TbFe2 3 +3798 2 +5 34 159Tb(Tb:LuFe2) (b) 3 +3693 2 +349 159Tb(Tb:GdFe2) (b) 3 +3822 2 +356 The field dependence of the measured total dipolar hyperfine parameter at of 159Tb:TbFe2, is shown on Figure 4- b; the variation is linear with a slope of (9.5 + 0.1) MHz/T. The field dependence of the line-width is given in Figure 5. The data shows that like for TmFe2 the anisotropy is, as predicted by the computation, very large and the terbium ion remains close to the <111> direction even at high fields. A detailed investigation of the zero field hyperfine splitting of Tb3+ in various RFe2 hosts has been done by W M de Azevedo (de Azevedo, 1986 and de Azevedo et al 1985). Except in the TmFe2 host, the NMR spectra from 165Ho in the RFe2 hosts are poorly resolved because the transverse relaxation time is comparable to the minimum pulse separation available on the spectrometer. ANISOTROPY IN RFe2 INTERMETALLICS 81 Moreover, the NMR lines broaden rapidly with the applied field hence the large uncertainty on the hyperfine parameters for the holmium given in table 3. In the HoFe2 compound it was not possible to obtain the complete 7-lines spectrum therefore it was not possible to obtain a value for the dipolar hyperfine parameter. The quadrupolar parameter Pt obtained from the spacing between adjacent lines for this compound is (55 + 5) MHz. Figure 5. The behavior of the line width of the central line of the NMR spectra of 169Tm:TmFe2, 159Tb:TbFe2 and 165Ho: Ho(TmFe2). For holmium in the TmFe2 host the spectrum is well resolved in fields up to 5 T. However, as can be seen in Figure 5, the line-width increases sharply with the applied field. If one assumes that the line broadening is mainly due to a spread of the orientations of the Ho3+ moments with respect to the crystallites axis then, since it has been established that the Tm3+ moment remains along the <110> direction, one can conclude that the applied field does, as predicted by the computation, achieve some rotation of holmium moment away from the <100> direction. The slope of at versus field is (8.3 + 0.2) MHz/T, significantly less than the expected 8.9 MHz/T; this is almost certainly due to a gradual quenching of the holmium moment as it is rotated away from its preferred direction. 5. Acknowledgements The experimental part of this work was done during the time the author was at the University of Manchester (UK) supported by grants from the UK Engineering and Physical Science Research council (EPSRC). The author is also grateful to Prof MAH McCausland and Dr D StP Bunbury for comments, discussions and support. C. CARBONI 82 Appendix The Steven operator equivalent pertaining to a site of cubic symmetry. [ ] 222404 )1(3)1(625)1(3035 +++−−+−= JJJJJJJJO zz ( )44 2 14 4 JJO −+ += [ ] [ ] )1((60)1(40)1(5 294)1(525)1(1057)1(3105231 222 222460 6 +−+++− ++−++−+−= JJJJJJ JJJJJJJJJO zzzz [ ]( ) ( ) [ ]38)1(1138)1(11 444414424146 −+−+++−+−= −+−+ JJJJJJJJJJO zz In the notation of Lea Leak and Wolf: 4 6 0 66 4 4 0 44 O21OO O5OO −= += 6. References AL ASSADI, K.F., MACKENZIE, I.S. and McCAUSLAND, M.A.H. 1984. Domain and domain-wall nuclear magnetic resonance of 265Ho in (Gd,Y)Fe2, J Phys F: Met. Phys 14: 525. BUNBURY, D.ST.P., CARBONI, C. and McCAUSLAND, M.A.H. 1989. Field dependence of the hyperfine splittings in holmium hydroxide J. Phys. Condens. Mat. 1: 1307. BUSHOW, K.H.J. 1977. Intermetallic compounds of Rare Earth and 3-d transition metals, Rep Prog Phys 40: 1179-1256. CARBONI, C., MACKENZIE, I.S. and McCAUSLAND, M.A.H. 1989. Nuclear Magnetic Resonance at microwave frequencies, Hyperfine Interactions 51: 1139. CARBONI, C., BUNBURY, D.ST.P. and McCAUSLAND, M.A.H. 2009. Crystal-field anisotropy of TmFe2. To be published. DE AZEVEDO, W.M., MACKENZIE, I.S. and BERTHIER, Y.J. 1985. NMR of 159TbFe2 and TbFe3. J Phys F. Met. Phys. 15: L243. DE AZEVEDO, W.M. 1986. The hyperfine interaction of terbium in magnetically ordered compounds, PhD thesis, University of Manchester, UK. DUMESNIL, K., DUFOUR, C., MANGIN, PH. and ROGALEV, A. 2002. Magnetic springs in exchange coupled DyFe2-YFe2 superlattices: an element selective X-ray magnetic circular dichroism study, Phys. Rev. B, 65: 094401-094406. DE LA FUENTE., C., ARNAUDAS, J.I., CIRIA, M., DEL MORAL, A., DUFOUR, C., MOUGIN, A. and DUMESNIL, K. 2001. Magnetic and magnetoelestic behavior of epitaxial TbFe2-YFe2 bilayers, Phys. Rev. B 63: 054417-054427. FUJIWARA, A., ASAKURA, K., HARADA, I., OGASAWARA, H. and KOTANI, A. 2005. XMCD study on electronic and magnetic states of rare earth 5d electrons in Laves compounds RFe2, Phys. Scr. T115: 113- 114. GERMANO, D.J. and BUTERA, R.A. 1981. Heat capacity of, and crystal-field effect in, the RFe2 intermetallic compounds, Phys. Rev. B, 24(7): 3912-3927 GORDEEV, S.N., BEAUJOUR, J-M L., BOWDEN, G.J., RAINFORD, B.D., DE GROOT, P.A.J., WARD, R.C.C., WELLS, M.R. and JANSEN, A.G.M. 2001. Giant magnetoresistence by exchange springs in DyFe2-YFe2 superlattices, Phys. Rev. Lett. 87: 186808-186812. ANISOTROPY IN RFe2 INTERMETALLICS 83 GRATZ, E., KOTTAR, A., LINDBAUM, A., MANTLER, M., LATROCHE, M., PAUL-BONCOUR, V., ACET, M., BARNER, C., HOLZAPFEL, W.B., PACHECO, V. and YVON, K. 1996. Temperature and pressure-induced structural transitions in rare-earth deficient RNi2 Laves phases, J. Phys.: Condens. Matter 8: 8351-8361. LEA, K.R., LEASK, J.M. and WOLF, W.P. 1962. The raising of angular momentum degeneracy of f-electron terms by cubic crystal field, Journal of the Physics and Chemistry of Solids, 23: 1381-1405. LEE, S.J., LANGE, R.J., CANFIELD, P.C., HARMON, B.N. and LYNCH, D.W. 2000. Optical and magneto- optical properties of RFe2 and RCo2, Phys. Rev. B, 61: 9669-9679. LI Y, CARBONI, C., ROSS, J.W., MCCAUSLAND, M.A.H. and BUNBURY, D.ST.P, 1996. Hyperfine fields at lanthanides nuclei in the intermetallic compounds RFe2, J Phys Condens. Matter 8: 865. MCMORROW, D.F., McCAUSLAND, M.A.H. and HAN, Z.P. 1989. Exchange and crystal field interaction of Ho3+ in GdAl2 : a single crystal NMR study, J Phys Condens. Matter 1: 10439-10458. McCAUSLAND, M.A.H. and MACKENZIE, I.S. 1979. Nuclear Magnetic resonance in Rare Earth metals, Adv. Phys. 28: 305. MOUGIN, A., DUFOUR, C., DUMESNIL, K. and MANGIN, P. 2000. Strain induced magnetic anisotropy in single crystal RFe2 (110) thin films, Phys. Rev. B, 62: 9517-9531. OSTER, J., WIEHL, L., ADRIAN, H. and HUTH, M. 2005. Magnetic and magnetoelastic properties of epitaxial (211)-oriented RFe2 thin films J. Mag. Mag. Mater., 292: 164-177. PARLEBAS, J.C., ASAKURA, K., FUJIWARA, A., HARAD, I. and KOTANI, A. 2006. X-ray magnetic circular dichroism at rare-earth L23 absorption edges in various compounds and alloys, Physics Reports, 431-1: 1-38. PAUL-BONCOUR, V. 2004. Comparison of the influence of hydrogen on the magnetic properties of RMn2 and RFe2 Laves phase compounds, Journal of Alloys and Compounds, 367: 185-190. SAWICKKI, M., BOWDEN, G.J., DE GROOT, P.A.J., RAINFORD, B.D., BEAUJOUR, J-M L., WARD, R.C.C. and WELLS, M.R. 2000. Exchange springs in antiferromagnetically coupled DyFe3-YFe2 superlattices, Phys Rev B 62: 5817-5820. STEVENS, K.W.H. 1952. Matrix Elements and Operator Equivalents Connected with the Magnetic Properties of Rare-Earth Ions, Proc. Phys. Soc. A 65: 209-215. Received 24 November 2008 Accepted 23 August 2009