Introduction

Nd–Fe–B based alloys have been widely used as the most powerful permanent magnets for the last 30 years. The replacement of Nd by Ce element has been suggested as an alternative to reduce the global cost of the permanent magnets1,2. Even though the intrinsic properties of the Ce2Fe14B phase (Js = 1.17 T, Ha = 30 kOe and Tc = 442 K) are relatively lower compared to those of Nd2Fe14B (2:14:1) phases3. Some experimental studies suggested that up to 30 at.% substitution of Ce for Nd in Nd-Fe-B alloys did not destroy the magnetic properties, but further increase of Ce content would reduce the hard magnetic properties significantly4,5. Particularly, several topics deserve more attention in the case of Ce-substituted Nd-Fe-B magnets. First, when the Nd is partially substituted by other rare earth (RE), the RE elements show their preference to enter into certain phase due to substitution energy variation. It was reported that the Y, Tb and Dy elements prefer to enter the 2:14:1 phase but avoid the Nd-rich phase, while the La and Ce ones exhibit the opposite preference6,7,8. Thus, it is crucial to explore the influence of the Ce elements on magnetic properties and to control their distribution in the Ce-substituted magnets. Second, the previous studies have shown that REFe2 phase is easy to precipitate in the Ce-Fe-B and Ce-containing Nd-Fe-B alloys9,10,11. The CeFe2 phase behaves as a soft magnetic phase below 230 (±2) K12 and its presence is harmful to the hard magnetic properties. On the other hand, REFe2 phase also brings some benefits for the sintered magnets. Zhang et al.10 reported that the REFe2 phase could improve the wettability of 2:14:1 phase when the magnets are sintered at temperature above the melting point of REFe2 phase. Our previous work also suggested that Ce atom prefer to enter into the 1:2 phase, which leads the element segregation and increases the Nd content in the (NdCe)2Fe14B phase8. Third, Capehart et al. reported a mixed-valent state with the coexistence of Ce3+ and Ce4+ atoms in the Ce2Fe14B phase based on the evidence provided by Ce L3-edge x-ray absorption near-edge structure measurement13. Their steric variation combined with different electron density give rise to the valence state dependent magnetic and metallurgical behaviors of Ce-Nd-Fe-B alloys14,15. The last but not the least, an abnormal increase of coercivity has been observed in the Ce-Nd-Fe-B alloys with around 25% Ce-substitution, and several researchers have attributed it to the phase separation2,16,17,18. However, to the best of our knowledge, no microstructural evidence of the phase separation has been reported yet. While, the recent study on the Ce-substituted single crystal indicated a disagreement19, which causes more doubts on the phase separation.

In this work, aiming at understanding the element segregation and phase separation in the Ce-substituted alloys, the microstructure and the magnetic properties of the directly cast Ce-substituted Nd-Fe-B based alloys have been investigated. The segregation of RE elements, the roles of REFe2 phase, the mixed valence state of Ce, and the phase separation of 2:14:1 phase have been carefully investigated and discussed. In particular, the phase separation is confirmed directly by novel microstructural analysis. Here, directly cast (Nd1−xCex)25Fe40Co20Al4B11 alloys are chosen because their grain sizes are in micrometer scale, which can be easily characterized by scanning electron microscopy and energy dispersive spectrometry. In addition, those cast alloys have high coercivities, which is beneficial for revealing the relationships between the magnetic properties and microstructures20,21. The present work may help to further understand the structure and magnetic properties of Ce-containing Nd-Fe-B alloys.

Results

Phase constitution and microstructure

The XRD results of selected (Nd1−xCex)25Fe40Co20Al4B11 samples with various Ce substitutions for Nd are showed in Fig. 1a. As discussed in our previous work21, the x = 0 sample is composed of Nd2(FeCoAl)14B phase (2:14:1 phase) with a space group of P42/mnm and Nd1+ε(FeCo)4B4 phase (1:4:4 phase) with a space group of P42/n. With substitution of Ce for Nd, the 1:4:4 phase is replaced by the (NdCe)(FeCo)2 phase (1:2 phase) with a space group of Fd\(\overline{{\rm{3}}}\)m. With increasing Ce content, the relative intensity of the 1:2 phase increases, indicating that the content of the 1:2 phase in the sample increases. Figure 1b shows the XRD patterns in the 2θ-range of 40–43°. The marked dotted lines in Fig. 1b show that the peaks shift to higher angles with the increasing Ce content. Refined values of cell parameters of 2:14:1 and 1:2 phases are plotted with x in Fig. 1c and d, respectively. The lattice parameters of the 2:14:1 phase and the 1:2 phase decrease monotonically, but not linearly, with the increasing Ce content. For the 2:14:1 phase, the values of a, c, and V decrease slowly from x = 0 to x = 0.4 at first, and then drop more quickly for x > 0.4. For the 1:2 phase, an opposite change tendency is observed, and the values of a and V decrease quickly from x = 0 to x = 0.4 but decrease slowly between x = 0.5 and 0.7. From these data, it could be expected that the magnetic properties would have changed accordingly since the strength of the magnetic coupling and anisotropy are correlated to the distance between magnetic atoms, which is dependent on the cell parameters.

Figure 1
figure 1

XRD patterns for the as-cast (Nd1−xCex)25Fe40Co20Al4B11 alloys (a), the enlarged patterns in the range of 40–43° (b), and the calculated cell parameters of 2:14:1 phase (c) and 1:2 phase (d) changed with Ce content x.

The back-scattered SEM images for the selected samples are shown in Fig. 2a–f. Three different contrasts, marked as area A, B, and C in all figures, can be observed. For the x = 0 sample (Fig. 2a) without Ce substitution, as discussed in our previous work20, area A is represented for 2:14:1 phase, B is for Nd-rich phase while C (dark gray) for 1:4:4 phase. For the Ce-substituted alloys, according to the EDS results shown in the supplementary materials as Table S1, the ratio of rare earth Nd + Ce content over transition metal Fe + Co content (RE/TMatom) is near to 1:7 for point A, 2:1 for point B, and 1:2 for point C, which should correspond to 2:14:1 phase, RE-rich phase, and 1:2 phase, respectively. It is found that when the alloy was substituted by Ce, the 1:4:4 phase disappeared while the 1:2 phase appeared. The area fraction of the 1:2 phase increases with the increasing Ce content, in agreement with the XRD results. The distribution of the 1:2 phase also changes with Ce content. For the x = 0.1–0.4 alloys, as displayed in Fig. 2b–d, the 1:2 phases are randomly distributed. In the case of x = 0.5–0.7, as illustrated in Fig. 2e and f, the 1:2 phase prefers to stay around the 2:14:1 phase, giving rise to a ‘core-shell’ like structure.

Figure 2
figure 2

Back-scattered SEM images for the as-cast (Nd1−xCex)25Fe40Co20Al4B11 alloys with (a) x = 0, (b) x = 0.2, (c) x = 0.3, (d) x = 0.4, (e) x = 0.5, and (f) x = 0.7, and the dependences of RCe/RE exp/RCe/RE nom (g) and RFe/TM exp/RFe/TM nom (h) ratios on Ce content x for different phases.

Based on the EDS results, the experimental Ce/RE and Fe/TM ratios in the 2:14:1, RE-rich, and 1:2 phases are compared with nominal Ce/RE and Fe/TM ratios, the RCe/RE exp/RCe/RE nom and RFe/TM exp/RFe/TM nom changing with Ce content x were plotted in Fig. 2g and h, respectively. In Fig. 2g, the line y = 1 is the standard line. For the 2:14:1 phase, all the values of RCe/RE exp/RCe/RE nom are lower than 1, which means that the Nd atoms prefer to enter into the 2:14:1 phase. By contrast, all the values of RCe/RE exp/RCe/RE nom in the 1:2 phase are higher than 1, indicating a Ce segregation. For the RE-rich phase, in the x = 0.1–0.3 alloys, the values of RCe/RE exp/RCe/RE nom are close to 1, but it starts to be lower for x > 0.3, suggesting the Nd segregation. The results thus indicate a strong RE element segregation in the alloys. Ce atoms prefer to enter into 1:2 phase, but expel from the 2:14:1 phase, and the distribution in RE-rich phase is dependent on the Ce content. Figure 2h presents the RFe/TM exp/RFe/TM nom in the various phases, and the y = 1 is the standard line. One could find that the segregation of TM also occurs. Fe atoms prefer to enter into 2:14:1 phase, but to expel from both RE-rich and 1:2 phases. On the other hand, the Fe/TM ratios in those three phases are independent of the Ce content.

The selected samples are characterized by TEM. Figure 3a illustrates a micrometer sized grain (labelled 1) and intergranular phases (labelled 3) for the x = 0 sample. The select area diffraction patterns (SADP) corresponding to areas 1 and 3 confirm the microstructure of 2:14:1 grains surrounded by amorphous RE-rich phase. Figure 3b shows the interface (labelled 2) between the 2:14:1 grain and the intergranular phase. The smooth interfacial structure is beneficial for the coercivity22. Similarly, the TEM image obtained from the x = 0.3 sample indicates that the micrometer sized grain (labelled 4) is surrounded by intergranular phase (labelled 6) (see Fig. 3c). However, according to the SADP image for the selected area 6, the intergranular phase could be indexed as nanocrystalline grains embedded in the amorphous matrix. Moreover, in the SADP image for selected area 4, except the main pattern of 2:14:1 phase with a \([\bar{2}2\bar{1}]\) crystal belt, a satellite pattern or, more concisely, a steak-line pattern, is also observed. Such a feature should be related to the phase separation2, as discussed later. The interfacial area labelled 5 in the HRTEM graph (Fig. 3d) shows the smooth grain boundary with the amorphous RE-rich phase, similar to that in the x = 0 sample. Figure 3e clearly shows that the structure of the intergranular phase for the x = 0.7 sample is more complicated than those of the x = 0 and x = 0.3 samples. The SADP image of area 9 suggest that the intergranular phase is mainly composed of nanocrystals and the 2:14:1 grains are surrounded by crystalline grains, consistently with the SEM results. Figure 3f is the HRTEM micrograph for the interfacial area 8, the right part of the graph shows an interplanar distance of 0.264 nm, corresponding to the [220] plane of 1:2 phase. It thus confirms that the surrounded phases are 1:2 phases, as discussed in Fig. 2, which form a ‘core-shell’ like structure with the 2:14:1 grains.

Figure 3
figure 3

TEM images for the selected as-cast (Nd1−xCex)25Fe40Co20Al4B11 alloys (a), (b) x = 0; (c,d) x = 0.3; (e,f) x = 0.7.

Magnetic properties and hyperfine structure

The demagnetization curves for the alloys with various Ce contents are shown in Fig. 4a. All alloys exhibit a single phase-like hard magnetic behavior, and the demagnetization curve shrinks rapidly after 40% of Nd is substituted. As shown in Fig. 4b, for the low Ce contents (x < 0.2), both saturation magnetization Ms and remanence Mr remain almost constant, but drop monotonically for the increasing Ce content (x > 0.2). Since all samples are magnetically isotropic, the Mr/Ms values all keep close to 1/2. It has been demonstrated that the volume fraction of 1:2 phase and the Ce/RE% of 2:14:1 phase increase with the increasing Ce content, which should be the reasons for the degradation of both Mr and Ms. The (BH)max, not shown here, also shows the same trend as the change of Ms. Interestingly, intrinsic coercivity Hci increases unexpectedly at x = 0.3, similar to those previously observed by various researchers2,16. It increases from 641 kA/m for x = 0.2 to 863 kA/m for x = 0.3, which is still higher than the 829 kA/m of the x = 0.1 sample. After x > 0.4, the rapid decrease of the coercivity could be attributed to the “core-shell” like structure discussed above22.

Figure 4
figure 4

The demagnetization curves (a) and magnetic properties (b) for the as-cast (Nd1−xCex)25Fe40Co20Al4B11 alloys, and Mössbauer spectra at 300 K for the selected as-cast (Nd1−xCex)25Fe40Co20Al4B11 alloys (c).

57Fe Mössbauer spectrometry was employed to investigate the hyperfine structures of the different alloys and to discriminate their Fe-based phases. The spectra registered at 300 K are shown in Fig. 4c and the refined values of hyperfine parameters are listed in Table 1. The hyperfine structures for all samples are quite complicated and have to be described by at least three main components, including one consisting of six sextets to describe the 2:14:1 phase, and two quadrupolar doublets assigned to the paramagnetic RE-rich phase and the 1:2 phase. The assignment of the six magnetic components is based upon the Wigner-Seitz analysis of the RE2Fel4B structure23, which is even more complicated in the presence of Co24. The refined values of hyperfine parameters of the CeFe2 phase are consistent with those previously reported11,25 and the remaining paramagnetic component was assigned to the RE-rich phase. From Table 1, the content of the 1:2 phase increases with the increasing Ce content at the expense of the 2:14:1 phase, while that of the RE-rich phase remains almost independent with Ce content. The 2:14:1 phase is the only magnetic phase in the alloys at 300 K and its magnetic moments are from 14 transition metal (TM) atoms Fe and Co, and 2 rare earth (RE) atoms Nd and Ce. It was reported that the magnetic moment of the Ce atom in the 2:14:1 structure was extremely low, ranging 0.1 ~ 0.5 μB11,15. According to the mean hyperfine field of the Fe atom (Table 1), and by using the rule of 15 T/μB23, the mean magnetic moment of Fe atoms can be estimated at 1.70 μB, 1.68 μB, 1.67 μB, and 1.55 μB (±0.01 μB) for the x = 0.1, 0.2, 0.3, and 0.5 alloys, respectively23. It is clear that, with Ce substitution for Nd, the magnetic moment exhibits a minor decrease from x = 0.1 to 0.3, then a rapid decrease afterward. Thus, due to the decrease of the Fe moment with Ce substitution and the substitution of magnetic Nd (3.2 μB) by near non-magnetic Ce3, the magnetic moment of the 2:14:1 phase decreases with the increasing Ce content, especially for x > 0.3.

Table 1 Refined values of the hyperfine parameters at 300 K for the selected as-cast (Nd1−xCex)25Fe40Co20Al4B11 alloys.

Discussions

To take the advantage of the RE elements doping in the Nd-Fe-B alloys, it is important to understand their partitioning in the main phase and the intergranular phase. Recently, Li et al.26 reported that Nd atoms tend to diffuse into the main phase during the Ce substitution in the strip cast (Nd1−xCex)30Fe69B (wt%) alloys, in agreement with our results in Fig. 2g.

The reason for Nd segregation in the 2:14:1 phase can be contributed to the substitution energy of RE (Esub) in 2:14:1 phase, which can be described as equation 1:

$${E}_{sub}={E}_{{(N{d}_{1-x}R{E}_{x})}_{2}F{e}_{14}B}-{E}_{N{d}_{2}F{e}_{14}B}-2x{E}_{RE}+2x{E}_{ND}$$
(1)

where E(Nd1−xREx)2Fe14B represents the total energy at 0 K. The negative Esub means that the substituted RE prefers to enter into 2:14:1 phase, while the positive Esub indicates that the substitution will be expelled from the 2:14:1 phase. As reported by Liu et al.6,7, the Esub for Y, Dy and Tb were negative, but for La it was positive (0.41 eV/atom). Since La and Ce elements are expected to have similar properties, and both La3+ and Ce3+ ionic radii are larger than that of Nd3+. Hence, the Esub for Ce should be positive, and then it will be expelled from the 2:14:1 phase, as evidenced from the present experimental results.

However, this study gives different results from that of Li et al.26. Firstly, Li et al. reported that the Ce atoms preferred to enter into the RE-rich phase, which is a priori conflicting with the Nd segregation suggested in this work. We suggest here that the existence of REFe2 phase be responsible for Ce segregation. It was established that the large atom radius prevent to obtain stable NdFe2, LaFe2, PrFe2 and YbFe2 phases27. Cannon et al.28 also reported that the NdFe2, PrFe2 and YbFe2 phases could only be formed under high pressure and temperature (80 kbars, 1200 °C), except for CeFe2 phase due to its mixed valence state. The small radius of Ce4+ ion makes it easier to form the CeFe2 phase, which explains the Ce segregation into the 1:2 phase and thus also enhances the Nd segregation into the 2:14:1 phase. Secondly, in the strip cast (Nd1−xCex)30Fe69B(wt.%) alloys, Li et al. found that the Ce/RE ratio in 2:14:1 phase depended linearly on the Ce content, but in our case, as shown in Fig. 2g, the Ce/RE ratio increases more quickly for x > 0.4 in the 2:14:1 phase while the opposite tendency is observed in the 1:2 phase. To explain our result, the increase of the occupation results from the decrease of Esub, which can be achieved effectively by decreasing the cell volume6. In the Ce-based or Ce-substituted 2:14:1 phase with mixed valence state, i.e. the coexisting Ce3+ and Ce4+ states13,17, as Ce4+ ion is smaller than Ce3+ ion, Esub will decrease when more Ce atoms are in the Ce4+ state, and consequently more Ce atoms will enter into 2:14:1 phase.

The changes of Ce/RE ratios in 2:14:1 and 1:2 phases lead to the variations of the cell parameters of both phases. As shown in Fig. 1c and d, the changes of the cell volume V can be correlated to the values of RCe/RE exp/RCe/RE nom displayed in Fig. 2g and h. For the 2:14:1 phase, Ce/RE ratio in the x = 0–0.4 alloys increases more slowly than that in the x = 0.5–0.7 alloys, and as a result, the V2:14:1 of x = 0–0.4 alloys decreases more slowly. For x > 0.4, the quicker decrease of V2:14:1 is related to the quicker increase of Ce/RE ratio and also to that more Ce atoms in Ce4+ state. Similarly, the Fe magnetic moments of the 2:14:1 phase estimated from the Mössbauer results also show the same tendency. The weak decrease of magnetization at the early stage of Ce substitution indicates that the Ce-substituted magnets may be suitable for application.

The most interesting result in the present work is the microstructural evidence for the phase separation. As illustrated in Fig. 3c, the steak line diffraction pattern of 2:14:1 phase in the x = 0.3 sample indicates a two-phase structure in the alloy, which should be responsible for the abnormal increase of coercivity. The mechanism for phase separation in such alloys was investigated by Alam et al.17 and they found that when the Ce substitution exceeded 25%, the chemical inhomogeneity occurred in (NdCe)2Fe14B phase due to the significant difference in Nd2Fe14B and Ce2Fe14B lattice volumes. In our work, the phase separation is observed in the x = 0.3 sample, which has an actual Ce/RE ratio close to 21%. The HRTEM micrograph in Fig. 5a shows that the crystal planes of the 2:14:1 phase are not perfectly aligned for the x = 0.3 sample, different from the typical ones of Nd2Fe14B crystal29,30. The periodic differences are marked with dashed line in Fig. 5a. It can be assumed that each period corresponds to one phase, and the strain between the two neighboring phases creates the dislocations at the interfaces. The difference between those two phases may result from the Ce/RE ratio difference, one is higher than 21% (with lower cell volume) while the other one is lower than 21% (with higher cell volume), but both of them have the same 2:14:1 structure and crystal orientation. Wu et al.31 also reported the similar phenomenon observed at the interfaces of PbTe/PbS phases. Here we found this structure in all 2:14:1 phases of the x = 0.3 sample. The fast Fourier transformation (FFT) graph, shown as inset 1, presents a typical pattern for the 2:14:1 phase and a satellite pattern. Those two patterns can be separated by using ‘apply mask’ in the DigitalMicrograph software into the 2:14:1 pattern (inset 2) and satellite pattern (inset 3), and their inverse FFT graphs are illustrated in Fig. 5b and c, respectively. In Fig. 5c, the boundaries of the two-phase structure are clearly observed and it proves that the dislocations in the boundaries of those two phases are responsible for the steak-line diffraction pattern. To the best our knowledge, this nanoscale periodic two-phase structure with a width around 2 nm has never been reported before.

Figure 5
figure 5

(a) HRTEM image and FFT image (inset 1) of the 2:14:1 phase in the as-cast x = 0.3 sample (the FFT image pattern can be divided into 2 patterns, i.e. a 2:14:1 main pattern (b inset 2) and a satellite pattern (b inset 3)); (b) the inverse FFT image for the main pattern in b inset 2; (c) the inverse FFT image for the satellite pattern in c inset 3; (d) the HRTEM image and FFT image (inset 4) of 2:14:1 phase in the as-cast x = 0.3 sample heat treated at 750 °C for 10 min.

As predicted by Alam et al.17, this inhomogeneous structure would disappear at high temperature. Thus, the x = 0.3 sample was heat treated at 750 °C for 10 min. The HRTEM micrograph for the heat treated sample exhibits perfectly aligned crystal planes (Fig. 5d) and the FFT graph presents the typical single crystal diffraction pattern (inset 4). It means that the two-phase structure disappears after heat treatment. This study well explains why all the reported phase separation phenomena were only observed in the directly quenched alloys2,16,18, but not in the single crystal or sintered samples. Hence, our work will make this topic more conclusive19,32,33.

In summary, the (Nd1−xCex)25Fe40Co20Al4B11 (x = 0–0.7) alloys were prepared by copper mold casting and their composition-microstructure-magnetic properties relationships have been studied in this work in order to better understand the element segregation and phase separation in the Ce-substituted Nd-(Fe,Co)-B based alloys. In particular, the element segregation, mixed Ce valence states, roles of the REFe2 phase, phase separation, and high coercive field in the Ce-substituted alloys have been discussed in details. The cast Ce-substituted alloys are composed of 2:14:1 phase, 1:2 phase and RE-rich phase. The Ce/RE ratio is found lower than the nominal ratio in the 2:14:1 phase but higher in the 1:2 phase, indicating that Ce atoms prefer to enter 1:2 phase and tend to avoid the 2:14:1 phase due to the positive Esub. The weak increase of the Ce/RE ratio in the 2:14:1 at the early stage (x = 0–0.4) contributes to a slow decrease of the cell parameters and Fe magnetic moment for the 2:14:1 phase. But when x exceeds 0.4, the rapid decrease of the cell parameters could be attributed to more tetravalent Ce atoms. It is interesting to find the preliminary microstructural evidence for the phase separation in this work, which should be the reason for the abnormal increase of coercivity in the as-cast x = 0.3 sample.

Methods

The (Nd1−xCex)25Fe40Co20Al4B11 (x = 0–0.7, atomic ratio) alloy ingots were prepared by arc-melting the mixture of pure metals Ce (99.9 wt.%), Nd (99.9 wt.%), Fe (99.9 wt.%), Co (99.9 wt.%), Al (99.9 wt.%), and B (99.9 wt.%) in an argon atmosphere. The ingots were re-melted 5 times to promote chemical composition homogeneity. Bulk rods with 2 mm in diameter were produced by injection casting the molten alloy into the copper mold without cooling water under the protection of Ar atmosphere. The selected x = 0.3 sample was annealed at 750 °C for 10 min in an argon atmosphere.

The magnetic properties were measured using physical property measurement system (PPMS-9, Quantum Design) equipped with a 9 T vibrating sample magnetometer (VSM). The microstructure was examined by transmission electron microscopy (TEM, Tecnai G2 F20 S-TWIN 200 kV) and scanning electron microscope backscattering electron mode (SEM, Nano430, FEI Co, BSE mode) equipped with an energy dispersive spectrometer (EDS). For the EDS results, at least 5 points were measured to get the mean values. The phase structure was determined by X-ray diffraction (XRD, Philips X-pert) with Cu-Kα radiation. The alloys were investigated by 57Fe Mössbauer spectrometry (MS), and the spectra were recorded at 300 K in a transmission geometry using 57Co/Rh γ-ray source mounted on an electromagnetic drive with a triangular velocity form. For both the XRD and MS measurement, the rod-shaped samples had to be crushed into fine powders to prevent from any preferential and thickness effects.