Complex exchange mechanism driven ferromagnetism in half-metallic Heusler Co$_{2}$TiGe: Evidence from critical behavior
Shubhankar Roy, Nazir Khan, Ratnadwip Singha, Arnab Pariari, Prabhat, Mandal

TL;DR
This study investigates the critical magnetic behavior of half-metallic Co$_{2}$TiGe, revealing complex exchange interactions and critical exponents indicating long-range ferromagnetic order with significant spin fluctuations.
Contribution
It provides the first detailed critical exponent analysis of Co$_{2}$TiGe, demonstrating the complex exchange mechanisms and deviation from classical models in this half-metallic Heusler alloy.
Findings
Critical temperature $T_{C}$=371.5 K identified.
Critical exponents suggest mean-field-like behavior with long-range interactions.
Evidence of sizeable critical spin fluctuations and magnetic inhomogeneity.
Abstract
We have investigated the critical phenomenon associated with the magnetic phase transition in the half-metallic full-Heusler CoTiGe. The compound undergoes a continuous ferromagnetic to paramagnetic phase transition at the Curie temperature =371.5 K. The analysis of magnetization isotherms in the vicinity of , following modified Arrott plot method, Kouvel-Fisher technique, and critical isotherm plot, yields the asymptotic critical exponents =0.495, =1.324, and =3.67. The self-consistency and reliability of the obtained exponents are further verified by the Widom scaling relation and scaling equation of states. The mean-field-like value of the critical exponent suggests long-range nature of the exchange interactions, whereas the values of the critical exponents and , imply sizeable critical spin fluctuations. The…
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8| Material/Model | Ref. | Technique | ||||
|---|---|---|---|---|---|---|
| Co2TiGe | This work | MAP | 371.5 | 0.495 | 1.325 | 3.677 |
| This work | KF | 371.5(1) | 0.495(2) | 1.324(4) | 3.675 | |
| This work | CI | 3.671(1) | ||||
| Ni | kaul | KF | 0.391(10) | 1.314(16) | 4.39(2) | |
| Gd80Au20 | sjpoon | KF | 0.44(2) | 1.29(5) | 3.96(3) | |
| Mean Field | Eugene | Theory | 0.5 | 1 | 3 | |
| 3D Heisenberg | Banerjee ; Campostrini | Theory | 0.365 | 1.386 | 4.8 | |
| 3D XY | Banerjee ; Guillou | Theory | 0.345 | 1.316 | 4.81 | |
| 3D Ising | Banerjee ; Guillou | Theory | 0.325 | 1.24 | 4.82 | |
| LR exchange: | ||||||
| = 3, =3, = 1.88 | MEFisher | RG- | 0.393 | 1.32 | 4.36 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Complex exchange mechanism driven ferromagnetism in half-metallic Heusler Co2TiGe: Evidence from critical behavior
Shubhankar Roy
Saha Institute of Nuclear Physics, HBNI, 1/AF Bidhannagar, Kolkata 700 064, India
Vidyasagar Evening College, 39, Sankar Ghosh Lane, Kolkata 700 006, India
Nazir Khan
Saha Institute of Nuclear Physics, HBNI, 1/AF Bidhannagar, Kolkata 700 064, India
Ratnadwip Singha
Saha Institute of Nuclear Physics, HBNI, 1/AF Bidhannagar, Kolkata 700 064, India
Arnab Pariari
Saha Institute of Nuclear Physics, HBNI, 1/AF Bidhannagar, Kolkata 700 064, India
Prabhat Mandal
Saha Institute of Nuclear Physics, HBNI, 1/AF Bidhannagar, Kolkata 700 064, India
Abstract
We have investigated the critical phenomenon associated with the magnetic phase transition in the half-metallic full-Heusler Co2TiGe. The compound undergoes a continuous ferromagnetic to paramagnetic phase transition at the Curie temperature =371.5 K. The analysis of magnetization isotherms in the vicinity of , following modified Arrott plot method, Kouvel-Fisher technique, and critical isotherm plot, yields the asymptotic critical exponents =0.495, =1.324, and =3.67. The self-consistency and reliability of the obtained exponents are further verified by the Widom scaling relation and scaling equation of states. The mean-field-like value of the critical exponent suggests long-range nature of the exchange interactions, whereas the values of the critical exponents and , imply sizeable critical spin fluctuations. The half-metallic itinerant character of Co2TiGe in the presence of magnetic inhomogeneity may result in such a strong deviation from the three-dimensional Heisenberg values (=0.369, =1.38 and =4.8) of the critical exponents towards the mean field values (=0.5, =1 and =3). The results suggest complex nature of exchange couplings that stabilize the long-range ferromagnetic ordering in the system and are consistent with the earlier theoretical studies on the exchange mechanism in Co2TiGe.
I Introduction
The discovery of ferromagnetism by Friedrich Heusler in Cu2MnAl, which does not contain any of the pure ferromagnetic elements, initiated the evergrowing research interests in the inter-metallic compounds Heusler . Subsequently, numerous experimental works Webster1 ; Webster2 ; Hamzic have been carried out and several theoretical models Price ; Kubler ; JLMor ; Kurtulus ; Sasioglu ; Thoene ; Trudel have been proposed to understand the microscopic interactions that lead to the long-range magnetic ordering in this class of materials. The fascinating half-metallic character Groot ; Galanakis makes them very promising candidates for the spintronics applications at room temperature Chappert . Moreover, recently, novel topological semi-metal states Chang have been discovered in these systems. The Co2-based ternary intermetallic Heusler compounds viz., Co2M^{\prime}$$Z, where is a transition metal and is a main group element, received special interest and became the most widely studied systems because of their very high Curie temperatures Trudel .
Co2TiGe crystallizes in cubic 21 structure (space group: -), which consists of four inter-penetrating face-centered-cubic (fcc) lattices. The crystallographic positions of Co atoms are (0,0,0) and (1/2, 1/2, 1/2) and that of the Ti and Ge atoms are (1/4, 1/4, 1/4) and (3/4, 3/4, 3/4), respectively. There are several theoretical and experimental reports exploring the structural, electronic, and magnetic properties of Co2TiGe Barth . The spin-resolved band structure calculations show that the majority spin-band has a metallic character, whereas the minority spin-band exhibits semiconducting behavior with a band gap of about 0.5 eV around the Fermi level Chang ; Barth . Further, the inclusion of the spin-orbit coupling (SOC), albeit very weak, in the band structure calculation shows novel topological Weyl semimetal state Chang . The temperature dependence of the resistivity measurement in Co2TiGe shows typical metallic behavior down to 2 K Barth . The origin of ferromagnetism in half-metallic Co2-based Heusler alloys is a rather complicated issue and remains one of the most interesting problems in modern magnetism. They exhibit Slater-Pauling-type behavior for the magnetization, where the saturation magnetization scales with the total number of valance electrons in the unit cell, leading to a characteristic integer magnetic moments. Magnetic Heusler alloys are traditionally believed to be ideal local-moment systems Telling and their exchange couplings can be described by a Heisenberg Hamiltonian. However, several theoretical reports claim the presence of complex exchange interactions of the localized magnetic moments in the Co2-based Heusler compounds Kubler ; Kurtulus ; Sasioglu ; Thoene ; Trudel . Both the short-range direct exchange coupling between the nearest neighbor spins and the coupling between further neighbor spins mediated by the Ruderman-Kittel-Kasuya-Yosida (RKKY) type long-range interaction, are believed to be responsible for the long-range magnetic ordering Trudel . It is concluded that inter-sublattice and intra-sublattice exchange interactions of different strength stabilize robust ferromagnetic ordering with remarkably high Curie temperatures in the Co2-based Heusler compounds. However, a recent theoretical study shows that in the Co2Ti system, the magnetic moment of the Co atoms and that of the Ti atoms are not completely localized, which is further supported by the disorder local-moment calculation Barth . Thus, the presence of localized magnetic moment is not prerequisite for the occurrence of half-metallic ferromagnetism. Furthermore, significant orbital magnetic moments in different Co2M^{\prime}$$Z compounds as probed by x-ray magnetic circular dichroism studies Telling ; Klaer ; Miyamoto , indicate that the SOC may result in magneto-crystalline anisotropy. So, anisotropic exchange interaction may also be involved in the magnetic ordering phenomenon in these systems. The values of the critical exponents associated with a continuous magnetic phase transition, reflect the nature of the exchange mechanism, spin symmetry or the magnetic anisotropy, and the effective dimensionality of spin-spin interactions in the system Fisher . Therefore, by experimentally investigating the critical phenomenon associated with the ferromagnetic to paramagnetic phase (FM-to-PM) transition in Co2TiGe, we elucidate the nature of exchange mechanism in a Co2-based Heusler alloy.
II Experimental details
Co2TiGe compound was prepared by arc melting of stoichiometric amounts of its constituents in a highly purified argon atmosphere. The obtained ingots were then sealed in an evacuated quartz tube and annealed for about 3 weeks. Phase purity and the structural analysis of the sample were done using the powder x-ray diffraction (XRD) technique with Cu-Kα radiation in a Rigaku x-ray diffractometer (TTRAX III). The Rietveld profile refinement of the XRD pattern (Fig. 1) shows that the compound is single phase in nature with space group symmetry 21. The refined lattice parameters, ===5.818(2) Å, are in good agreement with that reported previously Barth ; Barth2 . The magnetization measurements were done using a Superconducting Quantum Interference Device-Vibrating Sample Magnetometer (SQUID-VSM) (MPMS 3, Quantum Design) in fields up to 5 T. The sample used for the magnetic measurements is of approximate dimensions mm3. To minimize the demagnetization effect the external magnetic field was applied along the longest sample direction. The magnetization data were taken over the temperature range from 366 to 377 K at 1.0 K interval. To achieve good thermal equilibrium, we have stabilized each temperature for 45 minutes. For each M(H) isotherm, the magnetic field was increased from 0 to 5T and then reduced to 0. We didn’t find any difference in between the increasing and decreasing field.
III Results and Discussions
Figure 2 shows the magnetic field dependence of dc-magnetization () at 2 K. In the inset, we have plotted as a function of temperature, measured at 500 Oe. The observed behavior is quite similar to the earlier reports Barth ; Barth2 . As mentioned before, the Co-based Heusler alloys, which are half-metallic ferromagnets (FMs), exhibit the Slater-Pauling-type behavior of the magnetization as given by /f.u. Here, is the total magnetic moment and is the total number of valence electrons in the unit cell of the compound. For Co2TiGe, the value of is 26. Therefore, according to the above relation, the total magnetic moment should be 2/f.u. From Fig. 2, the saturation magnetization is estimated to be /f.u., which implies that the compound obeys the Slater-Pauling rule. As shown in the inset of Fig. 2, the magnetization curve mimics a continuous or second order FM-to-PM phase transition. From the temperature derivative of this magnetization curve, the Curie temperature () is estimated to be 371 K.
A continuous magnetic phase transition exhibits critical behavior of different thermodynamic variables, governed by the critical fluctuation. The critical behavior is characterized by a set of static critical exponents. The values of these exponents depend only on the symmetry of the order parameter and the lattice dimensionality of the system, resulting in different universality classes of the continuous phase transition in uniform FMs. In the vicinity of , the spin-spin correlation length () diverges as, , which leads to universal scaling laws for the spontaneous magnetization, (0, ), and inverse initial susceptibility, (0, ). below , above , and the magnetization isotherm at are characterized by a set of static critical exponents , , and , respectively. These exponents are defined as follows Eugene :
[TABLE]
[TABLE]
[TABLE]
where is the reduced temperature and , , and are the critical amplitudes. The scaling theory predicts that in the close vicinity of the phase transition, the magnetic equation of state for the system can be expressed as
[TABLE]
where and are regular functions for temperature above and below , respectively. In terms of the renormalized magnetization, m$$\equiv$$\mid\varepsilon\mid-β, and the renormalized field, h$$\equiv$$\mid$$\varepsilon$$\mid-(γ+β), the above scaling equation can be rewritten as
[TABLE]
The above equation implies that for the right choice of , , and , the vs. isotherms will fall onto two separate branches of the scaling function: for isotherms above and for isotherms below . This is also an important criterion to check whether the set of the critical exponents are the same below and above .
To characterize the FM-to-PM phase transition in the present system, a series of magnetization isotherms, as shown in Fig. 3(a), has been measured in the vicinity of the approximate Curie temperature, which was determined from the inset of Fig. 2. To test whether the framework of the Landau mean-field theory of magnetic phase transition holds for Co2TiGe, we have constructed the conventional Arrott plot Arrott of vs. isotherms as shown in Fig. 3(b). For the mean-field values of the critical exponents ( and ), Arrott plot should generate a set of parallel straight lines and the one that passes through the origin of the plot, corresponds to the isotherm at the exact critical temperature. However, the Arrott plot for Co2TiGe does not constitute a set of parallel straight lines and shows significant non-linearity with downward curvature even at high magnetic fields. Therefore, the Landau mean-field theory fails to describe the phase transition in Co2TiGe and suggests the presence of significant critical fluctuations. According to the Banerjee’s criterion Banerjee , the positive values of the slope of the curves in the Arrott plot, indicate second order nature of the phase transition. In Fig. 3(b), all the vs. isotherms show concave downward curvatures, which yield positive values of the slopes and, thereby, confirm second-order FM-to-PM phase transition in Co2TiGe.
To determine the true critical exponents, we have reanalyzed the magnetization isotherms using the modified Arrott plot Arrott2 method, which is based on the following Arrott-Noaks equation of state:
[TABLE]
Here, and are constants. The critical exponents predicted for three-dimensional (3D) system by theoretical models such as 3D Heisenberg (, ), 3D XY (, ), and 3D Ising model (, ) have been used to construct modified Arrott plots as shown in Figs. 4(a), 4(b), and 4(c), respectively. We can see from Figs. 4(a)-4(c) that at the high-field region the constructions show quasi-straight lines, which appear to be parallel to each other. To determine which one among these three models generates the best modified Arrott plot with parallel isotherms, we have compared the slope of these straight lines (obtained from the linear fits to the isotherms at high-field region) to that obtained for the isotherm at , for each modified Arrott plot. Here, the slope is defined as . The obtained normalized slopes, , for each model have been plotted as a function of temperature in Fig. 5. It is evident that the values of the deviate progressively from the ideal value of 1 both below and above . However, among these three models, the deviation is minimum for the 3D Heisenberg model. Therefore, the modified Arrott plot corresponding to 3D Heisenberg model [Fig. 4(a)] has been used for further refinement of the exponents and using a rigorous iterative method as described below.
In Fig. 4(a), the intercepts of the linear fits (to the isotherms in high-field region) with the axes and ()1/γ, give reliable values of and , respectively. The power law fits to the data following Eqs. (1) and (2), yield new values of the exponents and , which are used to construct a new modified Arrott plot. During the fitting, was varied to get the best fit results. From this plot, new and ()1/γ are obtained by linear fits to the isotherms at high-field region. So we get a new set of and . This procedure is repeated to achieve self-consistency, i.e., when the new values of the critical exponents become almost equal to the values before that iteration.
The final modified Arrott plot, constructed using =0.495 and =1.325, is shown in Fig. 4(d). The normalized slopes of the linear fits to the isotherms of this modified Arrott plot, are also shown in Fig. 5. Here, one can see that both above and below , is close to 1, implying that all the isotherms form parallel straight lines with the isotherm at =371.5 K, which almost passes through the origin of the plot. Thus, the more accurate value of is estimated to be 371.5 K. In the low-field region, the small curvatures in the isotherms of Fig. 4(d) is due to the reorientation of the domains at low magnetic fields. The temperature dependence of and , which are obtained from Fig. 4(d), are shown in Fig. 6(a). The power law fit following Eq. (1) to the data yields =0.494(3) and =371.5(1) and fitting of Eq. (2) to the data gives =1.33(3) and =371.5(1). Theses values of the exponents and , are almost same to that obtained from the modified Arrott plot [Fig. 4(d)]. Therefore, the estimated values of , , and , are very reliable and intrinsic.
The accurate values of the exponents as well as can also be obtained using the Kouvel-Fisher technique Kouvel . According to this method, vs. and vs. should be straight lines with slope and , respectively, and their intercepts on the -axis give the value of . The linear fits to these plots [shown in Fig. 6(b)] yield =0.495(2), =371.48(4) and =1.324(4), =371.51(2).
The critical exponent can be obtained from the critical isotherms i.e., the isotherm at as well as from the Widom scaling relation widom :
[TABLE]
The plot at =371.5 K is shown in Fig. 7 and the inset shows the same plot on the log-log scale. Following Eq. (3), from the inverse of slope of the linear fit to the log() vs. log() curve, we obtain =3.671(1). Using the values of and obtained from Figs. 6(a) and 6(b), the Widom scaling relation, Eq. (7), gives and , respectively. These values of are very close to that obtained from the critical isotherms (Fig. 7). Hence the obtained values of critical exponents and are self-consistent and obey the Widom scaling relation very well.
The values of the critical exponents , , and , obtained from various methods, are presented in the Table I along with the theoretically predicted values for different models. From the Table I, at a glance we see that the values of the exponents for the present system do not belong to any of the conventional universality class. Therefore, we must check whether these exponents can generate the scaling equation of state as given by Eq. (5). In Fig. 8(a), scaled- vs. scaled- has been plotted using the values of critical exponents , , and , calculated from the Kouvel-Fisher plot. The inset shows the same plot on log-log scale to expand the low-field region. From Fig. 8(a), it is clear that all the isotherms collapse onto two separate branches: one above and another below . The reliability of the values of the exponents and are further verified by plotting vs. , as shown in Fig. 8(b). In this case, all the isotherms also fall on two independent branches of the scaling function. Furthermore, in the vicinity of the transition, the scaling equation of state can be rewritten in another form,
[TABLE]
where is a regular scaling function. The scaling Eq. (8) implies that all the isotherms now collapse into a single curve. Following Eq. (8), we have plotted vs. in the inset of Fig. 8(b). The inset shows that all the experimental data fall on a single curve passing through , which is located at zero on the horizontal axis of the plot. Therefore, the scaling behavior shown in Figs. 8(a) and 8(b), implies that all the interactions in the vicinity of get properly renormalized and the set of the critical exponents are same below and above the for the present compound.
The static critical exponents for Co2TiGe, obtained from aforementioned rigorous methods, obey the scaling relation and the scaling equation of state. Therefore, these exponents seem to truly characterize the critical behavior in the compound. However, a comparison of these exponents with that predicted by different theoretical models for 3D system, confirms that estimated critical exponents for Co2TiGe do not belong to any conventional universality class and fall between the values predicted by the 3D Heisenberg model and the mean-field theory. It should be mentioned that the asymptotic critical exponents, which are temperature independent in the asymptotic critical regime, truly characterize the critical behavior in a material belonging to a universality class. Usually, the asymptotic critical regime is \mid$$\varepsilon$$\mid 10*-2* for homogeneous magnets and slightly larger for disordered and amorphous systems kaul . The reduced temperature range, used to derive the critical exponents for Co2TiGe, is 1.3101.4 10*-2*. Therefore, we have performed the above analysis in the asymptotic critical regime. Indeed, we noticed that the estimated effective critical exponents ( and ) almost remain same in the studied reduced temperature rangepdbabu ; kaul . The observed discrepancy from the 3D Heisenberg values may result from two possibilities: (i) exchange interaction of extended-type, i.e., interactions beyond the nearest neighbor and (ii) magnetic inhomogeneity in the system. Hence, it is necessary to understand the nature as well as the range of interaction in this material. Renormalization-group analysis has shown that the isotropic long-range dipolar interaction can drive a crossover towards mean-field exponents at a certain value of the reduced temperature, which depends on the strength of the dipolar interactions kaul ; aaharony ; fisherme . The analysis also shows that the new crossover exponents are only slightly different from the Heisenberg ones and shifted towards the mean-field values. The general theory of crossover phenomenon gives the crossover temperature as \varepsilon_{co}$$\approx g_{d}$${}^{1/\phi_{d}}, where is the strenght of the dipolar intearction and is the crossover exponent KRied . can be estimated as g_{d}$$\approx$$\frac{0.87}{T_{C}}$$\frac{p^{2}\theta}{v} aaharony ; fisherme ; KRied , where , is a correction factor, and is the elementary cell volume in Å3. For the present material, =196.93 Å3 and for spin =1/2, can be taken as 0.68 KRied . Using these parameters, we get =2.42610*-5* for Co2TiGe. Taking the crossover exponent as the Hesisenberg exponent =1.386, we estimate the crossover temperature to be =4.7810*-4* kaul ; KRied . Only for \varepsilon$$<<$$\varepsilon_{co}, the dipolar inteaction results in a new regime with exponents shifted towards mean-field values. Therefore, the isotropic dipolar picture is not consistent with the present scenario since the estimated for Co2TiGe is well below the investigated asymptotic critical regime of 1.3101.4 10*-2*. However, the extended type of interaction can also arise from isotropic exchange interaction between spins involving itinerant electrons that decays spatially as , where is the effective dimensionality of the spin interaction and is a measure of the range of the exchange interaction Fisher . In a 3D system with isotropic spins, the Heisenberg model is realized for 2 such that decreases with faster than . Mean-field exponents are realized, when with 4.5. For 3/22, the system belongs to a different class with critical exponents that take intermediate values depending on . Given the dimensionality of lattice () and spin (), can be estimated using the renormalization-group approach which gives Fisher ,
[TABLE]
where and . Using Eq. (9), for and , we have obtained the experimentally estimated value of . However, for , scaling relationsEugene yield which is deviated from the 3D Heisenberg value towards the mean-field one but not as much deviated as observed experimentally (Table I). The value of suggests the extended nature of the exchange interaction, which could be due to the presence of RKKY-type long-range interaction between spins involving itinerant electron as argued by others and discussed before Trudel . The values of the asymptotic critical exponents for uniform ferromagnet Ni (Table I) are also shifted towards the mean-field values and the shift is attributed to a crossover to the fixed-point corresponding to the isotropic long-range exchange interactions kaul . But, in Co2TiGe the large deviation in the value of the critical exponent from the 3D Heisenberg can not be explained only on the basis of extended nature of the exchange interaction. It should be mentioned that the mean-field like values of has been reported for several magnetically inhomogeneous systems sjpoon ; nkhan ; rvenka . The amorphous FM Gd80Au20 also shows unusual values of the critical exponents (Table I), all being largely shifted towards the mean-field values, and are attributed to the dilution model for inhomogeneous FMs sjpoon ; hmuller . Magnetic inhomogeneities may lead to critical phenomena those could not be described on the basis of existing universality classes. However, unlike inhomogeneous FMs, the critical exponent shows quite small shift towards the mean-field one for Co2TiGe. Thus, all the obtained values of the critical exponents in Co2TiGe can not be accounted for the magnetic inhomogeneity alone. Therefore, we believe that the observed unconventional critical exponents in Co2TiGe comply with the complex nature of exchange interactions, involving both the short-range and beyond nearest neighbor spin-spin interactions, and also suggest some magnetic inhomogeneity in the system.
IV Conclusion
We elucidate the nature of exchange interactions in the half-mettalic full-Heusler ferromagnet Co2TiGe by comprehensively studying the critical behavior in the vicinity of the ferromagnetic to paramagnetic phase transition at =371.5 K. Although, magnetic Heusler alloys are regarded as an ideal local-moment system and their spin-spin interactions are modeled by a Heisenberg Hamiltonian, we observe a significant deviation in the values of critical exponents from the 3D Heisenberg towards the mean-field ones in Co2TiGe. In particular, the value of the critical exponent is very close to the mean-field one. This implies diminished critical spin-fluctuation in the vicinity of the transition temperature, which could be due to the presence of exchange interactions beyond the nearest neighbor spin as well as magnetic inhomogeneity in the compound. The results are consistent with the recently argued fact that both the short-range direct exchange coupling between nearest neighbor and the indirect coupling between further neighbor spins mediated by long-range RKKY interaction, stabilize the long-range ferromagnetic ordering. The present study provides experimental evidence of the complex nature of exchange interactions in the a Co2-based full Heusler ferromagnet, which does not belong to any known universality class predicted by theory.
V ACKNOWLEDGMENTS
We thanks Mr. Arun Kumar Paul for his help during sample preparation and measurements.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) F. Heusler, Verh. Deutsch. Phys. Ges. 5 , 219 (1903).
- 2(2) P. J. Webster, J. Phys. Chem. Solids 32 , 1221 (1971).
- 3(3) P. J. Webster, and K. R. A. Ziebeck, J. Phys. Chem. Solids 34 , 1647 (1973).
- 4(4) A. Hamzic, R. Asomoza, and I. A. Campbell, J. Phys. F 11 , 1441 (1981)
- 5(5) D. C. Price, J. Phys. F: Met. Phys. 8 933 (1978).
- 6(6) J. Kubler, A. R. William, and C. B. Sommers, Phys. Rev. B 28 , 1745 (1983).
- 7(7) J. L. Morán-Lópeza, R. Rodriguez-Albab, and F. Aguilera-Granjab, J. Magn. Magn. Mater. 131, 417-426 (1994).
- 8(8) Y. Kurtulus, R. Dronskowski, G. D. Samolyuk, and V. P. Antropov, Phys. Rev. B, 71 , 014425 (2005).
