In the present work we present a preliminary experimental and ab-initio study of a series of compounds with nominal composition Sm1-xMMxCo5-y-zFeyNiz (x = 0 – 0.7; y = 0.5 – 1.5; z = 0.5 – 1). Compounds with MM content up to 50% were successfully synthesized in CaCu5-type structure, retaining uniaxial magnetocrystalline anisotropy. Magnetization values are within the 85 to 50 Am2/kg range, with Ni content playing a significant role in stabilization of the structure at the cost of magnetization weakening. Theoretical calculations predict higher magnetization values, especially in the case of SmCo2.5Fe1.5Ni compound with a value close to that of SmCo5.
I. INTRODUCTION
High demand predictions for PM materials are difficult to be met by the current supply chain, mostly due to the scarcity of some critical raw materials.1–3 A safe solution to this may arise from development of new compounds via targeted modification of established materials, for example reducing Co content. In the same scheme, we can investigate the possibility of using as substitutes raw materials with abundance greater than their demand, like Ce and La, a path already utilized in the case of Nd-Fe-B based PM.4,5 Sm-Co, a system with excellent magnetic properties has been investigated towards reduction of both constituent elements, especially Co but also for replacing Sm with other rare-earth (RE) atoms.6–14 Specific properties are positively affected and regarding the supply cost, overall merit of the material is improved. Another important field towards development of new materials suitable for PM applications is the theoretical understanding of the underlying physics. Ab-initio methods are widely used, and especially DFT-based calculations are usually successful in RE-TM intermetallic compounds, however, there are still issues to be clarified.7,15–17
In this work we have examined the effect of Sm-Co replacement for Mischmetal-LaCe3 and Fe, Ni respectively on the fundamental properties of the SmCo5 intermetallic, towards a material with better overall merit, including the cost-criticality, for PM applications. We combine theoretical ab-initio calculations and experimental results, with the additional goal to improve the understanding of the fundamental physics of the system.
II. EXPERIMENTAL-COMPUTATIONAL DETAILS
Samples of nominal composition Sm1-xMMxCo5-y-zFeyNiz (x = 0 – 0.75; y = 0.5 – 1.5; z = 0.5 – 1) were prepared by Ar arc-melting, followed by thermal homogenization at 900-1200 K. Alloys were manually ground down to less than 50 μm using a sieve of 50 mesh. Structural characterization of the powders was carried out using X-ray diffraction (XRD) patterns recorded with a SIEMENS D500 diffractometer (Cu-Kα radiation). XRD patterns of magnetically aligned powder samples with the alignment direction normal to the sample holder were recorded in order to study the magnetic anisotropy of the compounds. The magnetization measurements were performed on magnetically oriented powders fixed in epoxy resin at room temperature using a PAR 155 vibrating sample magnetometer (VSM) in applied fields up to 2 T. The composition was determined by a scanning electron microscope (SEM), equipped with an electron microprobe analyzer (EDAX).
DFT (ab initio) calculations were performed by using VASP package,18 with the PAW version of the PBE pseudopotentials.19–21 All calculations were fully non-colinear, with spin-orbit coupling and aspherical corrections inside the PAW spheres and to Kohn-Sham potential being all explicitly considered. For the treatment of strongly correlated electrons, the rotationally invariant DFT+U approach was employed.22 A value of Ueff = U-J = 4.7 eV was used for the 4f electrons of Sm-Ce, along with a value of Ueff = 2.22 eV for 3d electrons of La, Co, Fe, Ni.23,24 Kinetic energy cutoff for the wave functions was 500 eV. Sm1-xMMxCo5-y-zFeyNiz structures were represented by using 96 atom, 2⨯2⨯4 supercells derived from the primitive cell of the hexagonal P6/mmm, CaCu5-type structure of SmCo5. It has been previously computationally demonstrated that the substitution of a single Co atom with a Ni one is energetically preferable to occur on the 3g sites; if more Co atoms are substituted with Ni, then one substitution per unit cell occurs at a 3g site and the rest at 2c.24 Fe atoms favor the occupation of 3g sites in the CaCu5-type structure.10,16 Therefore, we considered substitution of Co atoms with Fe ones on the 3g sites and with Ni ones at the 2c sites. In all cases calculations simulated the stoichiometries from first principles, independently of the experimental results.
III. RESULTS AND DISCUSSION
In Fig. 1 X-ray powder patterns for all stoichiometries are presented; in the case of Sm0.5MM0.5Co4Ni a pattern of epoxy-oriented sample is also included. In almost all cases the desired CaCu5-type structure is formed, with obvious exception of Sm0.3MM0.7Co3FeNi compound. Phenomenological observations were confirmed by Rietveld analysis, Table I. We have assumed that Fe atoms replace Co in 3g positions and Ni atoms in 2c positions as in the case of ab-initio calculations and we have allowed for transition metal content which in all cases but Sm0.3MM0.7Co3FeNi was below 2 wt.%; in the latter case the 1:5 and 2:7 mass fractions were almost equal (41 and 49 wt.%) while excess transition metal accounts for the rest. The basal plane unit cell parameter is slightly smaller than the reported values for SmCo5 (a = b = 0.5002 nm; c = 0.3961 nm)25 while the perpendicular is slightly larger. This results in a larger c/a ratio. Materials retain the uniaxial character of the magnetocrystalline anisotropy, as observed by using epoxy-oriented samples. In Fig. 1 the corresponding plot for Sm0.5MM0.5Co4Ni is presented as an example, it is evident that only the (0 0 1) and (0 0 2) reflections are present in the oriented sample. The successful alloying of the materials was confirmed by scanning electron microscopy and EDAX analysis (Fig. 2). The refined composition as obtained is close to the starting composition, for example for spot nr. 10 (Fig.2) was found the stoichiometry Sm0.7MM0.3Co2.95 Fe1.04Ni0.99.
X-ray diffraction patterns (Cu Kα radiation) of Sm1-xMMxCo5-y-zFeyNiz samples. 2:7 phase peaks are marked with an asterisk.
X-ray diffraction patterns (Cu Kα radiation) of Sm1-xMMxCo5-y-zFeyNiz samples. 2:7 phase peaks are marked with an asterisk.
Experimental structural and magnetic properties of Sm1-xMMxCo5-y-zFeyNiz samples.
Nominal stoichiometry . | a = b (nm) . | c (nm) . | c/a . | MS (Am2/kg) . | μB/f.u. . |
---|---|---|---|---|---|
Sm0.5MM0.5Co4Ni | 0.4968(1) | 0.4004(1) | 0.806 | 57.0 | 4.486 |
Sm0.5MM0.5Co4Fe0.5Ni0.5 | 0.4987(1) | 0.4013(1) | 0.805 | 85.0 | 6.667 |
Sm0.5MM0.5Co3FeNi | 0.4976(1) | 0.4033(1) | 0.810 | 71.8 | 5.611 |
Sm0.7MM0.3Co3FeNi | 0.4973(1) | 0.4029(1) | 0.810 | 67.6 | 5.308 |
Sm0.3MM0.7Co3FeNi (1:5) | 0.4973(1) | 0.4041(1) | 0.813 | 72.9 | 5.669 |
(2:7) | 0.5063(2) | 2.4407(9) | |||
SmCo2.5Fe1.5Ni | 0.4964(1) | 0.4047(1) | 0.815 | 51.8 | 3.818 |
Nominal stoichiometry . | a = b (nm) . | c (nm) . | c/a . | MS (Am2/kg) . | μB/f.u. . |
---|---|---|---|---|---|
Sm0.5MM0.5Co4Ni | 0.4968(1) | 0.4004(1) | 0.806 | 57.0 | 4.486 |
Sm0.5MM0.5Co4Fe0.5Ni0.5 | 0.4987(1) | 0.4013(1) | 0.805 | 85.0 | 6.667 |
Sm0.5MM0.5Co3FeNi | 0.4976(1) | 0.4033(1) | 0.810 | 71.8 | 5.611 |
Sm0.7MM0.3Co3FeNi | 0.4973(1) | 0.4029(1) | 0.810 | 67.6 | 5.308 |
Sm0.3MM0.7Co3FeNi (1:5) | 0.4973(1) | 0.4041(1) | 0.813 | 72.9 | 5.669 |
(2:7) | 0.5063(2) | 2.4407(9) | |||
SmCo2.5Fe1.5Ni | 0.4964(1) | 0.4047(1) | 0.815 | 51.8 | 3.818 |
SEM image of Sm0.7MM0.3Co3FeNi sample; the twelve black crosses are the sampling points.
SEM image of Sm0.7MM0.3Co3FeNi sample; the twelve black crosses are the sampling points.
The magnetic properties of the materials were investigated by means of RT isotherm loops in fields up to 2 T in powder samples of size 50 μm or less with random orientation. The corresponding curves are presented in Fig. 3 and magnetization is included in Table I. In Fig. 3 we distinguish between x = 0.5 samples and the rest. In the first case, magnetization is mostly reduced by Ni content, the compound with less Ni (y = z = 0.5) presents the highest magnetization. For the other two samples where Ni content is the same the material with no Fe present lower magnetization than itscounterpart with 20% replacement of Co by Fe. These results are consistent with previous reports in the literature.7,16,24 Doping with Ni seems obscure due to its negative effect on magnetic properties, however, there are reports which connect the stability of the Fe-doped compound with Ni inclusion as well.26 The rest of the samples present unexceptional magnetization values, below 70 Am2/kg; Sm0.7MM0.3Co3FeNi is close to Sm0.5MM0.5Co3FeNi, as expected from the similar stoichiometry. The value for Sm0.3MM0.7Co3FeNi is presented as reference since the it is affected by sample’s structural inhomogeneity. Finally, SmCo2.5Fe1.5Ni compound has very low magnetization, although in this case the negative -in general-effect of MM in absent. The results of the calculated structural parameters, total magnetization per formula unit and formation enthalpy are presented in Table II. The multiple substitutions of Sm with the MM alloy and Co with Fe and Ni were found to lead to a 2.5-3.4% increase of the in-plane lattice constant and a 9.4-11.7% increase of the c/a ratio with respect to the calculated values for the SmCo5 compound. Moreover, the introduction of the substituents was found to lead to the presence of a small in-plane structural distortion (0.27-1.82%) which can be detrimental for the magnetic properties; a is not exactly equal to b within the hexagonal plane.27 The introduction of MM at 25% concentration was found to lead to the smaller value of in-plane distortion but also to a reduction of the total magnetization. The substitution of Sm with MM was found to lead to a reduction of the total magnetization by ∼1 μB/f.u. in all cases.
RT hysteresis loops of Sm1-xMMxCo5-y-zFeyNiz powder samples; x = 0.5 (a) and x = 0, 0.3, 0.7 (b).
RT hysteresis loops of Sm1-xMMxCo5-y-zFeyNiz powder samples; x = 0.5 (a) and x = 0, 0.3, 0.7 (b).
Lattice constant a, basal plane distortion, c/a ratio, magnetization per formula unit, m and formation enthalpy, H, for the considered stoichiometries, as calculated from the DFT simulations. In plane anisotropy was defined as (a-b)/a * 100%, a and b are the lattice constants.
Nominal stoichiometry . | a (nm) . | In plane distortion (%) . | c/a . | MS (μB/f.u.) . | H (eV/f.u.) . |
---|---|---|---|---|---|
Sm0.5MM0.5Co4Ni | 0.4917 | 0.98 | 0.847 | 6.093 | −0.826 |
Sm0.5MM0.5Co4Fe0.5Ni0.5 | 0.4939 | 0.85 | 0.843 | 7.459 | −0.681 |
Sm0.5MM0.5Co3FeNi | 0.4961 | 1.62 | 0.839 | 7.384 | −0.751 |
Sm0.7MM0.3Co3FeNi | 0.4942 | 1.82 | 0.840 | 7.524 | −0.825 |
Sm0.3MM0.7Co3FeNi | 0.4945 | 0.27 | 0.845 | 7.119 | −0.713 |
SmCo2.5Fe1.5Ni | 0.4961 | 1.08 | 0.830 | 8.336 | −0.840 |
SmCo5 | 0.5100 | 0.00 | 0.755 | 8.535 | −0.508 |
Nominal stoichiometry . | a (nm) . | In plane distortion (%) . | c/a . | MS (μB/f.u.) . | H (eV/f.u.) . |
---|---|---|---|---|---|
Sm0.5MM0.5Co4Ni | 0.4917 | 0.98 | 0.847 | 6.093 | −0.826 |
Sm0.5MM0.5Co4Fe0.5Ni0.5 | 0.4939 | 0.85 | 0.843 | 7.459 | −0.681 |
Sm0.5MM0.5Co3FeNi | 0.4961 | 1.62 | 0.839 | 7.384 | −0.751 |
Sm0.7MM0.3Co3FeNi | 0.4942 | 1.82 | 0.840 | 7.524 | −0.825 |
Sm0.3MM0.7Co3FeNi | 0.4945 | 0.27 | 0.845 | 7.119 | −0.713 |
SmCo2.5Fe1.5Ni | 0.4961 | 1.08 | 0.830 | 8.336 | −0.840 |
SmCo5 | 0.5100 | 0.00 | 0.755 | 8.535 | −0.508 |
Also, from Table III, it is evident that MM (LaCe3) tends to slightly destabilize the structure, as Ce atoms present lower values of energy per atom than Sm in their groundstate structures. A low energy per atom value in this context can be interpreted as a measurement of the likelihood of an atom to “prefer” to form its corresponding groundstate structure phase, rather than an alloy or an intermetallic compound. The introduction of Ni significantly increases the stability of the structures, as it leads to a reduction of the formation enthalpies with respect to SmCo5, due to its higher value of energy per atom in its groundstate structure, compared to Co, Table III.
Calculated energies per atom for Sm, La, Ce, Co, Fe and Ni in their respective groundstate structures within the DFT+U framework as used in this work.
Atom . | E/atom (eV) . | Atom . | E/atom (eV) . |
---|---|---|---|
Sm | −4.959 | Co | −5.222 |
La | −4.353 | Fe | −6.485 |
Ce | −5.523 | Ni | −3.855 |
Atom . | E/atom (eV) . | Atom . | E/atom (eV) . |
---|---|---|---|
Sm | −4.959 | Co | −5.222 |
La | −4.353 | Fe | −6.485 |
Ce | −5.523 | Ni | −3.855 |
The average values of magnetic moments per each atom type of the rare earth and transition metal elements are shown in Table IV. Strong antiferromagnetic alignment of the Ce atom is evident and can explain the decrease of the total magnetic moments by ∼1 μB/f. u. in the MM containing stoichiometries; Ce has a 4f electron with itinerant character which interacts with Co 3d electrons, modifying their weight at the Fermi level.28 The introduction of both transition metal substitutes leads to a slight decrease of the average atomic magnetic moments of Co atoms by 0.1-0.2 μB/atom. However, the addition of Fe is strongly beneficial for the hard magnetic properties as its magnetic moment appears ∼1.75 times larger than the one of Co. Its addition, though, is known to destabilize the structure, fact which agrees with its corresponding calculated energy per atom in its groundstate structure. This destabilization appears to be eliminated with the addition of Ni, at the cost of reducing the total magnetization, as Ni atoms present very low atomic magnetic moments. The presence of Fe at concentrations ≥ 1 per formula unit appears also to be associated to an increase of the magnetic moments of Ni atoms.
Average values of magnetic moments per atom, calculated from the DFT simulations. Negative values of atomic magnetization indicate antiferromagnetic alignment with respect to the positive ones.
. | Sm . | La . | Ce . | Co (2c) . | Co (3g) . | Fe . | Ni . |
---|---|---|---|---|---|---|---|
Sm0.5MM0.5Co4Ni | −0.295 | −0.267 | −1.191 | 1.646 | 1.615 | - | 0.228 |
Sm0.5MM0.5Co4Fe0.5Ni0.5 | −0.298 | −0.271 | −1.213 | 1.594 | 1.657 | 2.877 | 0.252 |
Sm0.5MM0.5Co3FeNi | −0.295 | −0.279 | −1.187 | 1.621 | 1.604 | 2.842 | 0.340 |
Sm0.7MM0.3Co3FeNi | −0.300 | −0.278 | −1.215 | 1.606 | 1.600 | 2.833 | 0.355 |
Sm0.3MM0.7Co3FeNi | −0.313 | −0.295 | −1.223 | 1.602 | 1.589 | 2.841 | 0.320 |
SmCo2.5Fe1.5Ni | −0.297 | - | - | 1.616 | 1.591 | 2.818 | 0.403 |
SmCo5 | −0.300 | - | - | 1.740 | 1.785 | - | - |
. | Sm . | La . | Ce . | Co (2c) . | Co (3g) . | Fe . | Ni . |
---|---|---|---|---|---|---|---|
Sm0.5MM0.5Co4Ni | −0.295 | −0.267 | −1.191 | 1.646 | 1.615 | - | 0.228 |
Sm0.5MM0.5Co4Fe0.5Ni0.5 | −0.298 | −0.271 | −1.213 | 1.594 | 1.657 | 2.877 | 0.252 |
Sm0.5MM0.5Co3FeNi | −0.295 | −0.279 | −1.187 | 1.621 | 1.604 | 2.842 | 0.340 |
Sm0.7MM0.3Co3FeNi | −0.300 | −0.278 | −1.215 | 1.606 | 1.600 | 2.833 | 0.355 |
Sm0.3MM0.7Co3FeNi | −0.313 | −0.295 | −1.223 | 1.602 | 1.589 | 2.841 | 0.320 |
SmCo2.5Fe1.5Ni | −0.297 | - | - | 1.616 | 1.591 | 2.818 | 0.403 |
SmCo5 | −0.300 | - | - | 1.740 | 1.785 | - | - |
Overall, according to the predictions of our ab initio simulations, the SmCo2.5Fe1.5Ni stoichiometry presents a large value of total magnetization which is comparable to the respective value of SmCo5, with the partial substitution of Co with both Fe and Ni leading to a favorable interplay between stability and magnetization and the additional benefit of the reduction of the content of Co, which is regarded as a critical material. Further substitution of the also critical Sm with the MM alloy can lead to a further reduction of cost by reducing raw material and processing expenses, but seems to inevitably lead to a drop of the total magnetization by ∼ 1 μB/f. u.
IV. CONCLUSIONS
In this work we present structural and magnetic properties of nominal composition Sm1-xMMxCo5-y-zFeyNiz (x = 0 – 0.75; y = 0.5 – 1.5; z = 0.5 - 1) from both experimental results and ab-initio DFT calculations. All studied stoichiometries with MM content up to 50% can be synthesized in bulk form with conventional metallurgy techniques.
Ni probably has a positive effect in structure stabilization as also depicted by ab-initio calculations, at the cost of weakening the overall magnetization. In both structural and magnetic properties ab-initio calculations and experimental results are in good agreement, calculated magnetization is larger than experimental results; suggested DFT parameters are suitable for RE-TM intermetallics. The suggested substitutional approach may provide PM suitable for some applications with reduced cost and criticality compared to basic SmCo5 system.
ACKNOWLEDGMENTS
S. Giaremis acknowledges support for the research work by the Hellenic Foundation for Research and Innovation (HFRI) under the HFRI PhD Fellowship grant (Fellowship Number: 962). Furthermore, this work was supported by computational time granted from the Greek Research & Technology Network (GRNET) in the “ARIS” National HPC infrastructure under the project NOUS (pr010034) and the Aristotle University of Thessaloniki (AUTh) HPC Infrastructure and Resources.
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.