Blog

First principles and mean field study on the magnetocaloric effect of YFe3 and HoFe3 compounds | Scientific Reports

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  13, Article number: 2876 (2023 ) Cite this article Neodymium Powder

First principles and mean field study on the magnetocaloric effect of YFe3 and HoFe3 compounds | Scientific Reports

In this work, the magnetothermal characteristics and magnetocaloric effect in YFe3 and HoFe3 compounds are calculated as function of temperature and magnetic field. These properties were investigated using the two-sublattice mean field model and the first-principles DFT calculation using the WIEN2k code. The two-sublattice model of the mean-field theory was used to calculate the temperature and field-dependences of magnetization, magnetic heat capacity, magnetic entropy, and the isothermal change in entropy ∆Sm. We used the WIEN2k code to determine the elastic constants and, subsequently, the bulk and shear moduli, the Debye temperature, and the density-of-states at Ef. According to the Hill prediction, YFe3 has bulk and shear moduli of roughly 99.3 and 101.2 GPa respectively. The Debye temperature is ≈ 500 K, and the average sound speed is ≈ 4167 m/s. In fields up to 60 kOe and at temperatures up to and above the Curie point for both substances, the trapezoidal method was used to determine ∆Sm. For instance, the highest ∆Sm values for YFe3 and HoFe3 in 30 kOe are approximately 0.8 and 0.12 J/mol. K, respectively. For the Y and Ho systems, respectively, the adiabatic temperature change in a 3 T field decreases at a rate of around 1.3 and 0.4 K/T. The ferro (or ferrimagnetic) to paramagnetic phase change in these two compounds, as indicated by the temperature and field dependences of the magnetothermal and magnetocaloric properties, ∆Sm and ∆Tad, is a second-order phase transition. The Arrott plots and the universal curve for YFe3 were also calculated and their features give an additional support to the second order nature of the phase transition.

Sukanya Ghosh, Soheil Ershadrad, … Biplab Sanyal

Satadeep Bhattacharjee & Seung-Cheol Lee

Subhasish Mandal, Kristjan Haule, … David Vanderbilt

The magneto-caloric effect (MCE) is a well-known interesting phenomenon that involves a change in the temperature upon applying/removing an external magnetic field. Its physics and applications are of continuous interest to researchers worldwide. [e.g.,1,2].

Some of the intermetallic compounds that exist in the rare-earth-iron system are, for example, RFe2, RFe3, R6Fe23 and R2Fe173,4,5. Much work has been dedicated to the magnetic properties of those compounds. The RT2 compounds, where T = Fe, Co, have been also investigated6,7,8,9. Roe et al.10 determined the magnetic structure of the iron-holmium binary system. Hoffer and Salmans11 measured the magnetization of the RFe3 compounds as function of temperature. On the application side, magnetic refrigerators, in particular those operating close to room temperature are a subject of extensive research. One of the reasons of such interest in the magnetocaloric-based refrigerators are the relatively higher efficiency and the lower impact on environment as compared to the conventional ones12,13. Rare-earth intermetallic compounds were a subject of many works due to their distinguished electronic, magnetic and magnetocaloric properties14,15,16,17,18. The RFe3 compounds, in particular, are of interest for both their magnetostrictive and magnetic properties19,20,21, however not much investigations were done on their magnetocaloric properties.

Our motivation, in this paper, is studying the magnetothermal and the magnetocaloric effect in YFe3 and HoFe3, using both ab initio calculation and the molecular field theory, in the two-sublattice model22,23,24. We report on both temperature and magnetic field dependences of magnetization, magnetic specific heat magnetic entropy, isothermal change in entropy and adiabatic change in temperature in the temperature interval of 0–800 K. We report on new calculation of the elastic constants, the bulk, the shear moduli and the electronic and lattice contributions to the total heat capacity and entropy using the well-known WIEN2k code package25,26,27,28,29. Another motivation of the present work is studying the order of the phase transition in these compounds in the light of the dependence of their properties on both temperature and field.

To calculate the elastic moduli of YFe3 we used the linear augmented plane-wave (LAPW) method by WIEN2K code30. This Package uses the energy approach calculation 31 that calculates the elastic constants by second-order derivative. The results of these calculations are then used to calculate the Debye temperature and the electronic heat capacity coefficient for \({\text{ YFe}}_{3}\) . These two parameters for HoFe3 are obtained from available published data (Persson, the materials project)32.

The Wien2K electronic structure code is based on the Density Functional Theory (DFT)33. It employs the Full-Potential Linearized Augmented Plane Wave (FPLAPW)34. For correlation and exchange potentials, the Wien2k code30 employs the Local Density Approximation (LDA) of Perdew and Wang35 and the Generalized Gradient Approximation (GGA) of Perdew, Burke, and Ernzerhof36. Core and valence states are calculated in a self-consistent manner. For the spherical part of the potential, the core states are treated fully relativistically, while the valence states are treated using the full potential. To reduce linearization errors in R and Fe spheres, local orbital extensions37 with a converged basis of approximately 1000 basis functions are used. For the self-consistent band structure calculations, we used the modified tetrahedron method for Brillouin zone integration. Self-consistent calculations were performed, and convergence was checked by varying the number of k points in the irreducible Brillouin zone up to 64 points The ab initio calculated elastic properties of different systems have been reported, e.g. for YFe538, PrX2 (X = Fe, Mn, Co) compounds39, Cr-based full-Heusler alloys40 and GdFe241. Also ab initio calculation on the effect of applying pressure on PtXSb and GdN systems were reported, respectively by Habbak et al.42 and Reham Shabara et al.43. Magnetic properties and electronic structure of ThCo4B were reported by Abu-Elmagd et al.44 and Sherif Yehia et al.45 have reported on the spin and charge density maps of SmCo5.

The Lagrangian theory of elasticity is used to describe the elastic properties. According to this theory, a solid is a homogeneous and anisotropic elastic medium. As a result, strains (defined as the fractional change in length) are homogeneous and can be represented using symmetric second rank tensors.

where ∂ is the Lagrangian strain and εij are homogeneous strain parameters and i and j indicate Cartesian components44.

When we cast a crystal structure's Bravais lattice vectors, under an isotropic pressure, in a matrix form (R), the small homogeneous deformation (strain) distorts the Bravais lattice vectors of this crystal and hence cause a distortion of the lattice (R’) expressed by multiplying with a symmetric deformation matrix i.e. (R’ = R*D), where D is:

where I is a unique matrix that represents the symmetric strain tensor.

Now, we express the total energy of a crystal (R’), under strain in terms of a power series of the Lagrangian strain (∂):

where E (E0) is the energy and V (V0) is the volume of strained system.

The current \({\text{RFe}}_{3}\) compounds (R = rare earth) crystallize in the \({\text{PuNi}}_{3}\) -type structure, with space group 194-P63/mmc (Hexagonal).

C11, C12, C13, C33, and C55 are the five independent elastic constants for a hexagonal symmetry. We need five different strains to determine these elastic constants because we have five independent elastic constants. The five distortions are discussed further below. The first distortion46 is as follows:

and it modifies the basal plane while keeping the z-axis constant. As a result, the symmetry of the strained lattice remains hexagonal, and the energy for this distortion can be calculated as follows:

The second type of distortion is a volume-conserved distortion and lead to orthorhombic symmetry and written as:

and the energy for this distortion can be obtained as:

The third strain we have used is given by:

This strain changes C lattice parameter and keeps the strained lattice's symmetry hexagonal, and the energy for this distortion can be obtained by

The fourth elastic constant, C55, is determined by a lattice deformation that results in a low-symmetry object. The deformation is written as follows:

and it leads to triclinic symmetry and the energy for this deformation can be written as:

Finally, the last strain we have used is volume-conserved, keeps the symmetry of the strained lattice hexagonal, and can be written as:

and the energy for this strain is given by

Isotropic elastic constants such as the bulk, shear, and Young moduli can be determined using appropriate averaging procedures. All the previous parameters including the Young modulus, E, and the Poisson ratio, ν, the bulk modulus, B, and the shear modulus, G are calculated by IRelast Package47 as available in the Wien2k Package48. As for the crystal structure of the two compounds: YFe3 crystallizes in the hexagonal structure with space group P63/mmc and HoFe3 crystallizes in the trigonal structure with \({\overline{\text{R}}\text{3m}}\) space group. All the information about the lattice parameters and the locations of the Y or Ho and Fe atoms in the unit cells is available32.

The total effective fields of the R and Fe sublattices are expressed as follows by molecular field theory (MFT):

In these equations \({\overline{\text{H}}}\) represents the applied field, \(\overline{\mu }_{{\text{F}}}\) the magnetic moment per Fe ion at temperature T in units of the Bohr magneton (\(\mu_{{\text{B}}}\) ), and \(\overline{\mu }_{{\text{R}}}\) the magnetic moment per rare-earth ion. The factor d converts the moment per \({\text{RFe}}_{3}\) unit from \(\mu_{{\text{B}}}\) to gauss:

where NA is Avogadro's number, ρ the density in \({\text{g}}/{\text{cm}}^{3}\) , and A is the \({\text{RFe}}_{3}\) formula weight in g/mol. With these definitions the fields: \({\overline{\text{H}}}\) ,\({\overline{\text{H}}}_{{\text{R}}}\) and \({\overline{\text{H}}}_{{\text{F}}}\) are specified in gauss, and the molecular field coefficients \({\text{n}}_{{{\text{RR}}}}\) ,\({\text{n}}_{{{\text{RF}}}}\) and \({\text{n}}_{{{\text{FF}}}}\) respectively, describing the R-R, R-Fe, and Fe–Fe magnetic interactions, are dimensionless.

We begin with the magnetic energy of a binary magnetic compound to calculate the magnetic specific heat.

The magnetic specific heat is calculated as follows:

The magnetic entropy is calculated from the numerical integration of the magnetic heat capacity as follows:

The total heat capacity \({\text{C}}_{{{\text{total}}}}\) is made up of three contributions: the lattice heat capacity \({\text{ C}}_{{\text{l}}}\) , the electronic heat capacity \({\text{C}}_{{\text{e}}}\) and the magnetic heat capacity \({\text{ C}}_{{\text{m}}}\) :

The lattice heat capacity is expressed as:

where \({\text{x}} = \theta_{{\text{D}}} /{\text{T}}\) and \(\theta_{{\text{D}}}\) is Debye temperature, which can be calculated using Eq. (23)49

where h is the Planck's constant, \({\text{k}}_{{\text{B}}}\) is the Boltzmann constant,\({\text{v}}_{{\text{m}}}\) is the average sound velocity, and n is the number of atoms per formula unit50. Equation (24) gives the average sound velocity in a polycrystalline material:

where \({\text{v}}_{{\text{l}}}\) and \({\text{v}}_{{\text{s}}}\) are, respectively, the longitudinal and transverse sound velocities, which can be calculated using the shear and the bulk moduli G and B, from Navier’s equation51:

The electronic heat capacity is given by:

where \(\gamma_{{\text{e}}} = \frac{{\pi^{2} {\text{k}}_{{\text{B}}}^{2} {\text{D}}\left( {{\text{E}}_{{\text{F}}} } \right)}}{3}\) represents the electronic heat-capacity coefficient, and \({\text{D}}\left( {{\text{E}}_{{\text{F}}} } \right)\) represents the electron density of states at the Fermi energy \({\text{E}}_{{\text{F}}}\) .

A magnetic conducting material's total entropy includes three contributions: lattice (\({\text{S}}_{{\text{l}}}\) ), electronic (\({\text{S}}_{{{\text{el}}}}\) ) and magnetic entropy (\({\text{S}}_{{\text{m}}}\) ). The total entropy can be calculated by numerically integrating the total heat capacity:

The lattice entropy is given by:

The electronic entropy is given, like the electronic heat capacity, by Eq. (26).

The MCE is inherent in all magnetic materials and is induced by the magnetic lattice coupling with the applied magnetic field. The MCE is distinguished by two parameters: the adiabatic temperature change ΔTad (T, ΔH) and the isothermal entropy change \(\Delta {\text{S}}_{{\text{m}}}\) (T, ΔH). In an isothermal process, the entropy change caused by a magnetic field variation from H1 to H2, in compounds undergoing a second order phase transition, is calculated using the following Maxwell relation52:

The adiabatic temperature change \(\Delta {\text{T}}_{{{\text{ad}}}}\) is calculated from the following equation:

It is clear from Eq. (31) that the total heat capacity is applied field-dependent. However if \({\text{C}}\left( {{\text{T}},{\text{H}}_{{\text{F}}} } \right)\) has a weak field dependence, we can rewrite Eq. (31) as follows53:

The Wien2K electronic structure code calculated the total DOS (density of states) for YFe3 and HoFe3 as shown in Figs. 1 and 2 respectively.

The spin-up and spin-down electronic density of states (DOS) for YFe3.

The spin-up and spin-down electronic density of states (DOS) for HoFe3.

The total DOS, at \({\text{E}}_{{\text{f}}}\) for these two compounds are 10.915 and 24.461 states/eV, respectively. Therefore the electronic heat capacity coefficients, as calculated from Eq. 26, are 0.0256 and 0.0575 J/mol. K2 respectively. The DOS of both compounds show that these two systems are metallic with fairly large density of states, at Fermi energy, in both the spin up and spin down configuration. This, of course, demands taking the electronic contribution to the heat capacity into consideration as we pointed out (Eq. 26).

Using the two-sublattice molecular field theory, the temperature dependence of magnetization of the rare earth and Fe sublattices, as well as total magnetization for \({\mathrm{YFe}}_{3}\) and \({\mathrm{HoFe}}_{3}\) , are calculated. Total magnetic moments for \({\mathrm{YFe}}_{3}\) and \({\mathrm{HoFe}}_{3}\) ,calculated using the ab initio method are 4.8 µB/f.u and 4.4 µB /f.u, respectively, which are in good agreement with available experimental values, i.e. 4.88 and 4.59 µB /f.u as reported by J. F. Herbst et al.22. Table 1 displays the magnetic moments calculated in the present work and those by Ref.22. Figures 3 and 4a display mean-field- calculated temperature-dependence of magnetization, in zero field for \({\mathrm{YFe}}_{3}\) and \({\mathrm{HoFe}}_{3}\mathrm{respectively}\) . The total magnetic moments of these two compounds, at very low temperatures, are in excellent agreement with the values referred to in Table 1. In addition we displayed, in Figs. 3 and 4b, both our mean field calculated magnetization and the experimental data extracted from Ref.22.

Total magnetization of YFe3 in zero field vs. temperature.

(a) Mean-field calculated total and sublattice magnetization of HoFe3, in zero field, vs. temperature. (b) Total magnetization of HoFe3 vs. temperature.

Ferromagnetic order for YFe3 is evident from Fig. 3 and the magnetic moment per Fe atom is about 1.6 µB in agreement with the value in Fig. 1 in Herbst work22. On the other hand, a ferrimagnetic order with a compensation point is found for HoFe3. As the magnetic moments of the Fe and Ho sublattices become equal and antiparallel, a cancellation of the total moment takes place close to 400 K. The continuous decrease in magnetization at Tc is a well-known feature of SOPT materials54.

The temperature dependence of the magnetization, in 1, 3, 5 and 7 T fields, for \({\mathrm{HoFe}}_{3}\) , is shown in Fig. 5. an increase in the Curie temperature and a decrease in the compensation temperature is evident, with increasing the field. This behavior is consistent with that found in ferrimagnetic compounds, as reported for example, by P. von Ranke55. We may notice that even as the field increases the magnetization still has a continuous drop at Tc. Again, this confirms the SOPT nature of the studied compounds.

Total magnetization of HoFe3 in different applied fields 1, 3, 5 and 7 T vs. temperature.

As previously stated, three contributions must be calculated to determine the total heat capacity of these materials. Figures 6 and 7 show, respectively, the magnetic heat capacity for \({\mathrm{YFe}}_{3}\) and \({\mathrm{HoFe}}_{3}\) in different applied fields. The features of these two figures, at the Curie temperature, show that the transition is a second order phase transition. In contrast, in the FOPT case, the peak in the heat capacity shifts to higher temperatures as the field increases53.

Magnetic heat capacity of \({\mathrm{YFe}}_{3}\) , in different fields, vs. temperature.

Magnetic heat capacity of HoFe3 in different fields, vs. temperature.

The calculated magnetic entropy, in fields in the range 0–6 T, is shown in Figs. 8 and 9 for YFe3 and HoFe3 respectively. The entropy saturates to its maximum at Tc and above, and is reduced by increasing the field. The maximum values of the calculated magnetic entropy are about 27.3 and 50.9 J/mol. K for YFe3 and HoFe3, respectively. These values are calculated using the equation:

The calculated magnetic entropy of \({\mathrm{YFe}}_{3}\) in two different fields vs. temperature.

The calculated magnetic entropy of \({\mathrm{HoFe}}_{3}\) in different fields vs. temperature.

where Jr and JFe are the total angular momenta for the rare earth and Fe atoms, respectively. The molecular-field results of the maximum magnetic entropy is in very good agreement with the results of the equation shown above. The gradual increase of the magnetic entropy with temperature until it saturates at T ≥ Tc instead of the sudden increase found in FOPT materials show that the transition in these two compounds is a SOPT one.

The elastic properties of YFe3 are calculated by using IRelast Package as available in the Wien2k Package47, the output of the program is shown in Table 2.

The Debye temperature is calculated from the mean sound velocity. The bulk and shear moduli of \({\mathrm{YFe}}_{3}\) , using the Hill model, are 99.337 and 101.211 GPa, respectively. The calculated \({\theta }_{D}\) for \({\mathrm{YFe}}_{3}\) is 500.5 K. For the HoFe3 compound, the Debye temperature is 357.85 as reported by the Materials Project site (Kristin Persson)32.

Figures 10 and 11 show the isothermal entropy change, for different magnetic fields, for \({\mathrm{YFe}}_{3}\) and \({\mathrm{HoFe}}_{3}\) respectively. The \(\Delta \mathrm{Sm}\) curve for \({\mathrm{HoFe}}_{3}\) exhibits two peaks: the first is a broad peak below the compensation temperature, and the second, smaller peak, has its maximum at the ferrimagnetic-paramagnetic phase transition at \({T}_{C}=590 K\) . These two features correspond to the inverse and direct MCE effects, respectively. For ferromagnetic \({\mathrm{YFe}}_{3}\) there is only one peak at a temperature around its Tc (545 K). The temperature and field dependences of ∆Sm are those of SOPT materials. In particular, the curves at different fields have their maxima at Tc. In FOPT materials, the peak shifts to higher temperatures as the field increases53.

Temperature dependence of isothermal change in entropy for YFe3 in fields: 2, 3, 4, 5 and 6 T.

Temperature dependence of isothermal change in entropy for HoFe3 in fields of 1, 2 and 3 T.

We have calculated the adiabatic change in temperature for different field changes, as shown in Figs. 12 and 13, for YFe3 and HoFe3 systems, respectively. It is clear that the former compound has a higher cooling rate than the latter in agreement with reported literature56. For example, cooling rates about 1.3 and 0.4 K/T are achieved for these two compounds respectively for a field change of 3 T. The curves in each of Figs. 12 and 13 have their maxima centred at the Curie temperature. This is a known feature of compounds exhibiting second order phase transition and treated via the mean-field theory53.

The field-dependence of the adiabatic change in temperature vs. temperature for \({\mathrm{YFe}}_{3}\) .

The field-dependence of the adiabatic change in temperature vs. temperature for \({\mathrm{HoFe}}_{3}.\)

The RCP is a figure-of-merit for MCE materials (e.g.56). It is defined as RCP (T) = ∆Tad (max) × δTFWHM. Table 3 displays the RCP (T) for benchmark materials e.g. Gd, Gd-based compounds and FeRh system. From this data, we may conclude that the RCP of YFe3 is comparable with those of well-known materials.

The Arrott plots for YFe3 are shown, in a temperature range around Tc, in Fig. 14. The positive slopes at those temperatures, below and above Tc, are indicative of second order phase transition. First order transitions exhibit negative or s-shaped slopes57,58. The straight line starting near the origin is calculated at a temperature close to Tc (Fig. 3). The features of the universal curve, described below, supports the presence of SOPT as well.

The Arrott plots, at different temperatures, for YFe3.

The universal curves59, for YFe3, are shown in Fig. 15 for field changes of 3, 4, 5 and 6 T. The curves are collapsed on each other especially at high temperatures i.e. θ > 0. The parameter θ is given by: θ = (T − Tc) / (Tr − Tc), where Tr is a reference temperature defined as the temperature at which the following condition holds: ∆SM (Tr) = 0.7 (∆SM)peak. The collapse of the ΔSM curves is also indicative of second order phase transition, which is different from the features encountered in first order phase transitions60.

The universal curve for YFe3.

It would be of interest to support our finding that the mean field theory is suitable for explaining the physical properties. Therefore, calculating some of the critical exponents61,62,63,64,65 would be beneficial. We have calculated the exponent n in the relation:

from the field dependence of ΔSm close to Tc. The relation between n, β and γ is:

n = 1 + (β−1)/(β + γ) [e.g. 65, which can be used to evaluate γ if β is calculated from the relation:

or calculate β, if γ is calculated from:\(\chi_{T} \sim t^{ - \gamma } ,H = 0, T > T_{c}\) , where t = (T−Tc)/T. We have calculated β for YFe3 from its M-T data, in zero field, by taking the logarithm of M ~ (−t)β and performing the calculation for temperatures lower than Tc. We have obtained β around 0.52 for this compound. Therefore, the percentage error between our result and the value 0.5 of the mean-field is around 4%. The factor γ turned out to be around 0.92 for YFe3. In addition, we have calculated the critical exponent δ from the relation: δ = 1 + γ/β66 and found that δ = 2.8, which is close to the MFT value of 3. From the above analysis, we conclude that mean-field theory has fairly produced the critical exponents for YFe3. We would also emphasize the fact that the mean-field theory has its own limitations near Tc.

We should mention that other models e.g. 3D-Heisnberg, 3D- Ising and Tri-critical mean field model [e.g.67] may be used to study the critical exponents, however a future experimental and/or theoretical work, dedicated to this task, may shed more light on the most proper model.

The two-sublattice molecular field model was used for calculating the thermomagnetic properties for \({\mathrm{YFe}}_{3}\) and \({\mathrm{HoFe}}_{3}\) . The temperature dependence of magnetization shows that \({\mathrm{YFe}}_{3}\) is ferromagnetic with a Curie temperature close to 537 K, while \({\mathrm{HoFe}}_{3}\) is a ferrimagnetic compound with a compensation point around 389 K and a Curie temperature close to 565 K. The bulk and shear moduli are calculated using the WIEN2K ab-initio electronic code for YFe3 system. Those moduli were used to calculate its Debye temperature (500.53 K). The Debye temperature of HoFe3 is obtained from the Materials Project site. Using the calculated DOS, at Fermi energy, the electronic heat capacity coefficient \({\gamma }_{e}\) was found to be 0.0256 and 0.0575 J/\({K}^{2}\) mole, for \({\mathrm{YFe}}_{3}\) and \({\mathrm{HoFe}}_{3}\) respectively. The isothermal change in entropy ∆Sm, for a field change of 3 T, is about 0.8 and 0.12 J/mole. K for these systems, respectively. The Y-system exhibits a larger adiabatic drop in temperature (1.3 K/T), for a 3 T field change, than the Ho-system (0.4 K/T). The relative cooling power RCP(T) of YFe3 is comparable to well-known MCE materials The temperature and field dependences of the magnetization, magnetic heat capacity, entropy, the MCE quantities together with the Arrott plots and universal curve, are all indicative of SOPT.

The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

- Tishin, A. M. Handbook of Magnetic Materials, edited by KHJ Buschow North Holland. In Chapter 4: Magnetocaloric Effect in the Vicinity of Phase Transitions, 395–524. https://doi.org/10.1016/S1567-2719(99)12008-0 (1999).

Pecharsky , VK & Gschneidner Jr , KA Magnetocaloric effect and magnetic refrigeration .J. Magn.Magn.Mater.200(1–3), 44–56.https://doi.org/10.1016/S0304-8853(99)00397-2 (1999).

Article  CAS  ADS  Google Scholar 

O’Keefe, T. J., Roe, G. J. & James, W. J. X-ray investigation of the iron-holmium system. J. Less Common Metals 15(3), 357–360. https://doi.org/10.1016/0022-5088(68)90196-3 (1968).

Article  CAS  Google Scholar 

Tharp, D. E., Long, G. J., Pringle, O. A., James, W. J. & Grandjean, F. A Mössbauer effect study of ErFe3. J. Magn. Magn. Mater. 104–107, 1477–1478. https://doi.org/10.1016/0304-8853(92)90670-J (1992).

Article  ADS  Google Scholar 

Ray, A. & Strnat, K. Easy directions of magnetization in ternary R2 (Co, Fe)17 phases. IEEE Trans. Magn. 8(3), 516–518. https://doi.org/10.1109/TMAG.1972.1067471 (1972).

Article  CAS  ADS  Google Scholar 

Ross, J. W. & Crangle, J. Magnetization of cubic laves phase compounds of rare earths with cobalt. Phys. Rev. 133(2A), 509–510. https://doi.org/10.1103/PhysRev.133.A509 (1964).

Article  ADS  Google Scholar 

Wertheim, G. K. & Wernick, J. H. Mössbauer effect in some iron-rare Earth intermetallic compounds. Phys. Rev. 125(6), 1937. https://doi.org/10.1103/PhysRev.125.1937 (1962).

Article  CAS  ADS  Google Scholar 

Marei, S. A., Craig, R. S., Wallace, W. E. & Tsuchida, T. Magnetic characteristics of some Laves phase systems containing Fe and Mn. J. Less Common Metals 13(4), 391–398. https://doi.org/10.1016/0022-5088(67)90033-1 (1967).

Article  CAS  Google Scholar 

- Wallace. W. E., Skrabek, E. A. In Rare Earth Research Volume 2, Proceedings of the Third Rare Earth Conference, edited by K. S. Vorres (Gordon and Breech, New York, 1964), p. 417. 431.

- Roe, G. J. The iron-holmium binary system and related magnetic properties. "Doctoral Dissertations". 1876. https://scholarsmine.mst.edu/doctoral_dissertations/1876 (1969).

- Hoffer. G. I., Salmans, L. R. 7th Rare Earth Research Conf. Cornado, California 371 (1968).

Földeàki, M., Chahine, R. & Bose, T. K. Magnetic measurements: a powerful tool in magnetic refrigerator design. J. Appl. Phys. 77(7), 3528–3537. https://doi.org/10.1063/1.358648 (1995).

Article  ADS  Google Scholar 

Li, L. & Yan, M. Recent progresses in exploring the rare earth based intermetallic compounds for cryogenic magnetic refrigeration. J. Alloys Compd. 823, 153810. https://doi.org/10.1016/j.jallcom.2020.15381 (2020).

Article  CAS  Google Scholar 

Zhang, Y. Review of the structural, magnetic and magnetocaloric properties in ternary rare earth RE2T2X type intermetallic compounds. J. Alloy. Compd. 787, 1173–1186. https://doi.org/10.1016/j.jallcom.2019.02.175 (2019).

Article  CAS  Google Scholar 

De Oliveira, N. A., von Ranke, P. J. Theoretical aspects of the magnetocaloric effect. Phys. Rep. 489(4–5), 89–159. doi:https://doi.org/10.1016/j.physrep.2009.12.006 (2010).

Nikitin, S. A., Talalaeva, E. V., Chernikova, L. A. & Andreenko, A. S. Magnetocaloric effect in compounds of rare-earth metals with iron. Soviet J. Exp. Theor. Phys. 38, 1028 (1974).

ADS  Google Scholar 

Nagy, A., Hammad, T., Yehia, S. & Aly, S. H. Thermomagnetic properties and magnetocaloric effect of TmFe2 compound. J. Magn. Magn. Mater. 473, 324–330. https://doi.org/10.1016/j.jmmm.2018.10.050 (2019).

Article  CAS  ADS  Google Scholar 

Tian, G. et al. Large reversible magnetocaloric effect of light rare-earth intermetallic compound Pr5Si3. J. Alloy. Compd. 496, 517–520. https://doi.org/10.1016/j.jallcom.2010.02.093 (2010).

Article  CAS  Google Scholar 

Clark, A. E. Magneto-strictive rare earth-Fe2 compounds. Handb. Ferromagn. Mater. 1, 531–589. https://doi.org/10.1016/S1574-9304(05)80122-1 (1980).

Article  ADS  Google Scholar 

Campbell, I. A. Indirect exchange for rare earths in metals. J. Phys. F: Met. Phys. 2(3), L47. https://doi.org/10.1088/0305-4608/2/3/004 (1972).

Article  CAS  ADS  Google Scholar 

Craig, R. S., Wallace, W. E. & Kevin-Smith, H. Hydrogenation of intermetallic compounds: thermodynamic aspects. Sci. Technol. Rare Earth Mater. 55, 353. https://doi.org/10.1016/B978-0-12-675640-1.50021-0 (1980).

Article  Google Scholar 

Herbst, J. F. & Croat, J. J. Magnetization of RFe3 intermetallic compounds: molecular field theory analysis. J. Appl. Phys. 53(6), 4304–4308. https://doi.org/10.1063/1.331207 (1982).

Article  CAS  ADS  Google Scholar 

Herbst, J. F. & Croat, J. J. Magnetization of R6Fe23 intermetallic compounds: molecular field theory analysis. J. Appl. Phys. 55(8), 3023–3027. https://doi.org/10.1063/1.333293 (1984).

Article  CAS  ADS  Google Scholar 

Aly, S. H. & Yehia, S. A computer simulation study on a two-sublattice magnetic system. J. Magn. Magn. Mater. 213(3), 383–388. https://doi.org/10.1016/s0304-8853(99)00805-7 (2000).

Article  CAS  ADS  Google Scholar 

Shafiq, M., Ahmad, I. & Jalali-Asadabadi, S. Theoretical studies of strongly correlated rare-earth intermetallics RIn3 and RSn3 (R= Sm, Eu, and Gd). J. Appl. Phys. 116(10), 103905. https://doi.org/10.1063/1.4894833 (2014).

Article  CAS  ADS  Google Scholar 

Blaha, P. et al. WIEN2k: An APW+lo program for calculating the properties of solids. J. Chem. Phys. 152, 074101. https://doi.org/10.1063/1.5143061 (2020).

Article  CAS  PubMed  ADS  Google Scholar 

Debye, P. On the theory of specific heats.Ann.Phys.344(14), 789-839 (1912).

Article  MATH  Google Scholar 

- Kittel, C. Elementary Statistical Physics. Courier Corporation (2004).

Kittel, C., McEuen, P. (1996) Introduction to Solid State Physics, vol. 8. New York: Wiley.

Andersen, O. K. Linear methods in band theory. Phys. Rev. B 12(8), 3060. https://doi.org/10.1103/PhysRevB.12.3060 (1975).

Article  CAS  ADS  Google Scholar 

Stadler, R. et al. Ab initio calculations of the cohesive, elastic, and dynamical properties of CoSi2 by pseudopotential and all-electron techniques. Phys. Rev. B 54, 1729. https://doi.org/10.1103/PhysRevB.54.1729 (1996).

Article  CAS  ADS  Google Scholar 

- Persson, K. A. https://materialsproject.org/

Kohn, W. & Sham, L. J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 140(4A), A1133. https://doi.org/10.1103/PhysRev.140.A1133 (1965).

Article  MathSciNet  ADS  Google Scholar 

- Singh, D. J., Nordstrom, L. Planewaves, Pseudopotentials, and the LAPW method, 2nd Edition. Springer (2006).

Perdew, J. P. & Wang, Y. Accurate and simple analytic representation of the electron-gas correlation energy. Phys. Rev. B 45(23), 13244. https://doi.org/10.1103/PhysRevB.45.13244 (1992).

Article  CAS  ADS  Google Scholar 

Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77(18), 3865. https://doi.org/10.1103/PhysRevLett.77.3865 (1996).

Article  CAS  PubMed  ADS  Google Scholar 

Singh, D. J. Ground-state properties of lanthanum: Treatment of extended-core states. Phys. Rev. B 34(8), 6388. https://doi.org/10.1103/PhysRevB.43.6388 (1991).

Article  ADS  Google Scholar 

Mohammad, F. A., Yehia, S. & Aly, S. H. A first-principle study of the magnetic, electronic and elastic properties of the hypothetical YFe5 compound. Physica B: Condens. Matter 407(13), 2486–2489. https://doi.org/10.1016/j.physb.2012.03.050 (2012).

Article  CAS  ADS  Google Scholar 

Shabara, R. M. & Aly, S. H. A first-principles study of elastic, magnetic, and structural properties of PrX2 (X=Fe, Mn, Co) compounds. J. Magn. Magn. Mater. 423, 447–452. https://doi.org/10.1016/j.jmmm.2016.09.075 (2017).

Article  CAS  ADS  Google Scholar 

Aly, S. H. & Shabara, R. M. First principles calculation of elastic and magnetic properties of Cr-based full-Heusler alloys. J. Magn. Magn. Mater. 360, 143–147. https://doi.org/10.1016/j.jmmm.2014.02.030 (2014).

Article  CAS  ADS  Google Scholar 

Elalfy, G. M., Shabara, R. M., Aly, S. H. & Yehia, S. First-principles study of magnetic, electronic, elastic and thermal properties of GdFe2.". Comput. Condens. Matter 5, 24–29. https://doi.org/10.1016/j.cocom.2015.10.001 (2015).

Article  Google Scholar 

Habbak, E. L., Shabara, R. M., Aly, S. H. & Yehia, S. Investigating half-metallicity in PtXSb alloys (X=V, Mn, Cr, Co) at ambient and high pressure. Physica B: Condens. Matter 494, 63–70. https://doi.org/10.1016/j.physb.2016.04.009 (2016).

Article  CAS  ADS  Google Scholar 

Shabara, R. M., Yehia, S. & Aly, S. H. Pressure-induced phase transitions, electronic and magnetic properties of GdN. Results Phys. 1, 30–35. https://doi.org/10.1016/j.rinp.2011.08.001 (2011).

Article  CAS  ADS  Google Scholar 

Abu-Elmagd, M. S. M., Aly, S. H. & Yehia, S. Magnetic properties and electronic structure of Tℎ Co4B. Model. Numer. Simul. Mater. Sci. 2(3), 51–59. https://doi.org/10.4236/mnsms.2012.23006 (2012).

Article  CAS  Google Scholar 

Yehia, S., Aly, S. H. & Aly, A. E. Electronic band structure and spin-density maps of SmCo5. Comput. Mater. Sci. 41(4), 482–485. https://doi.org/10.1016/j.commatsci.2007.05.004 (2008).

Article  CAS  Google Scholar 

Jamal, M., Kamali Sarvestani, N., Yazdani, A. & Reshak, A. H. Mechanical and thermodynamical properties of hexagonal compounds at optimized lattice parameters from two-dimensional search of the equation of state. RSC Adv. 4(101), 57903–57915. https://doi.org/10.1039/C4RA09358E (2014).

Article  CAS  ADS  Google Scholar 

Jamal, M., Bilal, M., Ahmad, I. & Jalali-Asadabadi, S. IRelast package.J. Alloy.Compd.735, 569–579.https://doi.org/10.1016/j.jallcom.2017.10.139 (2018).

Article  CAS  Google Scholar 

- Jamal, M. IRelast. http://www.wien2k.at/ (2019).

Foster, A. S., Sulimov, V. B., Lopez Gejo, F., Shluger, A. L. & Nieminen, R. M. Structure and electrical levels of point defects in monoclinic zirconia. Phys. Rev. B 64(22), 224108. https://doi.org/10.1103/PhysRevB.64.224108 (2001).

Article  CAS  ADS  Google Scholar 

Anderson, O. L. A simplified method for calculating the Debye temperature from elastic constants. J. Phys. Chem. Solids 24(7), 909–917. https://doi.org/10.1016/0022-3697(63)90067-2 (1963).

Article  CAS  ADS  Google Scholar 

- Schreiber, E., Anderson, O. L., Soga, N., Bell, J. F. Elastic constants and their measurement, 747–748. https://doi.org/10.1115/1.3423687(1975).

Morrish, A. H. The Physical Principles of Magnetism. New York: Wiley. https://doi.org/10.1002/9780470546581 (2001).

Smith, A. et al. Materials challenges for high performance magnetocaloric refrigeration devices. Adv. Energy Mater. 2(11), 1288–1318. https://doi.org/10.1002/aenm.201200167 (2012).

Article  CAS  Google Scholar 

Lyubina, J. Magnetocaloric materials for energy efficient cooling. J. Phys. D: Appl. Phys. 50, 053002. https://doi.org/10.1088/1361-6463/50/5/053002 (2017).

Article  CAS  ADS  Google Scholar 

Von Ranke, P. J., Alho, B. P., Nobrega, E. P. & de Oliveira, N. A. Understanding the inverse magnetocaloric effect through a simple theoretical model. Physica B 404(19), 3045–3047. https://doi.org/10.1016/j.physb.2009.07.009 (2009).

Article  CAS  ADS  Google Scholar 

- Gschneidner Jr., K. A., Pecharsky, V. K. Intermetallic Compounds: Principles and Practice Vol. 3, Ch.25, Wiley. https://doi.org/10.1002/0470845856.ch25 (2002).

Law, J. Y. et al. A quantitative criterion for determining the order of magnetic phase transitions using the magnetocaloric effect. Nat. Commun. 9(1), 1–9. https://doi.org/10.1038/s41467-018-05111-w (2018).

Article  CAS  Google Scholar 

Foldeaki, M. et al. Comparison of magnetocaloric properties from magnetic and thermal measurements. J. Appl. Phys. 82(1), 309–316. https://doi.org/10.1063/1.365813 (1997).

Article  CAS  ADS  Google Scholar 

Franco, V. & Conde, A. Scaling laws for the magnetocaloric effect in second order phase transitions: from physics to applications for the characterization of materials. Int. J. Refrig. 33(3), 465–473. https://doi.org/10.1016/j.ijrefrig.2009.12.019 (2010).

Article  CAS  Google Scholar 

Bonilla, C. M. et al. Universal behavior for magnetic entropy change in magnetocaloric materials: an analysis on the nature of phase transitions. Phys. Rev. B 81(22), 224424. https://doi.org/10.1103/PhysRevB.81.224424 (2010).

Article  CAS  ADS  Google Scholar 

Magnetic critical scattering, Malcolm F. Collins, Theory of critical phenomena: critical exponents, table 3.1 on p. 12 (Oxford University Press, 1989).

- Long Range Order in Solids, Robert M. White, and Theodore H. Geballe, ch.1: Phase transitions and order parameters (Academic Press, 1979).

- Magnetism in the Solid State: an introduction, Peter Mohn, ch.5, Springer 2006.

- Quantitative Theory of Critical Phenomena, George A. Baker, Jr. ch.1 (Academic press 1990).

El Ouahbi, A., Yamkane, Z., Derkaoui, S. & Lassri, H. Magnetic properties and the critical exponents in terms of the magnetocaloric effect of amorphous Fe40Ni38Mo4B18 alloy. J. Superconductivity Novel Magn 34, 1253. https://doi.org/10.1007/s10948-021-05832-y (2021).

Article  CAS  Google Scholar 

Widom, B. "Equation of state in the neighborhood of the critical point. J. Chem. Phys. 43, 3898. https://doi.org/10.1063/1.1696618 (1965).

Article  ADS  Google Scholar 

Bally, M. A. A. et al. Magnetocaloric properties and analysis of the critical point exponents at PM-FM phase transition. Results Phys. 28, 104546. https://doi.org/10.1016/j.rinp.2021.104546 (2021).

Article  Google Scholar 

Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB). The research leading to these results received funding from [ASRT, Egypt] under Grant Agreement No [6701]. ASRT is the second affiliation of this research.

Department of Engineering Mathematics and Physics, Higher Institute of Engineering, El Shourok Academy, Cairo, Egypt

Department of Physics, Faculty of Science, Helwan University, Cairo, Egypt

Tarek Hammad, Ahmed Abdel-Kader & Sherif Yehia

Physics Department, Faculty of Science, Taibah University, Madinah, Saudi Arabia

Department of Physics, Faculty of Women, Ain Shams University, Cairo, Egypt

Department of Physics, Faculty of Science, Damietta University, New Damietta, Egypt

Samy H. Aly & Fatema Z. Mohammad

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

M.S.M.A.-E.: Performing the ab initio calculation, writing, reviewing. T.H.: Reviewing calculation and manuscript. A.A.- K.: Assisted in carrying out the MCE calculation using Mathematica, preparing figures. N.E.-S.: Carrying out the trapezoidal calculation, reviewing the manuscript. S.Y.: reviewing all the calculation, writing the original manuscript. SH.A.: Suggesting the systems to work on, reviewing the calculation and the original manuscript. F.Z.M.: performing the magnetocaloric calculation using Mathematica, reviewing the ab initio calculation, reading the manuscript.

Correspondence to Mohammed Said Mohammed Abu-Elmagd.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Abu-Elmagd, M.S.M., Hammad, T., Abdel-Kader, A. et al. First principles and mean field study on the magnetocaloric effect of YFe3 and HoFe3 compounds. Sci Rep 13, 2876 (2023). https://doi.org/10.1038/s41598-023-29676-9

DOI: https://doi.org/10.1038/s41598-023-29676-9

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

First principles and mean field study on the magnetocaloric effect of YFe3 and HoFe3 compounds | Scientific Reports

Cerium Powder Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.