The Anisotropic Magnetocaloric Effect and Size-Dependent Magnetic Properties of Iron Particles

We present a theoretical study on the anisotropic magnetocaloric effect and the size-dependent magnetic properties of Fe particles of radii in the range 25–150 Å. An observable increase has been found in the magnetization, of the low radii (25–75 Å) particles, by reducing the temperature to 4 K. The anisotropic isothermal change in entropy ΔSm has been calculated by taking the difference between maximum ΔSm along the easy [100] and hard [111] directions. The maximum anisotropic ΔSm is 0.015 J/kg K for a field change of 500 Oe along the [100] direction. The ΔSm temperature dependence exhibits a table-like plateau for small radii (25–75 Å) and in low fields below 300Oe. This enhances the relative cooling power (RCP) of the Fe element to be 8.11 J/kg for particles of 25 Å radius. Also, the calculation of anisotropic ΔTad was performed along the easy axis and showed an increase in the maximum value around 37% relative to the experimental conventional value.

The fast development of nanometer-scale electronics demands the development of future cooling systems that can operate in complex geometries in the nanoscale range. An increase in the saturation magnetization in transition elements has been reported, when particle size is reduced to the nano-scale, at very low temperatures [15][16][17].
Research efforts, both experimental and theoretical, have been reported on the conventional MCE in the 3d transition elements at low and high fields, e.g., [14,16,18]. On the other hand, the anisotropic MCE (AMCE) is of interest because an increase in the entropy change could be achieved by rotating the crystal, along a specific crystalline direction, without changing the magnetic field intensity. This property was studied first by Kuzmin and Tishin in DyAlO 3 [17].
The application of magnetocaloric materials, either ordinary or anisotropic, in magnetic refrigeration has been the subject of many papers [14,[18][19][20][21][22][23]; in particular, search for materials with magneto-volume coupling and first-order phase transition is currently an active research area. Other applications of MCE materials are in the domain of medicine, e.g., in treating hyperthermia and in drug delivery [15,[24][25][26][27][28], where nanoparticles could be used.
In the present work, we report on the anisotropic sizedependent magnetic properties and the anisotropic magnetocaloric effect of crystalline cubic Fe, for spherical particles of radii in the 25-150-Å range, and at temperatures up to T c . A statistical mechanics-based model is used, in the current work, for calculating the size-dependent magnetic and magnetocaloric properties of Fe. This model has been successfully used in studying a host of magnetic properties of magnetic systems of different crystalline symmetries. Examples of those systems are hexagonal, tetragonal, and mixed-anisotropy systems [29][30][31][32][33]. The aforementioned magnetic properties of those systems include magnetization curves, magnetic susceptibility, energy, torque, and probability landscapes.

Model and Analysis
The interplay between different magnetic energies dictates the magnetic properties of a given magnetic system and the total energy is subsequently used in evaluating magnetothermal and magnetocaloric properties using the well-known relations of statistical mechanics.
Many publications provide experimental values of the temperature dependence of the anisotropy constants: K 1 , K 2 , and K 3 of Fe [34][35][36][37][38][39][40]. There is almost an agreement between these references on the sign and order of magnitude of K 1 ; however, this trend is not found regarding the value and the sign of K 2 and K 3 . Therefore, we preferred to use the Callen and Callen power law [41,42] to calculate the temperature-dependent anisotropy constants. The Callen and Callen law has been checked and verified for elements, e.g., Fe and Co [41][42][43][44]. The values of K 1 and K 2 , at very low temperatures, used in our model were obtained from the reference by Getzlaff [40]. Those  where, for cubic systems, E a can be expressed as [45]: where: We assume that the saturation magnetization vector is oriented at an angle θ with the z-axis and its projection on the x-y plane makes an angle φ with the x-axis which is chosen to be along the [100] direction of the cubic crystal. Hence, the M s direction cosines are: α x = sinθ cos φ, α y = sinθ sinφ, and α z = cosθ. Using the dimensionless quantities: 3 , and ε z = VβH M S , we rewrite Eq. (1) as: where V, β, and ψ are the volume, 1/k B T, and the angle between H and M s respectively. The quantity ε z is the ratio between the Zeeman and the thermal energies. As a result, the value of ε z represents the relative strength of these two quantities [45]. Iron particles are assumed to be spherical in shape with radius r. The classical partition function is given by: and is subsequently used to evaluate magneto-thermal and magnetocaloric quantities.
The isothermal change in the magnetic entropy is given by: where M(H,T) is the isothermal magnetization along a specific crystalline direction.
Equation (6) could be casted to take the following form [8,[44][45][46][47][48][49]: In the reference by Pecharsky et al. [49], the authors discuss the limitations that should be imposed on Tand H in order for Eq. (7) to be valid. The same reference also discusses and comments on the analysis given in another reference [50] on the same subject.
The isothermal magnetization vs. magnetic field for Fe particles (r = 150 Å) at H = 800 Oe at room temperature. Our model (solid curves) experimental data (★ easy, · mid-hard, ○ hard directions) The adiabatic change in temperature for a specific change in the external magnetic field is given by many references, e.g., [4,7,8,12]: However, one may write Eq. (8) in the following approximate form if the heat capacity is weakly fielddependent [8,12]

M(T, H) T ) H dH
A figure-of-merit called the relative cooling power (RCP) has been introduced to assess the performance of magnetocaloric materials. The RCP is determined by multiplying the magnitude of the maximum ∆S m , at a definite field, times the difference in temperature, at full width at half maximum, of the ∆S m vs. temperature curve as shown in the following equation [8] : The isothermal magnetization curves at 300 K for Fe particle of different sizes at A 125 Å, B 100 Å, C 75 Å, D 50 Å, and E 25 Å magnetization Ms is about 1713 emu/cm 3 ( 220.17 emu/g), i.e., very close to the experimental value of 217.2 emu/g as reported by Coey [52] and Kittel [53]. The sequence of the easy [100], mid-hard [110], and hard [111] directions is in agreement with the well-known data on cubic Fe crystals, e.g., [53,54]. We have studied the effect of particle size on the magnetization curves for particles in the size range of 25-150 Å at room temperature. Figure 2A-E show the dependence of the magnetization curves on particle size in the range 25-125 Å. We noticed a relative decrease in the anisotropy between the [100] and [111] directions as the particle size is reduced (e.g., Fig. 2C). For a particle size of 25 Å, as shown in Fig. 2E, the magnetization curves along the three directions become identical and saturation (~ 1700 emu/cm 3 ) is far from being achieved in a field up to 800 Oe. Figure 3 displays the dependence of magnetization on particle size (r = 25-150 Å), at 300 K and in 800 Oe field. The magnetization drops considerably as the particle size is reduced below about 50 Å. This behavior is consistent with the experimental results on Fe [55][56][57] and other particles like Gd [29]. Figure 4A-D shows the isothermal magnetization curves parallel to the easy [100] axis for particles with radii in the 25-to 150-Å range and at temperatures 600, 400, 150, and 4 K, respectively. The field required to achieve The isothermal magnetization curves along [100] axis at A 600 K, B 400 K, C 150 K, and D 4 K for Fe particles of different sizes in the 25-150-Å range saturation increases as the particle size decreases, for these three temperatures 600, 400, and 150 K. In addition, for particle radii in the range 25-75 Å and at 4 K, a significant increase in the saturation magnetization has been achieved, e.g., for a radius of 25 Å, the magnetization is raised to ~ 1700 emu/cm 3 as shown in Fig. 4D. This behavior is expected because thermal energy is insufficient, at low temperatures, to reduce the magnetization of extremely small particles, whereas this effect becomes more significant at high temperatures. A common feature of super-paramagnetic single-domain particles is that at high temperatures, the thermal energy overcomes the magnetic anisotropy energy barrier [29]. The most remarkable feature of superparamagnetism is shown in Fig. 5, namely the   [58], and Gd element [29]. In contrast, Fig. 6

Anisotropic Magnetocaloric Effect (AMCE)
The anisotropic MCE effect is manifested in the difference between the isothermal entropy changes calculated along different crystallographic directions in a field of constant magnitude [59,60]. Figure 7A-D display the AMCE, using the trapezoidal method, in the field range 50-500 Oe, and at temperatures between 900 and 1050 K, for Fe particles of size 150 Å. The anisotropic ΔS anis can be calculated by taking the difference between ΔS m along the easy and hard directions. A large ΔS anis is achieved in the low field region 50-500 Oe because the entropy of the spins is changed by application of a small field, in cluster and superparamagnetic materials more easily than in paramagnetic systems [61]. Table 1 shows the maximum value of ΔS m at different fields 50, 100, 300, and 500 Oe for a radius of 150 Å. Figure 8 displays the temperature dependence of ΔS m for r = 150 Å in a field of 30 kOe. The maximum value is 8.36 J/kg K, along the easy axis, in comparison with 7.95 J/kg K for the bulk Fe [16]. So, the anisotropic MCE is enhanced over the ordinary MCE. Figures 9 and 10 display the size-dependent ΔS, in a 1kOe field along the [100] direction for radii in the range 25-150 Å. A plateau feature [51] is found for 25-Å particles almost in the whole temperature range. Table 2 shows the RCP values for radii 50, 75, 100, and 150 Å in a field of 1 kOe. The highest value of RCP is 8.11 J/kg for a radius of 50 Å.
This behavior was reported previously in nano-sized Gd element by Miller [62] and Belo et al. [21]. However, the previous studies have not taken into account the calculation of anisotropic ΔT ad and only the conventional ΔT ad was presented experimentally as reported by [8,15]. Based on Eq. (9), the anisotropic ΔT ad is calculated along the easy direction [100] for applied field of 30 kOe. The heat capacity data was obtained from the experimental work presented in [8]. The anisotropic ΔT ad has been calculated, using our model, with maximum value of 8 K as compared to the experimental value of 5 K as shown in Fig. 11. The increase in the value of anisotropic ΔT ad relative to ordinary MCE confirms the importance of applying a field along a specific crystallographic direction for obtaining a high adiabatic temperature change at the same field intensity.

Conclusions
The The statistical mechanics model proved to be fairly suitable for calculating the magnetothermal, magnetic sizedependent, and anisotropic magnetocaloric properties of the iron element. Data Availability This manuscript has no associated data or the data will not be deposited. The reader may contact me through my email: (eldocahmed4@gmail.com).
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:// creat iveco mmons. org/ licen ses/ by/4. 0/.