Next Article in Journal
Multi-Channel Expression Recognition Network Based on Channel Weighting
Previous Article in Journal
Advanced Assessment of Biomass Materials Degradation in Pneumatic Conveying Systems: Challenges and Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermodynamic Modeling of Formation Enthalpies of Amorphous and Crystalline Phases in Zr, Nd, and Ce-Substituted Fe-Si Systems

by
Mieszko Kołodziej
1,2 and
Zbigniew Śniadecki
1,*
1
Institute of Molecular Physics, Polish Academy of Sciences, M. Smoluchowskiego 17, 60-179 Poznań, Poland
2
NanoBioMedical Centre, Adam Mickiewicz University, Wszechnicy Piastowskiej 3, 61-614 Poznań, Poland
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(3), 1966; https://doi.org/10.3390/app13031966
Submission received: 19 January 2023 / Revised: 31 January 2023 / Accepted: 31 January 2023 / Published: 2 February 2023
(This article belongs to the Section Applied Physics General)

Abstract

:
The alloys that crystallize in a tetragonal ThMn12-type (space group I4/mmm) structure and are based on Fe and rare earth elements are believed to have a potential to plug the performance gap between ferrite and Nd-based magnets. Nevertheless, the progress is hindered by their poor structural stability, compared with other phases competing during the synthesis process, e.g., Th2Zn17-type. In this work, the enthalpies of the formation (and other thermodynamic parameters) of various phases in (Zr, Nd, Ce)-Fe-Si systems were calculated, with paramount focus on the Fe-rich compositions. We compared and discussed the stability range and stabilization routes for amorphous phases, solid solutions, and intermetallics. The beneficial influence of Zr and Si on the crystallization of intermetallic compounds was confirmed, simultaneously being valid for other phases. Among all of the analyzed Fe-rich phases, the lowest values for enthalpy of the formation of the amorphous phase and solid solution were determined for ZrFe10Si2 (−17.5 and −18.2 kJ/mol, respectively). Moreover, substitution by elements with a large atomic radius is indicated as a method for the introduction of topological disorder, giving possibility for the synthesis of metastable phases (even amorphous) and the utilization of more sophisticated synthesis routes in the future.

1. Introduction

High performance permanent magnets are currently mostly based on rare earth elements (REE). One can easily refer to the Sm- or Nd-based compounds [1,2,3,4,5] with the highest energy product BH|max|. The main shortcomings of the Nd2Fe14B (space group P42/mnm) compound is its low Curie temperature (312 °C [6]) and strong deterioration of its anisotropy field with temperature. Therefore, to enhance the anisotropy field and coercivity of modern 2:14:1 magnets, substantial amounts of Dy and other critical REE are used [7,8]. The criticality of such heavy REE forced the development in the field of REE free or at least REE lean permanent magnets. Summarizing, we can distinguish two important branches of research lately, namely: (i) decrease in REE content with the preservation of magnets’ magnetic performance and the (ii) search for novel REE-free materials. When considering the first approach, the tetragonal ThMn12-type phase [9] (space group I4/mmm) is one of the candidates. It is believed to fill the gap between cheap ferrite and expensive 2:14:1 magnets [10]. One of the main issues in the case of materials based on the ThMn12-type phase is their structural stability. Unfortunately, the ThMn12-type phase is less stable for the most REE containing Fe-based compounds than, e.g., hexagonal 2:17-type phase (space group P63/mmc) [11]. Therefore, additional stabilizing elements are necessary to change the energy landscape and to increase the probability of 1:12-type phase formation. There are various types of substitution/doping mechanisms decreasing the relative enthalpy of the formation of 1:12-type phase and simultaneously allowing for its synthesis [12]. Silicon was mentioned as a stabilizer many years ago by Buschow [13] and de Mooij [14], and its role was later confirmed by other authors [15,16,17,18,19]. The formation of the 1:12-type phase was recently confirmed for the ZrFe10Si2 compound [20] and in some other systems, where Zr was substituted with REE [21]. In this paper, we focus on the formation enthalpy calculations of Fe-Si-based alloys with various substitutions (Zr, Ce, and Nd), with an emphasis put on the Fe-rich compositions.
The main aim of the paper was to determine the structural stability of various phases in (Ce, Zr, and Nd)-Fe-Si systems. We compared the enthalpy of the formation of hypothetical amorphous phases, solid solutions, and intermetallic compounds of chosen stoichiometry to answer more detailed questions, including the following: (i) Is it possible to synthesize metallic glass (to have an amorphous precursor for further crystallization of 1:12-type phase) with particular compositions and how broad is the glass forming range? (ii) Is it possible to synthesize NdFe12 or NdFe10Si2 (not synthesized up to date) through the formation of off-stoichiometric alloys (to prepare structurally disordered precursors for further treatment or to restrain the formation of competing 2:17-type phase). To determine the enthalpies of formation, mismatch entropy, and glass forming ability parameter, we utilized a thermodynamic approach based on the semi-empirical Miedema’s model [22,23,24,25], which has been proven to be reliable in these types of calculations and matchless when dealing with the use of computing power. It should also be noted here that other thermodynamical models such as the cBΩ model [26,27,28] have also been applied to magnetic materials [29]. Later, a regular solution model was proposed by Takeuchi and Inoue [30]. To take into account the effect of the additive element on the properties of the binary alloy, asymmetric methods have been proposed [31], which in turn, in many cases, generate inaccuracy by the experience-based choice of asymmetric constituents. Therefore, the asymmetric geometric model proposed by Y. Ouyang et al. [32], independent on the choice of asymmetric element, was successfully implemented before [25,33,34,35] and is utilized in this paper.

2. Methods

The thermodynamic properties are calculated on the basis of Miedema’s semi-empirical model [23,24], also called the macroscopic-atom model. The atoms are taken in simplified form, as building blocks with the boundaries of the Wigner–Seitz cell. The enthalpy is the main thermodynamic quantity that is discussed throughout the paper and its variation is triggered by the contact of two blocks of atoms of different elements. Such an alloying process is thermodynamically described by three quantities that differentiate each element from the periodic table: (i) molar volume (V), (ii) potential similar to work function of an electron (φ), and (iii) density at the boundary of the Wigner–Seitz cell (nws). The interfacial enthalpy ΔHinter(A in B) for solving one mole of A atoms in an excess of B is defined as follows:
H i n t e r ( A   i n   B ) = V A 2 / 3 1 2 ( 1 n w s A 1 / 3 + 1 n w s B 1 / 3 ) { P ( φ ) 2 + Q ( n w s 1 / 3 ) 2 } ,
where P and Q are empirical constants that depend on the alloying elements. To determine the possibility of the formation of various phases, one should compare their formation enthalpies. The formation enthalpy of the amorphous phase is expressed as the following:
Δ H a m = Δ H c h e m + Δ H t o p o ,
where ΔHchem is the chemical contribution equal to the following:
H c h e m = c A c B ( c B s H i n t e r ( A   i n   B ) + c A s H i n t e r ( B   i n   A ) ) ,
while cA, cB and csA(= c A V A 2 3 c A V A 2 3 + c B V B 2 3 ), csB are the volume and surface fractions of each atom type. ΔHtopo is a topological contribution reflecting the disorder in amorphous state and strongly dependent on the melting temperature.
Δ H t o p o = 3.5 ( c A T m , A + c B T m , B ) 10 3 .
The formation enthalpy of the solid solution is defined as follows:
Δ H s s = Δ H c h e m + Δ H e l a s t + Δ H s t r u c t .
The first of the contributions unique for a solid solution is the elastic misfit enthalpy Equation (6), connected with the atom size mismatch:
H e l a s t ( A   i n   B ) = 2 K A G B ( W A W B ) 2 4 G B W A + 3 K A W B ,
where:
W A = V A + α ( φ A φ B ) n w s A ,   and   α = 1.5 V A 2 3 1 n w s A 1 3 + 1 n w s B 1 3 .
K is a bulk modulus, G is a shear modulus, and WA Equation (7) and WB are the corrected molar volumes (because of electron transfer from B to A and from A to B). The structural contribution ΔHstruct originates from the valence and the crystal structure of the solvent and the solute atoms. As a structural term can be only roughly estimated, according to Bakker [24], even for binary transition metal systems, we omitted it from the analysis. Moreover, it should be underlined that the structural contribution for the analyzed transition metal and rare earth element substitutions should not provide much energy difference, according to the structural enthalpies of binaries [24], especially when focussing on Fe-rich compositions. The structural enthalpy for solving one mole of Zr (four valence electrons) in Fe (eight valence electrons) is equal to 0 kJ/mol, while for solving Ce or Nd (both ~ 3 valence electrons) in Fe it is equal to 2 kJ/mol. The emphasis was put on the calculations of enthalpies of formation of amorphous phases and solid solutions, complementary to the vast number of literature data on intermetallic compounds, so as to compare the formation enthalpies of all phases that compete during the solidification process. The thermodynamic quantities determined for three sub-binary systems must be extrapolated to obtain the values for each ternary system. A geometric model was chosen for this purpose, as discussed in the introduction section. Apart from enthalpy of formation, the effect of atomic radius differences was determined by mismatch entropy (Sσ/kB) calculations, on the basis of the relation given by Mansoori et al. [36], while the glass forming ability was additionally estimated using the glass forming ability parameter ΔPHS [22], defined as the product of chemical enthalpy and normalized mismatch entropy ( H c h e m S σ k B ). A more detailed description of the utilized method can be found in [25]. Basic parameters used for the calculations are listed in Table 1.

3. Results and Discussion

The enthalpies of formation of various phases in (Ce, Zr, Nd)-Fe-Si systems were calculated. The emphasis in the discussion is put on the compositions with a high Fe content, due to the arguments put forward in the introduction. This computational paper is intended to provide a keystone for further experimental works on the stabilization of 1:12-type hard magnetic phase. Substituting elements were chosen on the basis of a vast number of references with the intention to improve the structural stability of the 1:12-type phase in comparison with the other competing phases, bearing in mind that the superior aim is to retain the hard magnetic properties. We also check the possibility of the formation of amorphous phase in order to assess the potential of (Ce, Zr, and Nd)-Fe-Si alloys for a more sophisticated synthesis route i.e., vitrification by rapid quenching and further nanocrystallization. It can be beneficial when designing soft/hard magnetic exchange spring magnet composites, comprising of a soft magnetic Fe matrix and hard magnetic nanocrystals [3].
The contour maps of the enthalpies of formation of amorphous phase in Zr-Fe-Si, Ce-Fe-Si, and Nd-Fe-Si (Figure 1) are qualitatively similar. The most negative values are reached for Fe-free compositions. Large negative interfacial enthalpies between (Zr, Ce, and Nd) and Si are responsible for such behavior. Bearing in mind that the large atomic radii difference of alloy constituents promotes glassy state formation, the significant difference in (Zr, Ce, and Nd) and Si atomic radii, which are equal to 1.6–1.82 Å and 1.17 Å [37], respectively, is also fundamental here. The Fe atomic radius with its intermediate value of 1.28 Å plays a minor role. This gives us a first impression of the limited or even destructive role of Fe in the process of the formation of structural disorder in the Fe-rich alloys. The enthalpies of the formation of the amorphous phase in the Zr-Fe-Si system are more negative than those for Ce- and Nd-containing alloys. The most negative values of H a m , obtained for Zr-Si, Ce-Si, and Nd-Si, are equal to −72.1 kJ/mol, −56.2 kJ/mol, and −53.0 kJ/mol for Zr47Si53, Ce43Si57, and Nd43Si57, respectively. A similar role for Zr has been already shown for (Mn, Zr)-Co-Ge alloys [38], where the enthalpies of formation of Zr-containing amorphous phases are much more negative, mainly due to the strongly negative values of the interfacial enthalpies of Zr-Co and Zr-Ge binaries. As can be seen in Figure 1, that a simultaneous decrease in Zr and Si content leads to a gradual increase in H a m . Positive values for H a m are reached for a small number of, mainly Fe-rich, Zr-Fe-Si compositions. For Ce-Fe-Si and Nd-Fe-Si systems, this range is extended to all Ce- and Nd-rich alloys (most compositions with Ce and Nd content extending 95 at.%). Focusing on compositions of interest, namely (Zr, Ce, Nd)Fe12−xSix (0 ≤ x ≤ 2) and (Zr, Ce, Nd)2Fe17-xSix (0 ≤ x ≤ 2), crucial differences are visible again when comparing Zr-containing and Ce- and Nd-containing alloys. For nominal compositions of the ThMn12-type phase, so ZrFe12, CeFe12, and NdFe12, we found H a m = −2.0 kJ/mol, 7.2 kJ/mol, and 6.5 kJ/mol, respectively. For Si-substituted compositions, which are believed to have a lower enthalpy of formation of the ThMn12-type phase, namely ZrFe10Si2, CeFe10Si2, and NdFe10Si2, we found H a m = −17.5 kJ/mol, −8.9 kJ/mol, and −9.3 kJ/mol, respectively. At first, CeFe12 and NdFe12 compositions will unlikely form an amorphous phase, as their enthalpies of formation are highly positive, whereas the enthalpy of formation of the amorphous ZrFe12 phase is slightly negative and theoretically exhibits some driving force for vitrification. Substitution of Fe by Si strongly influences enthalpy values, which are negative also for Ce and Nd. The most negative value of H a m among the analyzed compositions was found for the ZrFe10Si2 alloy, confirming the beneficial influence of the Zr and Si substitutions. One should bear in mind that kinetic parameters were not taken into account in the calculations and the formation of the amorphous phase would be more probable with specific methods (e.g., for high cooling rates). Moreover, the determined H a m values did not include any information on other phases (intermetallics and solid solution) that were competing during solidification process. To have a full picture, one needs to compare all of them. Therefore, another paragraph is devoted to the formation enthalpy of solid solutions.
The contour maps of the enthalpies of formation of solid solutions in Zr-Fe-Si, Ce-Fe-Si, and Nd-Fe-Si (Figure 2) are similar to the respective contour maps calculated for the amorphous phases (presented in Figure 1). The most negative values were reached for low Fe-content and for Zr-containing alloys. When we take a closer look at (Zr, Ce, Nd)Fe12−xSix (0 ≤ x ≤ 2) and (Zr, Ce, Nd)2Fe17-xSix (0 ≤ x ≤ 2) compositions, rare earth alloys exhibited higher values of H s s (as already reported also for H a m ). The enthalpy of formation of solid solution H s s was equal to −2.8 kJ/mol, 5.9 kJ/mol, and 5.9 kJ/mol, for ZrFe12, CeFe12, and NdFe12, respectively. Si-substituted compositions exhibited lower values of H s s than Si-free counterparts mentioned above. H s s is equal to −18.2 kJ/mol, −10.3 kJ/mol, and -10.0 kJ/mol, for ZrFe10Si2, CeFe10Si2, and NdFe10Si2, respectively. All of the H a m and H s s values listed above are gathered in Table 2 along with the data for 2:17-type phases (believed to be the main competing phases). The H a m and H s s values calculated for Zr2Fe17 were more negative than those determined for the ZrFe12 phase. Inverse dependence was observed for Ce and Nd substitutions, where both enthalpies were higher for 2:17-type phases. Si substitution plays a crucial role in all cases, reducing the enthalpy values significantly. To have a clear comparison of the H a m and H s s , the difference between both values (ΔHam-ss) is presented in Figure 3 in the form of a contour map. This gives us a more reliable way to show the phase formation preferences for each composition. The ΔHam-ss data specific for 1:12 and 2:17 compositions are shown in Table 2. The positive values point out compositions for which the formation of the solid solution is thermodynamically preferred over an amorphous phase. This is the case for all 1:12-type compositions gathered in Table 2. In turn, for all 2:17-type phases, except Ce2Fe17, the amorphous phase is preferred. If we take into account the influence of the structural enthalpy term (exact values discussed in the Calculations section), this preference should even be slightly strengthened. The most negative value of ΔHam-ss = −1.3 kJ/mol was calculated for the Nd2Fe17 composition. To have a more overarching view, one should focus on the general trends visible in Figure 3. To improve the glass forming ability, or to at least destabilize the solid solution and to introduce structural disorder, the content of the elements with a large atomic radius, and preferably with a large negative heat of mixing with Fe (e.g., Zr), should be further increased at the expense of the Fe and Si content. The change in ΔHam-ss sign (for about 15 at.% of Zr, Ce, and Nd) marks the change in phase preference (Figure 3).
To extend the analysis, the composition dependence of the normalized mismatch entropy was also calculated (Figure S1). Most of the compositions exhibited Sσ/kB values higher than 0.1. This is the minimum value for the achievement of considerable glass forming ability [30]. Normalized mismatch entropy reflects the effect of atomic radii mismatch between the main constituents and, as expected due to the large atomic radius difference between Zr, Ce, Nd, and Si, reaches the highest values for Fe-free compositions. The lowest value of Sσ/kB among the alloys of interest (Table 2) was determined for ZrFe12 and was equal to 0.11. The highest mismatch reflected in the Sσ/kB value of 0.44 was reached for Ce2Fe15Si2 and Nd2Fe15Si2. As expected, the results indicate a beneficial effect of the presence of largely mismatched atomic radii elements on the possibility of the formation of a disordered structure.
To combine both, the enthalpy of formation and the mismatch entropy, the glass forming ability parameter ∆PHS (Figure S2) was also calculated as a product of Sσ/kB and chemical enthalpy (basic term in the expression for enthalpy of formation for both, amorphous phase and solid solution). The obtained results were comparable to those of the formation enthalpy of the amorphous phase. The most negative values indicating the highest glass forming ability were reached again for the Fe-free compositions. Large negative interfacial enthalpies and an atomic radii difference between (Zr, Ce, and Nd) and Si are believed to be crucial here, as already stated in the case of H a m . However, ∆PHS reached the most negative values for the 1:12-type phases, −5.8 kJ/mol for CeFe10Si2 and −6.0 kJ/mol for NdFe10Si2. To conclude the first part of our research, one should expect that the substitution of Fe simultaneously by the elements with larger and smaller atomic radii, as in the case of (Zr, Ce, Nd) and Si substitutions, will promote the vitrification process. The minimum content of Zr, Ce, and Nd (equal to 15 at.%), determined on the basis of H a m s s , should be taken as an estimation only, as kinetic factors of synthesis process can be decisive. We can be more precise when comparing systems with various substitutions for similar compositions, and the data gathered in Table 2 can be treated as a guide. Nevertheless, qualitative and quantitative agreement of the calculated enthalpies with existing data can be found. For example, when one compares the H a m s s contour map of the Nd-Fe-Si system with the experimental data reported in [39] (Nd90-xFexSi10 amorphous alloys were successfully synthesized in the range of 20 ≤ x ≤ 50), it is evident that the very broad glass forming range confirmed here for the low Si content converged well with the literature.
Planning future experimental works, we should compare the calculated values of the enthalpies of formation of amorphous alloys and solid solutions to those of intermetallic compounds apparently also competing during the synthesis process. There is experimental evidence on the formation of the phase of interest (ThMn12-type) in (Zr, Ce)Fe10Si2 and similar systems. Sakurada et al. [40] showed that the ThMn12-type phase crystallized (along with relatively small amount of α-Fe) in Nd1-xZrxFe10Si2 (x = 0.25, 0.5 and 0.75), indicating the atomic radius of the rare earth site element as a controlling factor. Moreover, it has been reported that the changes in the cooling rate in a rapid quenching process influenced the phase constitution. The ThMn12-type structure has been stabilized in the samples quenched with a roll velocity of 10 m/s for Sm, Gd, and Er substitutions, while for 30 m/s, Zr and Nd also favored its formation [40]. Other studies have confirmed the stability of the tetragonal structure in (REE)Fe10Si2 for rare earth site elements having an atomic radius equal or smaller than 0.181 nm [20,41,42]. CeFe10Si2 and NdFe10Si2 preferably crystallize in the Th2Zn17-type structure [21]. The stabilization of the 1:12 structure in such alloys can be obtained through the substitution of larger (REE) atoms with smaller Zr [21,43,44]. To date, there is no direct evidence on the stability of 1:12-type phase in the alloys containing a combination of Nd, Zr, and Ce atoms. Nevertheless, it has been reported by Gabay and Hadjipanayis [18] for (Sm, Zr, and Ce)Fe10Si2 that the combination of three rare earth site substitutions resulted in the stabilization of a 1:12 structure with an improved hard magnetic performance (in comparison with stable alloys with the mixture of two elements on rare earth site) and simultaneously with a low content of critical rare earth elements.
A more exact comparison requires the consideration of the values of the enthalpies of formation of amorphous phases, solid solutions, and intermetallics. Because of the coexistence of many various types of phases, e.g., tetragonal 1:12-type (space group I4/mmm) or hexagonal (space group P63/mmc) and rhombohedral (space group R3m) 2:17-type phases in the region of interest in the phase diagram, and the limitations of Miedema’s model-based semi-empirical approach, the enthalpies of formation of intermetallic compounds were taken from theoretical works based on more sophisticated methods, e.g., ab initio calculations. The general conclusion, which can be drawn from the literature data, is the improvement of 1:12-type phase stability by the substitution of Fe by various stabilizing elements, as reported for SmFe12 and SmFe11(Co, Ti, V, and Ga) [45]. The lowering of the enthalpy of formation for the substituted SmFe11(Co, Ti, V, Ga) samples resembled the variation of enthalpies of formation of solid solutions reported in the present paper (Table 2). While the energy of formation for SmFe12 was approximately equal to 2 kJ/mol, it decreased for Co, Ti, V, and Ga substitutions and reached a value lower than -5 kJ/mol for SmFe11Ti [45]. The recent study by Landa et al. [46] showed that nickel can also play a role as a stabilizer in the 1:12-type structure. Comparing these values with the formation energies reported in the paper by Harashima et al. [11] shows us a general picture about the relative stability of various intermetallic phases. At first, the value of the formation energy of SmFe12 was determined to be equal to approx. 0.2 eV. Recalculation of this value to kJ/mol provided a much higher value than the 2 kJ/mol calculated in [45]. Nevertheless, the overall picture remained consistent. The formation energy of 1:12-type phase increased with increasing the atomic radius on the rare earth site. The energy of formation of ZrFe12 was strongly negative (approx. −0.25 eV). It also remained negative for Lu, Tm, Er, Ho, and Dy. A further increase in the atomic radius made the formation of the 1:12-type phase unfavorable. For NdFe12, the formation energy was highly positive and doubled the value determined for SmFe12. When comparing these values with the enthalpies of formation reported in Table 2, it is clear that the energy of formation of intermetallics always had a higher absolute value. It means that in the case of compositions of interest, the formation of intermetallic phases was preferred in respect to the solid solution and amorphous phase. Nevertheless, one has to bear in mind that any deviation from ideal stoichiometry and growing chemical disorder could reverse this situation. Finally, we compared the formation energy values for the competing 2:17-type and 1:12-type phases. As opposed to the enthalpies of formation of solid solutions (Table 2), the energies of formation of rhombohedral Th2Zn17 and hexagonal Th2Ni17 were 0.04–0.16 eV lower than the energy of formation of tetragonal ThMn12 phase for all of the analyzed rare earth site substitutions (Hf, Zr, Sc, Lu, Tm, Er, Ho, Dy, Y, Gd, Sm, Nd, Pr, and La) [11]. In conclusion, the utilization of Zr as the main rare earth site substitution had the highest impact on the lowering of the enthalpy of formation of all types of phases, intermetallics (as already reported in [11]), solid solutions, and amorphous. A further decrease is expected when substituting elements with a smaller atomic radius (e.g., Si, as discussed in the present paper) on Fe sites. The situation was slightly different when comparing various stable intermetallic phases containing the same chemical elements. In this case, the crystalline structure and hence varying interatomic distances played a more significant role. Therefore, the discussion that excludes mentioned factors is rather incomplete, but some general trends are still valid. The complex Zr-Fe-Si system with a vast number of stable intermetallic phases can be treated as an example [47]. Intermetallics with a relatively low content of Si and Zr (in comparison to Fe) are characterized by enthalpy of formation values equal to approximately −52 kJ/mol and −56 kJ/mol for Fe16Si7Zr6 and Fe3SiZr2, respectively [48]. A decreasing content of Fe (Fe4Si7Zr4) leads to a more negative formation enthalpy equal to approx. −68 kJ/mol [48]. We observed a similar trend for the formation enthalpy of a solid solution in the Zr-Fe-Si system (Figure 2), with slightly less negative values of enthalpy compared to the intermetallics. The enthalpy changes from −28 kJ/mol for Fe3SiZr2, through −31 kJ/mol for Fe16Si7Zr6, to −48 kJ/mol for the Fe4Si7Zr4 composition. In general, intermetallic phases are favored for stoichiometric (1:12 and 2:17) compositions, but it is believed that the further introduction of chemical disorder on rare earth and Fe sites can result in the formation of a solid solution (as in the case of high entropy alloys) or even amorphous phase. In this instance, the energy convex hull [12] can be used to determine the relative stability of the phases. The formation of a glassy structure can be also favored by the utilization of non-equilibrium synthesis methods.

4. Conclusions

In the present paper, we focused on the calculations of formation enthalpy of (Zr, Ce, and Nd)-Fe-Si systems. The most effort was put into the comparison of the formation enthalpies of various phases and the characterization of the mutual dependences between them for the Fe-rich compositions. As it was already elucidated, the ThMn12-type phase was, for most Fe-based compounds containing rare earth elements, less stable than other intermetallic compounds with a similar stoichiometry [11]. The formation of the 1:12-type phase has been lately confirmed, e.g., in ZrFe10Si2 compound [12], which was a starting point for our research. The beneficial influence of Zr was confirmed for intermetallic compounds, but was also apparent when comparing the enthalpy of formation of solid solution H s s , which was equal to −2.8 kJ/mol, 5.9 kJ/mol, and 5.9 kJ/mol for ZrFe12, CeFe12, and NdFe12, respectively. Further lowering of the H s s was observed after the substitution of Fe by Si; ZrFe10Si2, CeFe10Si2, and NdFe10Si2 exhibited lower values of H s s (−18.2 kJ/mol, −10.3 kJ/mol, and −10.0 kJ/mol, respectively) than their Si-free counterparts. The substitution of Nd and Ce had a similar impact on the formation enthalpies for all of the analyzed phases. Both elements were not favorable when aiming for the improvement of stability, but could be substituted to some extent (needed from the point of view of optimization of magnetic properties). The most negative value for H a m among the analysed compositions was reached for the ZrFe10Si2 alloy, confirming advantageous impact of Zr and Si atoms. In turn, the glass forming ability parameter was the lowest for the compositions containing Nd and Ce, indicating the great importance of large atomic radii substitutions on the formation of topological disorder. Intermetallic phases were favored for the analyzed stoichiometric compositions (1:12 and 2:17), but the substituting effect could be crucial for the formation of chemical and topological disorder and for the stabilization of solid solutions (growing entropy) or amorphous phase.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/app13031966/s1, Figure S1: Compositional dependences of normalized mismatch entropy Sσ/kB of ternary (Zr, Ce, Nd)-Fe-Si systems; Figure S2: Compositional dependences of glass forming ability parameter ΔPHS of ternary (Zr, Ce, Nd)-Fe-Si systems.

Author Contributions

Conceptualization, M.K. and Z.Ś.; data curation, M.K.; formal analysis, M.K. and Z.Ś; methodology, Z.Ś.; supervision, Z.Ś.; validation, M.K.; visualization, M.K. and Z.Ś.; writing—original draft, M.K., writing—review and editing, M.K. and Z.Ś. All authors have read and agreed to the published version of the manuscript.

Funding

M.K. was financially supported by the project “Środowiskowe interdyscyplinarne studia doktoranckie w zakresie nanotechnologii” (“Environmental interdisciplinary doctoral studies in nanotechnology”) No. POWR.03.02.00-00-I032/16 under the European Social Fund–Operational Programme Knowledge Education Development, Axis III Higher Education for Economy and Development, Action 3.2 PhD Programme.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Acknowledgments

M.K. was financially supported by the project “Środowiskowe interdyscyplinarne studia doktoranckie w zakresie nanotechnologii” (“Environmental interdisciplinary doctoral studies in nanotechnology”) No. POWR.03.02.00-00-I032/16 under the European Social Fund–Operational Programme Knowledge Education Development, Axis III Higher Education for Economy and Development, Action 3.2 PhD Programme.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Zhou, X.; Song, X.; Jia, W.; Xiao, A.; Yuan, T.; Ma, T. Identifications of SmCo5 and Smn+1Co5n−1-type phases in 2:17-type Sm-Co-Fe-Cu-Zr permanent magnets. Scr. Mater. 2020, 182, 1–5. [Google Scholar] [CrossRef]
  2. Tang, X.; Lia, J.; Sepehri-Amin, H.; Ohkubo, T.; Hioki, K.; Hattori, A.; Hono, K. Improved coercivity and squareness in bulk hot-deformed Nd–Fe–B magnets by two-step eutectic grain boundary diffusion process. Acta Mater. 2021, 203, 116479. [Google Scholar] [CrossRef]
  3. Abuin, M.; Turgut, Z.; Aronhime, N.; Keylin, V.; Leary, A.; Degeorge, V.; Horwath, J.; Semiatin, S.L.; Laughlin, D.E.; McHenry, M.E. Determination of pressure effects on the α→γ phase transition and size of Fe in Nd-Fe-B spring exchange magnets. Metall. Mater. Trans. A 2015, 46, 5002–5010. [Google Scholar] [CrossRef]
  4. Li, X.; Yao, Q.; Lu, Z.; Ma, L.; Du, Y.; Wang, J.; Zhou, H.; Rao, G. Atomic Mobilities in Liquid and fcc Nd-Fe-B systems and their application in the design of quenching Nd2Fe14B alloys. Metall. Mater. Trans. A 2021, 52, 2948–2958. [Google Scholar] [CrossRef]
  5. Haider, S.K.; Ngo, H.M.; Kim, D.; Kang, Y.S. Enhancement of anisotropy energy of SmCo5 by ceasing the coupling at 2c sites in the crystal lattice with Cu substitution. Sci. Rep. 2021, 11, 10063. [Google Scholar] [CrossRef]
  6. Hirayama, Y.; Miyake, T.; Hono, K. Rare-Earth lean hard magnet compound NdFe12N. JOM 2015, 67, 1344–1349. [Google Scholar] [CrossRef]
  7. Liu, X.B.; Altounian, Z. The partitioning of Dy and Tb in NdFeB magnets: A first-principles study. J. Appl. Phys. 2012, 111, 07A701. [Google Scholar] [CrossRef]
  8. Zhou, T.; Qu, P.; Pan, W.; Liu, R.; Li, M.; Ur Rehman, S.; Zhong, Z.; Xie, G. Sintered NdFeB magnets with Tb-Dy double-layer core/shell structure were fabricated by double alloy method and grain boundary diffusion. J. Alloys Compd. 2021, 856, 158191. [Google Scholar] [CrossRef]
  9. Solzi, M.; Pareti, L.; Moze, O.; David, W.I.F. Magnetic anisotropy and crystal structure of intermetallic compounds of the ThMn12 structure. J. Appl. Phys. 1988, 64, 5084–5087. [Google Scholar] [CrossRef]
  10. Coey, J.M.D. Permanent magnets: Plugging the gap. Scr. Mater. 2012, 67, 524–529. [Google Scholar] [CrossRef]
  11. Harashima, Y.; Fukazawa, T.; Kino, H.; Miyake, T. Effect of R-site substitution and the pressure on stability of RFe12: A first-principles study. J. Appl. Phys. 2018, 124, 163902. [Google Scholar] [CrossRef]
  12. Fukazawa, T.; Harashima, Y.; Miyake, T. First-principles study on the stability of (R, Zr)(Fe, Co, Ti)12 against 2-17 and unary phases (R=Y, Nd, Sm). Phys. Rev. Mater. 2022, 6, 054404. [Google Scholar] [CrossRef]
  13. Buschow, K.H.J. Structure and properties of some novel ternary Fe-rich rare-earth intermetallics (invited). J. Appl. Phys. 1988, 64, 3130–3135. [Google Scholar] [CrossRef]
  14. De Mooij, D.B.; Buschow, K.H.J. Some novel ternary ThMn12-type compounds. J. Less-Common. Met. 1988, 136, 207–215. [Google Scholar] [CrossRef]
  15. Solzi, M.; Xue, R.H.; Pareti, L. Magnetic anisotropy and first-order magnetization processes in Sm(Fe1−xCox)10M2 (M = Ti, Si) compounds. J. Magn. Magn. Mater. 1990, 88, 44–50. [Google Scholar] [CrossRef]
  16. Yücel, A.; Elmali, A.; Çakır, R.; Elerman, Y. Competing magnetic interactions in Dy1−xSmxFe10Si2 compounds. Solid State Commun. 2005, 135, 745–748. [Google Scholar] [CrossRef]
  17. Zhou, C.; Pinkerton, F.E.; Herbst, J.F. High Curie temperature of Co-Fe-Si compound with ThMn12 structure. Scr. Mater. 2015, 95, 66–69. [Google Scholar] [CrossRef]
  18. Gabay, A.M.; Hadjipanayis, G.C. Mechanochemical synthesis of magnetically hard anisotropic RFe10Si2 powders with R representing combinations of Sm, Ce and Zr. J. Magn. Magn. Mater. 2017, 422, 43–48. [Google Scholar] [CrossRef]
  19. Gjoka, M.; Sarafidis, C.; Psycharis, V.; Devlin, E.; Niarchos, D.; Hadjipanayis, G. Structure and magnetic properties of Sm1−xZrxFe10Si2 (x = 0.9–0.6) alloys. IOP Conf. Ser. J. Phys. Conf. Ser. 2017, 903, 12033. [Google Scholar] [CrossRef]
  20. Gabay, A.; Hadjipanayis, G. ThMn12-type structure and uniaxial magnetic anisotropy in ZrFe10Si2 and Zr1-xCexFe10Si2 alloys. J. Alloys Compd. 2016, 657, 133–137. [Google Scholar] [CrossRef] [Green Version]
  21. Barandiaran, J.M.; Martin-Cid, A.; Schonhobel, A.M.; Garitaonandia, J.S.; Gjoka, M.; Niarchos, D.; Makridis, S.; Pasko, A.; Aubert, A.; Mazaleyrat, F.; et al. Nitrogenation and sintering of (Nd-Zr)Fe10Si2 tetragonal compounds for permanent magnets applications. J. Alloys Compd. 2019, 784, 996–1002. [Google Scholar] [CrossRef]
  22. Bhatt, J.; Jiang, W.; Junhai, X.; Qing, W.; Dong, C.; Murty, B.S. Optimization of bulk metallic glass forming compositions in Zr–Cu–Al system by thermodynamic modeling. Intermetallics 2007, 15, 716–721. [Google Scholar] [CrossRef]
  23. De Boer, F.R.; Boom, R.; Mattens, W.C.M.; Miedema, A.R.; Niessen, A.K. Cohesion in Metals; De Boer, F.R., Pettifor, D.G., Eds.; North Holland: Amsterdam, the Netherlands, 1988; Volume 1. [Google Scholar]
  24. Bakker, H. Enthalpies in Alloys: Miedema’s Semi-Empirical Model; Materials Science Foundations 1; Trans Tech Publications Ltd.: Zurich, Switzerland, 1998. [Google Scholar]
  25. Śniadecki, Z.; Narojczyk, J.W.; Idzikowski, B. Calculation of glass forming ranges in the ternary Y-Cu-Al system and its sub-binaries based on geometric and Miedema’s models. Intermetallics 2012, 26, 72–77. [Google Scholar] [CrossRef]
  26. Varotsos, P.; Alexopoulos, K. Estimation of the migration enthalpy and entropy for cation vacancy motion in alkali halides with the NaCl-type structure. Phys. Rev. B 1977, 15, 2348–2351. [Google Scholar] [CrossRef]
  27. Varotsos, P.; Alexopoulos, K. Calculation of the formation entropy of vacancies due to anharmonic effects. Phys. Rev. B 1977, 15, 4111–4114. [Google Scholar] [CrossRef]
  28. Varotsos, P.; Alexopoulos, K. Negative activation volumes of defects in solids. Phys. Rev. B 1980, 21, 4898–4899. [Google Scholar] [CrossRef]
  29. Varotsos, C.; Lazaridou, M.; Alexopoulos, K.; Varotsos, P. Comments on “The Temperature and Pressure Dependence of Disaccommodation in a Manganese Zinc Ferrite Single Crystal”. Jpn. J. Appl. Phys. 1985, 24, 781. [Google Scholar] [CrossRef]
  30. Takeuchi, a.; Inoue, A. Calculations of mixing enthalpy and mismatch entropy for ternary amorphous alloys. Trans. Japan Instig. Met. 2000, 41, 1372–1378. [Google Scholar] [CrossRef]
  31. Toop, G.W. Predicting ternary activities using binary data. Trans. Metall. Soc. AIME 1965, 233, 850–855. [Google Scholar]
  32. Ouyang, Y.; Zhong, X.; Du, Y.; Jin, Z.; He, Y.; Yuan, Z. Formation enthalpies of Fe–Al–RE ternary alloys calculated with a geometric model and Miedema’s theory. J. Alloys Compd. 2006, 416, 148–154. [Google Scholar] [CrossRef]
  33. Śniadecki, Z. Glass-forming ability of Fe-Ni alloys substituted by group V and VI transition metals (V, Nb, Cr, Mo) studied by thermodynamic modeling. Metall. Mater. Trans. A 2020, 51, 4777–4785. [Google Scholar] [CrossRef]
  34. Śniadecki, Z. The influence of 3d and 4d transition metals on the glass forming ability of ternary FeCo-based alloys. Metall. Mater. Trans. A 2021, 52, 1861–1868. [Google Scholar] [CrossRef]
  35. Kołodziej, M.; Śniadecki, Z. The formation of structural disorder in FeNi-based alloys—Theoretical approach. Mater. Lett. 2022, 326, 132917. [Google Scholar] [CrossRef]
  36. Mansoori, G.A.; Carnahan, N.F.; Starling, K.E.; Leland, T.W., Jr. Equilibrium thermodynamic properties of the mixture of hard spheres. J. Chem. Phys. 1971, 54, 1523–1525. [Google Scholar] [CrossRef]
  37. Brandes, E.A.; Brook, G. (Eds.) Smithells Metals Reference Book, 7th ed.; Butterworth-Heinemann Ltd.: Oxford, UK, 1992. [Google Scholar]
  38. Gębara, P.; Śniadecki, Z. Structure, magnetocaloric properties and thermodynamic modeling of enthalpies of formation of (Mn,X)-Co-Ge (X = Zr, Pd) alloys. J. Alloys Compd. 2019, 796, 153–159. [Google Scholar] [CrossRef]
  39. Zhang, W.; Takeuchi, A.; Inoue, A. Amorphous Nd-Fe-Si thick ribbons and their hard magnetic properties. Mater. Trans. JIM 1997, 38, 1027–1030. [Google Scholar] [CrossRef]
  40. Sakurada, S.; Tsutai, A.; Sahashi, M. A study on the formation of ThMn12 and NaZn13 structures in RFe10Si2. J. Alloys Compd. 1992, 187, 67–71. [Google Scholar] [CrossRef]
  41. Stefański, P.; Wrzeciono, A. Structural and magnetic properties of RFe10Si2 compounds. J. Magn. Magn. Mater. 1989, 82, 125–128. [Google Scholar] [CrossRef]
  42. Ohashi, K.; Tawara, Y.; Osugi, R.; Shimao, M. Magnetic properties of Fe-rich rare-earth intermetallic compounds with a ThMn12 structure. J. Appl. Phys. 1988, 64, 5714–5716. [Google Scholar] [CrossRef]
  43. Gjoka, M.; Psycharis, V.; Devlin, E.; Niarchos, D.; Hadjipanayis, G. Effect of Zr substitution on the structural and magnetic properties of the series Nd1−xZrxFe10Si2 with the ThMn12 type structure. J. Alloys Compd. 2016, 687, 240–245. [Google Scholar] [CrossRef]
  44. Salazar, D.; Martín-Cid, A.; Garitaonandia, J.S.; Hansen, T.C.; Barandiaran, J.M.; Hadjipanayis, G.C. Role of Ce substitution in the magneto-crystalline anisotropy of tetragonal ZrFe10Si2. J. Alloys Compd. 2018, 766, 291–296. [Google Scholar] [CrossRef]
  45. Dirba, I.; Harashima, Y.; Sepehri-Amin, H.; Ohkubo, T.; Miyake, T.; Hirosawa, S.; Hono, K. Thermal decomposition of ThMn12-type phase and its optimum stabilizing elements in SmFe12-based alloys. J. Alloys Compd. 2020, 813, 152224. [Google Scholar] [CrossRef]
  46. Landa, A.; Söderlind, P.; Moore, E.E.; Perron, A. Thermodynamics and Magnetism of SmFe12 Compound Doped with Co and Ni: An Ab Initio Study. Appl. Sci. 2022, 12, 4860. [Google Scholar] [CrossRef]
  47. Lisenko, L.A.; Ban, Z.; Gladisevskii, E.I. Investigation of the system Zr-Fe-Si. Croat. Chem. Acta 1971, 43, 113–118. [Google Scholar]
  48. Cui, J.; Shen, Y.; Liu, X. Experimental investigation and thermodynamic modeling of the Fe−Si−Zr system. Calphad 2019, 65, 385–401. [Google Scholar] [CrossRef]
Figure 1. Compositional dependences of formation enthalpies of amorphous phase H a m of ternary (Zr, Ce, Nd)-Fe-Si systems. Key composition regions are marked in grey (2:17-type phases) and black (1:12 type phases). Black solid lines are ∆H = 0 indicators.
Figure 1. Compositional dependences of formation enthalpies of amorphous phase H a m of ternary (Zr, Ce, Nd)-Fe-Si systems. Key composition regions are marked in grey (2:17-type phases) and black (1:12 type phases). Black solid lines are ∆H = 0 indicators.
Applsci 13 01966 g001
Figure 2. Compositional dependences of the formation enthalpies of the solid solution H s s of ternary (Zr, Ce, and Nd)-Fe-Si systems. Key composition regions are marked in grey (2:17-type phases) and black (1:12-type phases). Black solid lines are ∆H = 0 indicators.
Figure 2. Compositional dependences of the formation enthalpies of the solid solution H s s of ternary (Zr, Ce, and Nd)-Fe-Si systems. Key composition regions are marked in grey (2:17-type phases) and black (1:12-type phases). Black solid lines are ∆H = 0 indicators.
Applsci 13 01966 g002
Figure 3. Compositional dependences of difference of two quantities, of formation enthalpies of amorphous phase and solid solution H a m s s of ternary (Zr, Ce, and Nd)-Fe-Si systems. Key composition regions are marked in grey (2:17-type phases) and black (1:12 type phases). Black solid lines are ∆H = 0 kJ/mol indicators (negative H a m s s values indicate a lower enthalpy of formation of amorphous phase).
Figure 3. Compositional dependences of difference of two quantities, of formation enthalpies of amorphous phase and solid solution H a m s s of ternary (Zr, Ce, and Nd)-Fe-Si systems. Key composition regions are marked in grey (2:17-type phases) and black (1:12 type phases). Black solid lines are ∆H = 0 kJ/mol indicators (negative H a m s s values indicate a lower enthalpy of formation of amorphous phase).
Applsci 13 01966 g003
Table 1. Basic parameters used in the calculations of enthalpies of formation [24,37].
Table 1. Basic parameters used in the calculations of enthalpies of formation [24,37].
Chemical
Element
φ
(V)
nws
(Arb. Units)
V
(10−6 m3/mol)
K
(1010 Pa)
G
(1010 Pa)
d
(Å)
Ce3.181.6921.622.3951.1973.64
Fe4.935.557.0916.838.1522.56
Nd3.191.7320.583.2681.4523.64
Si4.703.388.609.8883.9732.34
Zr3.452.8014.008.3353.4143.20
Table 2. Enthalpies of the formation of amorphous phase H a m , of solid solution H s s , their difference H a m s s , normalized mismatch entropy Sσ/kB, and glass forming ability parameter ∆PHS calculated for the chosen binary and ternary 1:12-type and 2:17-type compositions.
Table 2. Enthalpies of the formation of amorphous phase H a m , of solid solution H s s , their difference H a m s s , normalized mismatch entropy Sσ/kB, and glass forming ability parameter ∆PHS calculated for the chosen binary and ternary 1:12-type and 2:17-type compositions.
ΔHam (kJ/mol)ΔHss (kJ/mol)ΔHam-ss (kJ/mol)Sσ/kBΔPHS (kJ/mol)
ZrFe12−2.0−2.80.80.11−0.9
CeFe127.25.91.30.310.3
NdFe126.55.90.60.310.1
ZrFe10Si2−17.5−18.20.70.14−3.5
CeFe10Si2−8.9−10.31.40.36−5.8
NdFe10Si2−9.3−10.00.70.36−6.0
Zr2Fe17−4.6−3.6−1.00.14−1.6
Ce2Fe177.47.80.40.400.6
Nd2Fe176.57.8−1.30.400.1
Zr2Fe15Si2−16.6−15.7−0.90.17−4.0
Ce2Fe15Si2−5.1−5.0−0.10.44−5.3
Nd2Fe15Si2−5.7−4.7−1.00.44−5.6
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kołodziej, M.; Śniadecki, Z. Thermodynamic Modeling of Formation Enthalpies of Amorphous and Crystalline Phases in Zr, Nd, and Ce-Substituted Fe-Si Systems. Appl. Sci. 2023, 13, 1966. https://doi.org/10.3390/app13031966

AMA Style

Kołodziej M, Śniadecki Z. Thermodynamic Modeling of Formation Enthalpies of Amorphous and Crystalline Phases in Zr, Nd, and Ce-Substituted Fe-Si Systems. Applied Sciences. 2023; 13(3):1966. https://doi.org/10.3390/app13031966

Chicago/Turabian Style

Kołodziej, Mieszko, and Zbigniew Śniadecki. 2023. "Thermodynamic Modeling of Formation Enthalpies of Amorphous and Crystalline Phases in Zr, Nd, and Ce-Substituted Fe-Si Systems" Applied Sciences 13, no. 3: 1966. https://doi.org/10.3390/app13031966

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop