SQU Journal for Science, 2015, 20(1), 62-69 © 2015 Sultan Qaboos University 62 Investigation of the Direction of Spontaneous Magnetisation at 0 K in HoCo2 and TmFe2 Carlo Carboni* and Eman Al-Abri Department of Physics, College of Science, Sultan Qaboos University, P.O. Box: 36, PC 123, Al-Khod, Muscat, Sultanate of Oman. *Email: carlo@squ.edu.om. ABSTRACT: The direction of spontaneous magnetisation at 0 K across a phase orientation boundary in rare-earth laves phase compounds is computed as a function of the ratio of magnetic to electrostatic interactions. The phase boundary and range of parameters used for the computation correspond to the portion of the orientational phase diagram where experimental data can be obtained. In the series of pseudo-binary compounds (Gdx Ho1-x )Co2 and (Yx Ho1-x )Co2 the range of parameters covers the continuous transition from the <100> to the <110> orientations. The nmr measurements on this series show that the direction of spontaneous magnetisation is at intermediate orientation between the <100> and <110> directions. The material TmFe2, which is in the vicinity of a discontinuous orientational boundary, is also discussed. Keywords: Lanthanides; Crystal field; Laves phase; Phase transition. TmFe2و HoCo2( لمركبي okعنذ الصفر المطلق )دراسة االتجاه للتمغنط التلقائي ايمان العبريو كارلو كاربوني )انصفش انمطهق( عبش حذ تُجًٕ انطُس انمغىاطٕسٓ فٓ مشكباث انعىاصش األسضٕت K 0تم حساب اتجاي انتمغىظ انتهقائٓ عىذ دسجت انحشاسة ملخص: باعتباسي دانت فٓ وسبت انتفاعالث انكٍشبٕت إنّ انتفاعالث انمغىاطٕسٕت. حذ انطُس انمغىاطٕسٓ َمذِ –( laves phasesانىادسة َانتٓ نٍا أطُاس الفز ) رنك انجزء مه بٕاوٓ انطُس انتُجٍٕٓ انزْ ٔمكه استخالص وتائج تجشٔبٕت مىً. أما بانىسبت نسهسهت انمعامالث انتٓ تم استخذامٍا فٓ انحساب تخص انّ اال 〈 〉 فإن مذِ انمعامالث ٔغطٓ اإلوتقال انمتصم مه اإلتجاي . ( ) ( مثم pseudo – binaryانمشكباث انتٓ ٔطهق عهٍٕا) َُْ انمغىاطٕسٓ نٍزي انسهسهت مه انمشكباث تُضح أن اتجاي انتمغىظ انتهقائٓ ٔكُن متُسطا بٕه االتجإٌه ان قٕاساث انشوٕه انى 〈 〉 تجاي 〈 〉 َ 〈 〉 َانتٓ تقع عهّ مقشبت مه حذ تُجًٕ انطُس انمغىاطٕسٓ انمىقطع. تىاقش ٌزة انُسقً أٔضا انمادة ، انحقم انبهُسْ ، الثٕىٕذاث. اوتقال طُسْ ، طُس الفز: كلمات مفتاحية 1. Introduction ne of the characteristics of rare-earth ions is that the spin-orbit coupling in the 4-f shell is much larger than the interactions experienced by the ion when embedded in a crystal. As a result J, the total angular momentum of the 4f shell, remains, in a first approximation, a good quantum number in the solid [1]. The states of the ion subjected to the crystal-field and the magnetic interactions experienced in the solid are well described to first order by the simple effective electronic Hamiltonian [2]. (1) where Hz stands for the magnetic interaction of the 4f shell with the host crystal and Hcf is the crystal-field interaction which is the electrostatic interaction of the 4f shell with the electric field gradients due to the surrounding ions. In the Zeeman representation, the eigenstates of the electronic Hamiltonian are expressed as linear combinations of J M states. Since at room temperatures and below only the ground J-manifold is populated the index J is dropped and the electronic states are written | ⟩ ∑ | ⟩ O INVESTIGATION OF THE DIRECTION OF SPONTANEOUS MAGNETISATION 63 Each electronic state | ⟩ has a magnetic moment associated to the expectation value ⟨ | | ⟩ of the z- component of the angular momentum J tied to the corresponding aspherical distribution of the 4f electrons. The 4f shell and therefore the magnetic moment have a preferred orientation with respect to the electric field gradients at the ion’s site in the crystal. Different states have different preferred orientations. This feature is responsible for the wealth of magnetic structures and magnetic transitions displayed by rare-earth compounds. In particular, when a low-lying excited state has a strongly preferred orientation conflicting with that of the ground state, a change in the orientation of spontaneous magnetization can occur as the temperature is increased and the excited state becomes populated [3]. We note in passing that the gadolinium ion Gd 3+ being in an S state, the crystal-field interaction is zero in first order and therefore, in first order, Gd 3+ has no preferred orientations. At zero-temperature, the direction of spontaneous magnetization in magnetic rare-earth compounds is determined solely by the ground state of the 4f shell. The nature of the ground state depends of course on the parameters in the Hamiltonian in equation 1; as the parameters are varied, the direction of spontaneous magnetisation may change. For ions at a site of cubic symmetry the Hamiltonian can be expressed in the dimensionless form given in equation 7 below which it depends on two parameters only. Al-Abry et al. [4, 5] have done a systematic investigation of the direction of spontaneous magnetization at zero temperature in the two-parameter space for all the values of J found in the ground state of lanthanide ions. In the present paper we discuss the direction of spontaneous magnetization in two compounds, HoCo2 and TmFe2. In these compounds the rare-earth ion is at a site of cubic symmetry and the parameters for these compounds place them near a boundary between two regions in the parameter space with different orientations of spontaneous magnetization for the ground state. The crystal-field Hamiltonian is expressed as a linear combination of the Steven’s operators [6]. ∑ ∑ The are the crystal-field parameters. At a site of cubic symmetry there are only four non-vanishing terms in the expansion: and are the non-zero parameters. Furthermore there is a relation between the two fourth order and between the two sixth order terms; therefore only two parameters are required to describe the crystal field [7]. The crystal-field Hamiltonian is written (2) with and Following Lea, Leask and Wolf [6] the Hamiltonian is recast in the form { ( ) } (3) where W is an energy scale factor and the unitless parameter x gives the ratio of the fourth to sixth order crystal-field terms. ( | |) F4 and F6 are numerical factors depending on the particular value of J [7]; the magnetic interaction is (4) where BR is a vector that represents the sum of all the effective magnetic fields acting on the rare-earth ion; these can arise because of exchange, the dipolar field or an externally applied field. gJ is the Landé g-factor and B is the Bohr magneton. Following McMorrow et al. [8] we introduce the vector   )5( CARLO CARBONI and IMAN AL-ABRI 64 and write the electronic Hamiltonian { ( ) } )6( This Hamiltonian involves three parameters, W, x and the vector  For the purpose of the present work which is concerned with the direction of spontaneous magnetization, it is more convenient to recast the Hamiltonian in a two- parameter dimensionless form by dividing by the energy factor W. The Hamiltonian becomes { ( ) } )7( The magnitude of the vector parameter represents the relative strength of the magnetic to the electrostatic interactions; the balance between the relative strengths of |R| and x determines the direction of spontaneous magnetization [4, 5]. Given the values of x and |R| and the orientation (of R, the direction of spontaneous magnetization which is the direction of the angular momentum J can be in principle computed by diagonalising the Hamiltonian 7. However, for arbitrary values of (the computed direction of J will not in general be parallel to the input direction of R. To be the direction of spontaneous magnetization, the value of () must satisfy the two conditions that the computed direction of J is parallel to R, and that the energy of the state is at a minimum with respect to () [4, 5]. 2. Results and discussion In this work the evolution of the direction of spontaneous magnetization in the parameter space was computed as in reference [4] by varying systematically the parameters x and |R| and finding for each value of (x, |R|) which values of the polar angles () give the lowest energy for the ground state and the computed orientation of J parallel to the input  Table 1 gives a summary of the values of W, x and R for HoCo2 and TmFe2 derived from data published by various authors. The different experimental techniques used by the authors lead to a scatter in the values of the parameters. In most cases such an uncertainty in the parameters is not a problem. However, in the cases of Ho 3+ in HoCo2 and Tm 3+ in TmFe2, because these compounds are close to a boundary between regions of different directions of spontaneous magnetization in the {R, x} parameter space, a small variation in the values of the parameters can lead to an erroneous direction of spontaneous magnetization. Table 1. The parameters for HoCo2 and TmFe2 derived from published data. a: Gignou et al.[9], b: Aubert et al. [10], c: Castets and Gignoux [11], d: Germano and Butera [12], e: Martin et al.[13] derived values of B4 and B6 from the data of Atzimony et al.[14], f: Yanovsky et al. [15], g and h derived by Bleaney [18] from Rhyne and Koon [16] and Koon and Rhyne [17], i: Bleaney [18]. Ion J Compound W (K) x |R| ref Ho 3+ 8 HoCo2 0.60 -0.4687 46.6 a 0.60 -0.4665 45.8 b 0.30 -0.4839 60.2 c Tm 3+ 6 TmFe2 0.59 0.555 99 d 0.58 0.649 89 e 0.54 0.721 86 f 0.59 0.353 78 g 0.61 0.627 76 h 0.44 0.64 105 i  The relevant portions of the orientational phase diagram for the two ions Ho 3+ (J = 8) and Tm 3+ (J = 6) are given in Figures 1-A and 1-B. The published parameters for the two compounds given in table 1 are also shown in the figures. In the notation of Al-Abry et al (4, 5), the boundary in the case of J = 6 is of the type Da where the transition is discontinuous across the boundary. At a critical value of the parameters the two orientations (110) and (111) have the same energy but all the orientations in the plane containing these two directions have a higher energy. There is a sudden change of direction of spontaneous magnetization at the critical value. For J = 8 the boundary is of the type C where there are two critical values for the parameters; between these two values the direction of spontaneous magnetization is at an intermediate orientation between the (100) and the (110) axis. There is a continuous rotation of the direction of spontaneous magnetization. Figures 1-C and 1-D show the computed directions of spontaneous magnetization across the boundary for the two cases. INVESTIGATION OF THE DIRECTION OF SPONTANEOUS MAGNETISATION 65 Figure 2. The variation of across the boundary. A: for J = 8 and B for J = 6. The computed values of 〈 〉 corresponding to figures 1-C and 1-D are shown in figure 2-A and 2-B. As one might expect, in figure 2-A for the case of Ho 3+ , the value of 〈 〉 varies continuously across the boundary and the first derivative is continuous. For Tm 3+ , 〈 〉 is continuous through the boundary but there is a discontinuity with a change of sign in the first derivative. We note that for J = 8 the change in the value of 〈 〉 across the boundary is about 3%,whereas in the case of J = 6 the variations of 〈 〉 in the vicinity of the boundary are of the order of 0.02%. Figure 1. A and B: the relevant areas of the orientational phase diagrams for J = 6 and J = 8; the letters in brackets refer to the references in table 1. C and D are the computed direction of spontaneous magnetisation across the phase boundary. C D A 40 50 60 70 80 7.7 7.6 7.5 7.4 │R│ < J z > J = 8 x = -0.466 <110> <100> B J = 6 │R│= 80 5.926 5.922 5.920 5.918 5.916 5.914 5.912 5.924 0.6 0.62 0.64 0.66 0.68 <111> <110> < J z > x CARLO CARBONI and IMAN AL-ABRI 66 We consider first the case of HoCo2. The direction of spontaneous magnetization was investigated experimentally by 165 Ho:Ho 3+ NMR. The holmium nucleus has spin 7/2, the NMR spectrum presents seven quadrupole-split lines. The frequency at of the central line of the spectrum is directly related to the expectation value 〈 〉 which can be computed from the Hamiltonian in equation 7 [1, 2] 〈 〉 (8) Where a0 is the free-ion dipolar hyperfine coupling constant (a0 = 6497 MHz for holmium); a” is the contribution arising from the direct interaction of the nucleus with the surroundings in the crystal and a” is a small contribution that can be treated as a free parameter. Because of the large value of a0, any small change in the value of 〈 〉 results in a large shift in the NMR frequencies. In materials where the value of 〈 〉 depends strongly on the orientation of the ion with respect to the crystal axis, the NMR spectrum is a powerful tool to determine the direction of magnetization [19]. In practice it is difficult to control experimentally the parameters |R| and x. However, in the case of HoCo2, the value of |R| can be decreased by doping the compound with non-magnetic yttrium or increased by doping with gadolinium which has a larger spin than holmium and, as mentioned in the introduction, does not contribute to the anisotropy because Gd 3+ is in an S state. A series of compounds (Ho1-yGdy)Co2 (y = 0.0, 0.05, 0.1, 0.2, 0.3, 0.4, 0.97) was prepared by arc-melting. Assuming for the sake of computation that the value of |R| is 46.6 for pure HoCo2 (table 1) the concentrations of gadolinium chosen correspond to |R| = 46.6, 48.1, 49.8, 54.0, 58.1, 62.0, 85.4 respectively. To obtain lower values of |R| a compound doped with 5% of yttrium corresponding to the value of |R| = 44.3 was prepared. The value of the molecular field experienced by the Ho 3+ ion in the pseudo-binary compounds was estimated following McMorrow et al [8]. It is expected that these values of |R| take the compound across the orientational phase boundary. The NMR spectra were taken on the 2-8 GHz spectrometer at the University of Manchester (UK) [20] at liquid helium temperatures. A preliminary outcome of some of the still unpublished measurements on some of the (Ho,Gd)Co2 compounds was presented at the XIII Congress Ampere [21]. In the concentrated HoCo2 compound and in the (Ho,Gd)Co2 compounds with low gadolinium concentration two overlapping spectra are observed. For concentrations of gadolinium above 50% a single spectrum is observed which is consistent with the <100> orientation of spontaneous magnetisation. The NMR spectrum of 165 Ho:HoCo2 is shown in figure 3. The presence of two spectra indicates that there are two inequivalent sites with slightly different strengths of the exchange interaction. In the transition region this results in a significant shift of the NMR spectrum. The explanation for the two sites has not yet been unequivocally established. The two obvious possible causes, an impurity in the starting materials for the sample preparation or closure domains or domain walls were discarded by repeating the measurements with specimens of HoCo2 prepared at other institutions including a single crystal on loan from Laboratiore Louis Neel, Grenoble. An 8 T field was applied to the single crystal specimen. In all experimental conditions two spectra are observed in HoCo2. We will not attempt here to give any plausible explanations for the presence of the two sites. Figure 3. The observed spectrum for HoCo2 is the superposition of two 7-lines spectra. For each compound we have derived the values of at of each of the two overlapping spectra. In figure 4 we have plotted all the measured values of at as a function of concentration. The conclusion than can be drawn from the figure 5 is that two spectra are observed only for compounds with a value of │R│ that places them in the transition region. Therefore, the pure compound HoCo2 is inside the transition region. HoCo2 polycrystal T = 4.2 K Ba = 0 6000 6200 6400 6600 6800 7000 Frequency (MHz) E c h o a m p li tu d e ( a r b . u n it s) INVESTIGATION OF THE DIRECTION OF SPONTANEOUS MAGNETISATION 67 Figure 4. Summary of the frequency of the central line of all the spectra observed. The dashed line is a guide to the eye; the vertical line corresponds to the HoCo2 compound. In order to obtain a comparison between computation and experiment we have derived the values of using equation 8. We have considered a” as a free parameter for each compound and allowed a” to vary linearly with concentration from 392 MHz in the concentrated compound to 356 MHz in the compound with 40% gadolinium. We have also assumed that the low frequency spectrum corresponds to the values of │R│ estimated for the specimen preparation and that the high frequency spectrum corresponds to a site with a value of │R│ increased by 9. The values of thus obtained are compared to the values computed using the parameters of Aubert et al (10). The figure shows a very good agreement between the measurements and the computation. Figure 5. The values of derived from the measurements are compared to the computed values (solid line). Turning now to TmFe2 we note that some of the parameters given in table 1 give the <111> as direction of spontaneous magnetisation for the ground state, while others give the <110>. The parameters of Germano and Butera (12) give a reorientation of the spontaneous magnetisation as a function of temperature from the <110> at low temperatures to <111> above 55 K. As noted in Figure 3, there is only a small (0.02%) change in across the phase boundary. Therefore, nmr cannot be used to ascertain the direction of spontaneous magnetisation at low temperatures. The direction of magnetisation can be obtained from the 57 Fe Mössbauer spectrum since the iron sites are not magnetically equivalent [14]. A detailed study of the 57 Fe Mössbauer spectrum by Al-Sariri (22) confirmed that the direction of magnetisation in TmFe2 is <110> at all temperatures. The specimen used for the Mössbauer measurements by Al-Sariri was a crushed single crystal obtained from University of Birmingham (UK). 20 40 60 80 10 F r e q u e n c y ( M H z ) 6700 6600 6500 6400 6300 % of Gd % of Y 40 50 60 70 80 7.7 7.6 7.5 7.4 < J Z > │R│ CARLO CARBONI and IMAN AL-ABRI 68 3. Conclusion This work has shown that a boundary in the orientational phase diagram can place constraints in the range of the values for the parameters of the Hamiltonian. These constraints make it possible to refine the values of the parameters. The nmr data shows that HoCo2 is just inside the intermediate region on the <110> side of the boundary. Therefore, the parameters that describe best the Ho 3+ ion in this compound are those of Abert et al [4]. For TmFe2 most of the experimental data available would indicate that the <110> is the orientation of spontaneous magnetization at all temperatures; therefore the parameters given by Yanovsky et al. [15] are probably incorrect. However, the analysis given by Sariri, based on unpublished high field 169 Tm:TmFe2 nmr [23] and in some pseudo-binary alloys would indicate that magnetostriction is not negligible. An extra term needs to be included in the crystal-field Hamiltonian in order to take the magnetostriction distortion into account. 4. Acknowledgements The experimental part of this work was done during the time the author C C was at the University of Manchester (UK) supported by grants from the UK Engineering and Physical Science Research council (EPSRC). The authors are also grateful to Prof MAH McCausland and Dr D StP Bunbury for comments, discussions and support. References 1. McCausland, M.A.H. and Makenzie, I.S. Nuclear Magnetic Resonance in Rare earth Metals, Taylor and Francis Ltd., London, 1980. 2. Bunbury, D.St.P., Carboni, C., McCausland, M.A.H., Cone, R. and Mroczkowsky, S. and Ross, J.W. The hyperfine splitting of the ground state and first excited state of holmium in yttrium hydroxide. J. Phys. Condens. Matter. 1995,7, 6125-6136. 3. Koehlera, M.R., Garleab, V.O., McGuirec, M.A., Jiaa, L., Keppensa, V. Spin reorientation and magnetoelastic coupling in Tb6Fe1−xCoxBi2 (x = 0, 0.125, 0.25, and 0.375) alloy system. J. Alloys and Compounds, 2014, 615, 514– 520. 4. Al Abry, E. 2003 M.Sc. Thesis, The orientation of the direction of spontaneous magnetization of rare-earth ions at a site of cubic symmetry. Department of Physics, College of Science, Sultan Qaboos University, Sultanate of Oman. 5. Carboni, C., Al-Abry, E. and Arafin, S. The direction of spontaneous magnetisation of lanthanide ions at a site of cubic symmetry. SQU J. Science, 2007, 12(1), 53. 6. Stevens, K.W.H. Matrix elements and operator equivalents connected with the magnetic properties of rare earth ions, Proc. Physi. Soci., A, 1952, 65, 209. 7. Lea, K.R., Leask, M.J.M. and Wolf, W.P. The raising of angular momentum degeneracy of f-electron terms by cubic crystal fields. J. Phys. Chem. Solids, 1962, 23, 1381-1405. 8. McMorrow, D.F., McCausland, M.A.H., Han, Z.P. Exchange and crystal field interaction of Ho 3+ in a GdAl2 single crystal: an NMR study. J. Phys Condens. Matter, 1989, 1, 10439-10458. 9. Gignoux, D., Givord, F., and Lemaire, R. Magnetic properties of single crystals of GdCo2, HoNi2 and HoCo2 Phys. Rev B, 1975, 12, 3878-3884. 10. Aubert, G., Gignoux, D., Givord, F., Lemaire, R. and Michelutti, B. Magnetization reorientation in HoCo2, Solid State Commun, 1978, 25, 2, 85–87. 11. A Castets, D Gignoux and B Hennion, New aspects of the magnetic excitations in HoCo2. Phys. Rev. B, 1982, 25, 337-348. 12. Germano, D.J. and Butera, R.A. Heat capacity of, and crystal field effect in, the RFe2 intermetallic compounds, Phys. Rev. B, 1981, 24(7), 3912-3927. 13. Martin, K.N., de Groot, P.A.J., Rainford, B.D., Wang, K., Bowden, G.J., Zimmermann, J.P. and Fanghor, H. Magnetic anisotropy in the cubic Laves REFe2 intermetallic compounds. J. Phys. Condens. Matter, 2006, 18, 459- 478. 14. Atzmony, U., Dariel, M.P., Bauminger, E.R., Lebenbaum, D., Nowik, I. and Ofer, S. Spin-orientation diagrams and magnetic anisotropy of rare-earth-iron ternary cubic laves compounds. Phys. Rev., 1973, B 7, 4220-4232. 15. Yanovsky, R., Bauminger, E.R., Levron, D., Nowik, I. and Ofer, S. Crystalline fields, exchange interactions and spin reorientations in TmxHo1−x Fe2 systems, studied by a 170 Yb Mössbauer probe. Solid State Commun. 1975, 17, 1511-1514. 16. Rhyne, J. and Koon, N.C. Magnetic excitations in HoFe2. J. Appl. Phys., 1978, 49, 2133-2135. 17. Koon, N.C. and Rhyne, J. Excited state spin waves in ErFe2. Solid State Commun., 1978, 26, 537-540. 18. Bleaney, B., Bowden, G.J., Cadogan, J.M., Day, R.K. and Dunlop, J.B. A Mossbauer study of the cubic Laves phase intermetallic compound TmFe2. J. Phys. F, Met. Phys., 1982, 12, 795-811. 19. Waind, P.R., Mackenzie, I.S. and McCausland, M.A.H. Anisotropy and asymmetry of the hyperfine spectrum of holmium in magnetically ordered RAl2 alloys. J. Phys F, 1983, 13, 1041-1056. 20. Carboni, C., Mackenzie, I.S. and McCausland, M.A.H. Nuclear Magnetic Resonnance at microwave frequencies, Hyperfine Interactions, 1989, 51, 1139-1144. INVESTIGATION OF THE DIRECTION OF SPONTANEOUS MAGNETISATION 69 21. Bunbury, D.StP., Carboni, C., McCausland, M.A.H. and Prakash, O. NMR studies of spin reorientation in (Ho,Gd)Co2. Unpublished Proc. XXIII Congress Ampere (Istituto Superiore di Sanita Roma), 1986. 22. Al Sariri, M. M.Sc. Thesis, Crystal Field and Magenic Anisotropy in TmFe2. Department of Physics, College of Science, Sultan Qaboos University, Sultanate of Oman, 2010. 23. Carboni, C. Anisotropy in RFe2 intermetallics: an nmr study. SQU J. Science, 2009, 14(1), 71-83. Received 14 October 2014 Accepted 2 January 2015