Next Article in Journal
The Influence of Diet on Oxidative Stress and Inflammation Induced by Bacterial Biofilms in the Human Oral Cavity
Next Article in Special Issue
Study of the Structure, Magnetic, Thermal and Electrical Characterisation of ZnCr2Se4: Ta Single Crystals Obtained by Chemical Vapour Transport
Previous Article in Journal
Micro-/Nano-Structured Ceramic Scaffolds That Mimic Natural Cancellous Bone
Previous Article in Special Issue
Independently Tunable Multipurpose Absorber with Single Layer of Metal-Graphene Metamaterials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Enhancement of Hard Magnetic Properties in Fraktal-Like Nano and Mesoscopic Grains

by
Grzegorz Ziółkowski
1,*,
Dariusz Chrobak
2,
Grażyna Chełkowska
1,
Ondrej Zivotsky
3 and
Artur Chrobak
1
1
Institute of Physics, University of Silesia in Katowice, 75 Pułku Piechoty 1A, 41-500 Chorzów, Poland
2
Institute of Materials Engineering, University of Silesia in Katowice, 75 Pułku Piechoty 1A, 41-500 Chorzów, Poland
3
Department of Physics, VŠB-Technical University of Ostrava, 17 listopadu 15/2172, 708 33 Ostrava-Poruba, Czech Republic
*
Author to whom correspondence should be addressed.
Submission received: 17 February 2021 / Revised: 6 March 2021 / Accepted: 12 March 2021 / Published: 16 March 2021

Abstract

:
The paper refers to Monte Carlo magnetic simulations for fractal-like nano and mesoscopic grains. The analyzed objects differed in the size, surface development, magnetic anisotropy and the spin values attributed to the system nodes inside the fractal. Such an approach allowed us to determine their magnetization processes as well as optimization characteristics in the direction to enhancement of hard magnetic properties. As it was shown, the size effects depend on the chosen value of magnetic anisotropy. In the case of fractals with ultra-high coercivity, the decreasing of their size leads to deterioration of coercivity, especially for the high surface to volume ratio. Opposite effects were observed for soft magnetic fractals when the nanostructure caused an appearance of the coercive field, and the maximum of energy product was predictably significantly higher than for conventional rare earths’ free permanent magnets.

1. Introduction

Magnetic materials are very important in modern technologies such as electromobility, magnetic cooling, or energy harvesting. Especially interesting are magnets with extreme soft or hard magnetic properties. It is well known that iron-based nanocrystalline alloys exhibit excellent soft magnetic properties (i.e., close to zero coercivity and high magnetic saturation) [1,2] that are widely applicable in electric motors and transformers. On the other side of the magnetic materials application spectra, there are hard magnets characterized by high value of the so-called energy product [3,4,5,6,7,8]. In this area, the Fe-RE-B (RE (rare earths), typically Nd, Nd/Dy, Nd/Tb) type of alloys are considered the best, showing the energy product in order of 400 kJ/m3. It is interesting that the classical RE-free permanent magnets reveal values of this parameter no higher than 30 kJ/m3 (for example—ALNICO alloys [9,10]). For economic reasons, many scientific teams have worked on searching for new RE-free/reduced hard magnetic materials that can fill the gap between these two kinds of magnets [11]. One promising research direction is nanotechnology, allowing us to obtain nanosized magnetic systems and enabling the introduction of an additional magnetic anisotropy by the interaction between nano-objects [12,13,14,15]. Surface and size effects are also important. In this context, objects with fractal-like geometry seem to be very interesting [16,17]. In some cases, procedure of fractal reproduction can be considered as a model for early-stage crystal growth, allowing generating objects with a specific surface shape. Recently, we undertook a preliminary study into the magnetization processes of such objects, i.e., 3D nanosized fractals [18]. It turned out that it is possible to control magnetic characteristics by a specific combination of magnetic properties in volume, on the surface, and by the shape of the fractal objects. The results presented in this work motivated us to significantly broaden the range of analysis, including the mesoscopic size of the objects and optimization of the energy product, that is desired for permanent magnets applications.
In this paper we have focused on simulations of magnetization processes of nano and mososcopic objects characterized by different shape and size, and different magnetic properties of the objects in volume and surface. In particular, we proposed the analysis of fractal-like objects with different ratios between magnetic moments placed in the volume and on the surface. Such objects can be considered as an equivalent of crystal grains formed during the early stages of crystallization, for example, from amorphous or liquid state. We have analyzed the magnetic characteristics of the fractals that are described by magnetically soft, hard (coercivity about 1T), and ultrahigh-coercive volume (coercivity about 8T, experimentally observed and reported in [19]) as well as different values of surface magnetic anisotropy. Such an approach allowed us to determine magnetic hysteresis loops for the studied systems, and additionally, the approach made it possible to suggest a way for designing relatively good hard magnetic materials towards reduction or elimination of RE elements in their content.

2. Simulation Procedure

Generally, we used the Monte Carlo Metropolis method as the simulation procedure [20]. Recently, we developed a computational approach which allows the analysis of magnetization processes of hard magnetic systems containing nano and mesoscopic 3D objects [21,22]. We introduced the following two elements: (i) modification of the cluster Wolff algorithm [23] by an additional factor taking into account information (or configuration) entropy originating from a difference in magnetic properties of system nodes [21] and (ii) formulation of scaling rules for micromagnetic simulation of mesoscopic objects [22].
All simulated systems in this work were prepared based on the diffusion limited aggregation algorithm (DLA), described in detail in [24,25,26,27,28]. Basically, in this method, a single particle (seed) has been placed in the center of the system, and then the new particles are introduced into the system starting from the edges. The new particle moves randomly and, with some probability, it may stick to another stationary particle encountered in the system. Moreover, the algorithm was extended by a growing procedure (existing particles may extend to adjacent fields) similar to the grain growing procedure observed in the real systems. For example, the DLA can simulate chemical deposition, MBE deposition or rapid solidification from liquid state [29,30,31,32]. By controlling the attachment probability and the growth speed, it is possible to obtain the systems with a constant number of particles (7000), but with a different degree of the surface development, i.e., the ratio of the number of particles on the surface (NS) to the volume (NV), as was shown in Figure 1.
The generated fractals were embedded in 3D space divided by 40 × 40 × 40 nodes, and digitization of the system lies in assigning magnetic properties to each node. In particular, the nodes are described by a “spin vector” S, the definition of which is represented by the node magnetic moment, and, therefore, it can have any value and any direction. The size of the fractals were defined by a spacial distance between the nodes for the initial systems and the scaling factors. Based on the three presented systems, i.e., NS/(NS + NV) equal to 0.3, 0.5 and 0.7 (further called surface ratio), a series of simulation systems with different anisotropy values of the nodes on the surface, KS = 0 eV, 5 × 10−4 eV, 5 × 10−3 eV (with easy magnetization axes perpendicular to the surface), as well as in the volume KV = 0 eV, 5 × 10−5 eV and 5 × 10−4 eV (with easy magnetization axes parallel to the z-axis) were prepared. The initial distance between the nodes were 0.28 nm in order to simulate material with similar to α-Fe density. Moreover, the value of the spin varied from 0.1 to 1, and the scaling factor n (see the next section) varied from 1 to 100, which allows us to analyze nano and meso-sclaled objects.
The main simulation procedure consists of a series of the following repeated steps:
  • Select a random i-th node for further analysis and choose a random number r1 ∊ [0,1).
  • If r1 < Pcl (i.e., cluster analysis probability), find a cluster around selected node and change cluster directions by θ angle, otherwise only the direction of the selected node is changed. In order to simulate multiphase ferromagnetic systems, it is necessary to use the disorder-based cluster MC method (which we described in [21]) with the adding probability equal to:
    P i j a d d = 1 e x p E i j c o u p l i n g k B T e x p α S i l o c
    where E i j c o u p l i n g is the direct exchange coupling energy between the spins attributed to nodes i and j, S i l o c is the local configuration entropy of anisotropy (calculated in the defined sphere around the i-th node) and α is the factor responsible for strengthening and weakening of the entropy impact on the adding probability (can be randomly selected from a selected range).
  • Calculate the difference of the system energy ΔE before and after the change in step 2. The energy is calculated based on the 3D Heisenberg model:
    E = i , j J i j S i S j i K i s ^ i n ^ i 2 g μ B μ 0 i H i S i + D i , j S i S j 3 S i e i j S j e i j r i j 3
    where Jij is the exchange parameter, Si is the spin vector on site i, Ki is the anisotropy constant (per site), n ^ i is the versor of the easy magnetization axis, s ^ i is the versor of the spin, g is the Lande factor, µB is the Bohr magneton, µ0 is the vacuum permeability, Hi is the magnetic field on site i, D is the dipolar constant, eij is the directional versor between the i-th and j-th nodes and rij is the distance between them.
  • Choose a random number r2 ∊ [0,1) and if r2 > e x p E k B T restore all changes made in step 2, otherwise accept the new configuration.
It is worth noting that the basic energy equation presented in Equation (2) is extended by the scaling rules (described in details in [22]), in order to model a large-scale mesoscopic systems, with an acceptable level of computing resource consumption. The main idea lies in the assumption that one node in the rescaled system can represent a finite volume consisting of n × n × n original (unscaled) nodes. In this situation, the parameters appearing in Equation (2) should be recalculated using the formulation of scaling rules in [22]. Table 1 summarizes the main simulation parameters and the scaling rules depending on the scaling factor n.

3. Results and Discussion

Magnetization processes of the objects in question were studied by carrying out a set of simulations resulting in the so-called reverse magnetization curves, i.e., gradual remagnetization of the object in opposite direction of external magnetic field after full magnetic saturation. The reverse magnetization curves were simulated for all the examined objects. Simultaneously, spin configuration was recorded in order to show coercivity mechanisms responsible for the remagnetization process. Figure 2 depicts such curves for the surface ratio equal to 0.7, S = 1 and for full ranges of the other simulation parameters.
One can see that in some cases the reverse magnetization process occurs in a rapid jump, which is obviously caused by a coherent rotation of the whole spins in the fractal. In contrast to this, when the surface anisotropy is high, a graduate change of magnetization was observed. The influence of the anisotropies and scaling factor on the coercive field HC is interesting. For the fractals with ultra-high volume anisotropy, the size enlargement leads to increase of HC while for the other cases, this effect is not very strong and even opposite. For deeper analysis it is worth presenting dependencies of the scaling factor on HC for all studied cases, as shown in Figure 3. In the case of the soft magnetic fractals (i.e., KV = 0), the appearance of relatively high coercivity requires a simultaneous occurrence of their low nanometric size and high value of surface magnetic anisotropy.
In a contrast to this, for the fractals with ultra-high coercivity (i.e., KV = 5 × 10−4 eV and KS = 5 × 10−4 eV) the increasing scaling factor results in the significant increase of coercive field, especially for the surface ratio equal to 0.7. An explanation of this behavior can be found in the recorded spin configurations. Figure 4 and Figure 5 depict magnetic moments of the nodes placed on a central plane (KV = 5 × 10−4 eV, KS = 5 × 10−4 eV and KV = 0, KS = 5 × 10−3 eV, respectively) of the analyzed 3D fractals for the surface ratio 0.7 and with the applied external magnetic fields around the coercive field. It should be emphasized that the surface anisotropy can play different roles dependently on the KV values. When the volume anisotropy is high, the magnetic moments of the developed surface can change their directions in lower external magnetic field due to different angles between the KS vector and the field H; e.g., some of the surface magnetic moments are nearly parallel to the field. If magnetic coupling between volume and surface is strong enough, the anisotropy field of the whole fractal decreases. Therefore, the reverse magnetization process is initialized on surface, as shown in Figure 4. Furthermore, enlargements of the fractals to mesoscopic size is attributed to a relative decrease of the NS/NV ratio, which obviously leads to reducing the meaning of the surface magnetic moments and the increase in HC is observed (see Figure 3, middle-top graph). Relatively high surface anisotropy can cause the completely opposite effect when the volume of the fractal is magnetically soft (see Figure 3, bottom-right graph). In this situation, only magnetic moments on the surface are a source of anisotropy and via exchange interactions they can influence magnetic moments in the volume. Indeed, Figure 5 depicts such a case in which the reverse magnetization process is initiated by the soft magnetic nodes in the volume of the fractal. The effect of interaction between the volume and the surface is clearly visible in Figure 5d when the spins are aligned not perfectly co-linear, as was observed for ultra-high coercive objects (see Figure 4d).
The current part of our study can be summarized as follows. Unexpectedly, enhancement of hard magnetic properties for ultra-high coercive fractal objects is impossible in the application of higher development of the surface. It is clear that real nanosized magnets cannot occur without significant contribution of magnetic surface anisotropy. Therefore, such materials should be designed as mesoscopic grains, accounting for the fact that high contribution of the surface caused deterioration of their hard magnetic properties. For the remaining kind of fractals, i.e., with magnetic hard and soft properties, a decrease in their size together with high surface anisotropy can cause a significant increase in the coercive field.
The indicated possibility to create hard magnetic properties in soft magnetic fractal is extremely intriguing and worthy of further investigations. It is not a new idea; e.g., the well-known ALNICO alloys contain nanosized iron grains in the form of elongated cigars. Reviewing the literature, one can find a value of iron surface anisotropy in the order of 10−4 eV [33,34], so the use in our simulation of KS = 5 × 10−4 eV is a realistic value. Unfortunately, for S = 1 we did not observe any occurrence of coercivity in analyzed fractals. The reason for this could be to high magneto-static energy of the volume in a comparison with the surface anisotropy energy. The only way to obtain the desired enhancement of hard magnetic properties effect is to decrease the S values. For this purpose, we chose the fractal with the surface ratio equals 0.7, KS = 5 × 10−4 eV, KV = 0 eV and S ranging from 0.1 to 1. For example (i.e., for S = 0.2 and 0.5), Figure 6 and Figure 7 show the simulated reverse magnetization curves and the resulting HC as a function of the scaling factor, respectively. As it was expected, the decrease of the fractal magnestostatic energy (the third term in Equation (2)) caused the complexity of the reverse magnetization process that, in the case of n = 1, leads to an appearance of coercivity.
Analyzing the full variability of S, one can notice an increase of coercive field with decreasing S value. Possible applications as permanent magnets also require high values of the magnetic remanence and, consequently, the energy product |BH|max. Let us estimate the |BH|max values, taking density and magnetic induction as for α-Fe (i.e., induction in saturation equal to 2 T for S = 1). With this assumption, it is possible to plot dependencies of the S value on the coercive field and the estimated energy product, as shown in Figure 8 for n = 1 and n = 10.
For both quantities, the presented curves reveal maxima at different spin values. In particular, for S = 0.7 and n = 1 the energy product |BH|max equals 85 kJ/m3, which is very promising considering that, for the best rare earths free permanent magnets (the ALNICO alloys), this parameter is about 30 kJ/m3. This means that a specific nanostructure with highly developed surface may fill the gap between conventional and rare earths-based permanent magnets. The open question is how to obtain such a structure and how to control magnetic moment of atoms in iron-based magnetic compounds. One can imagine applying rapid crystallization (flash annealing, mold suction casting, melt-spinning techniques, etc.) to control the microstructure as well as alloying addition (for example, some metaloids) to control magnetic moments attributed to Fe atoms.
In order to complete this study, the simulated hysteresis loop as well as the spin configurations for the fractal exhibiting maximum of |BH|max are presented in Figure 9 and Figure 10, respectively.
When the volume of the fractal is magnetically soft, the reverse magnetization process starts from the volume inside the system, similarly to the other soft magnetic fractal systems (see Figure 5). In these cases, the resulting non-zero anisotropy field appears as a competition between the volume magnetostatic energy and the magnetic anisotropy energy on the surface.

4. Concluding Remarks

Summarizing the presented simulations, one can state that the fractal like magnetic objects are very interesting regarding the possibility to optimize their hard magnetic properties. As it was shown, a combination of the fractal characteristics (i.e., size, surface development, magnetic parameters in volume and on surface) can lead to a change of coercivity mechanism. This effect might cause either improvement or deterioration of the fractal hard magnetic properties. In the case of ultra-high coercive objects, nanosized structure and fractal-like shape is not favorable because the surface magnetic moments make the reverse magnetization process more “easy”. When the anisotropy field of the volume is not extremely high, the change of coercivity mechanism may result in the observed increase of coercivity. Such a desired effect requires low dimensional fractals with a highly developed surface as well as relatively high surface magnetic anisotropy.
Quite interesting results have been obtained in the case of fractals with magnetically soft volume. Despite the fact that the presented analysis is only a rough estimation, one can define the system parameters that resulted in significant enhancement of hard magnetic properties, filling the gap between the conventional and rare-earth-based permanent magnets.

Author Contributions

Conceptualization, G.Z., A.C.; methodology, G.Z., A.C., D.C.; software G.Z.; validation, O.Z., G.C.; formal analysis, O.Z., G.C.; writing—original draft preparation, G.Z., A.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Science Centre in Poland with the grant 2015/19/B/ST8/02636. D. Chrobak acknowledges support from the National Science Centre of Poland, grant no. 2016/21/B/ST8/02737. O. Zivotsky acknowledges support from the project SP2021/64 under Ministry of Education, Youth and Sports of the Czech Republic.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McHenry, M.E.; Willard, M.A.; Laughlin, D.E. Amorphous and nanocrystalline materials for applications as soft magnets. Prog. Mater. Sci. 1999, 44, 291. [Google Scholar] [CrossRef]
  2. Herzer, G.; Varga, L.L. Structure and coactivity of nanocrystalline Fe-Si-B-Nb-Cu alloys. J. Magn. Magn. Mater. 2000, 215/216, 506–512. [Google Scholar]
  3. Luborsky, F.E. Permanent Magnets in Use Today. J. Appl. Phys. 1966, 37, 1091–1094. [Google Scholar] [CrossRef] [Green Version]
  4. Sagawa, M.; Fujimura, S.; Yamamoto, H.; Matsuura, Y.; Hiraga, K. Permanent magnet materials based on the rare earth-iron-boron tetragonal compounds. IEEE Trans. Magn. 1984, 20, 1584–1589. [Google Scholar] [CrossRef] [Green Version]
  5. Croat, J.J.; Herbst, J.F.; Lee, R.W.; Pinkerton, F.E. Pr-Fe and Nd-Fe based materials: A new class of high performance permanent magnets. J. Appl. Phys. 1984, 55, 2078–2081. [Google Scholar] [CrossRef]
  6. Gutfleisch, O. Controlling the properties of high energy density permanent magnetic materials by different processing routes. J. Phys. D Appl. Phys. 2000, 33, 157–172. [Google Scholar] [CrossRef]
  7. Coey, J.M.D. Perspective and Prospects for Rare Earth Permanent Magnets. Engineering 2020, 6, 119–131. [Google Scholar] [CrossRef]
  8. Ziółkowski, G.; Randrianantoandro, N.; Chrobak, A.; Klimontko, J.; Kędziołka-Gaweł, M.; Haneczok, G. Influence of Transition and Rare Earth Elements on Magnetic Properties of Fe-Nb-B-M (M = Ni, Ag, Gd, Tb) Bulk Nanocrystalline Alloys. Acta Phys. Pol. A 2012, 121, 1266–1269. [Google Scholar] [CrossRef]
  9. Rao, A.S. Alnico permanent magnets an overview. In Proceedings of the Electrical/Electronics Insulation Conference, Chicago, IL, USA, 4–7 October 1993. [Google Scholar] [CrossRef]
  10. Zhou, L.; Tang, W.T.; Ke, L.; Guo, W.; Poplawsky, J.D.; Anderson, I.E.; Kramer, M.J. Microstructural and magnetic property evolution with different heat-treatment conditions in an alnico alloy. Acta Mater. 2017, 133, 73–80. [Google Scholar] [CrossRef] [Green Version]
  11. Coey, J.M.D. Permanent magnets: Plugging the gap. Scr. Mater. 2012, 67, 524–529. [Google Scholar] [CrossRef]
  12. Yuan, X.H.; Zhao, G.P.; Yue, M.; Ye, L.N.; Xia, J.; Zhang, X.C.; Chang, J. 3D and 1D calculation of hysteresis loops and energy products for anisotropic nanocomposite films with perpendicular anisotropy. J. Magn. Magn. Mater. 2013, 343, 245–250. [Google Scholar] [CrossRef]
  13. Pellicelli, R.; Solzi, M.; Pernechele, C.; Ghidini, M. Continuum micromagnetic modeling of antiferromagnetically exchange-coupled multilayers. Phys. Rev. B 2011, 83, 054434. [Google Scholar] [CrossRef]
  14. Pellicelli, R.; Solzi, M.; Neu, V.; Häfner, K.; Pernechele, C.; Ghidini, M. Characterization and modeling of the demagnetization processes in exchange-coupled SmCo5/Fe/SmCo5 trilayers. Phys. Rev. B 2010, 81, 184430. [Google Scholar] [CrossRef]
  15. Sadovnikov, A.V.; Grachev, A.A.; Gubanov, V.A.; Odintsov, S.A.; Martyshkin, A.A.; Sheshukova, S.E.; Sharaevskii, Y.P.; Nikitov, S.A. Spin-wave intermodal coupling in the interconnection of magnonic units. Appl. Phys. Lett. 2018, 112, 142402. [Google Scholar] [CrossRef]
  16. Komogortseva, S.V.; Iskhakova, R.S.; Felk, V.A. Fractal Dimension Effect on the Magnetization Curves of Exchange-Coupled Clusters of Magnetic Nanoparticles. J. Exp. Theor. Phys. 2019, 128, 754–760. [Google Scholar] [CrossRef]
  17. Zou, G.; Xiong, K.; Jiang, C.; Li, H.; Li, T.; Du, J.; Qian, Y. Fe3O4 Nanocrystals with Novel Fractal. J. Phys. Chem. B 2005, 109, 18356–18360. [Google Scholar] [CrossRef]
  18. Ziółkowski, G.; Chrobak, A. Magnetization processes of irregular dendrite structures-A Monte Carlo study. Phys. B 2020, 577, 411745. [Google Scholar] [CrossRef]
  19. Chrobak, A.; Zółkowski, G.; Randrianantoandro, N.; Klimontko, J.; Chrobak, D.; Prusik, K.; Rak, J. Ultra-high coercivity of (Fe86−xNbxB14)0.88Tb0.12 bulk nanocrystalline magnets. Acta Mater. 2015, 98, 318–326. [Google Scholar] [CrossRef]
  20. Metropolis, N.; Rosenbluth, A.W.; Teller, A.H.; Teller, E. Equation of State Calculations by Fast Computing Machines. J. Chem. Phys. 1953, 21, 1087. [Google Scholar] [CrossRef] [Green Version]
  21. Chrobak, A.; Ziółkowski, G.; Granek, K.; Chrobak, D. Disorder-based cluster Monte Carlo algorithm and its application in simulations of magnetization processes. Comput. Phys. Commun. 2019, 238, 157–164. [Google Scholar] [CrossRef]
  22. Chrobak, A.; Ziółkowski, G.; Chrobak, D.; Chełkowska, G. From atomic level to large-scale Monte Carlo magnetic simulations. Materials 2020, 13, 3696. [Google Scholar] [CrossRef]
  23. Wolff, U. Collective Monte Carlo Updating for Spin Systems. Phys. Rev. Lett. 1989, 62, 361. [Google Scholar] [CrossRef]
  24. Witten, T.A.; Sander, L.M. Diffusion-Limited Aggregation: A Kinetic Critical Phenomenon. Phys. Rev. Lett. 1981, 47, 1400–1403. [Google Scholar] [CrossRef]
  25. Sander, L.M. Fractal growth processes. Nature 1986, 322, 789–793. [Google Scholar] [CrossRef] [Green Version]
  26. Witten, T.A.; Sander, L.M. Diffusion-limited aggregation. Phys. Rev. B 1983, 27, 5686. [Google Scholar] [CrossRef] [Green Version]
  27. Sander, L.M.; Cheng, Z.M.; Richter, R. Diffusion-limited aggregation in three dimensions. Phys. Rev. B 1983, 28, 6394–6396. [Google Scholar] [CrossRef]
  28. Bourke, P.D. Constrained Diffusion Limited Aggregation in 3 Dimensions. Comput. Graph. 2006, 30, 646–649. [Google Scholar] [CrossRef]
  29. Levi, A.; Kotrla, M. Theory and Simulation of Crystal Growth. J. Phys. Condens. Matter 1996, 9, 299–344. [Google Scholar] [CrossRef]
  30. Hadj-Achour, M.; Brutin, D. Fractal pattern formation in nanosuspension sessile droplets via evaporation-spreading on a glass substrate. J. Colloid Interface Sci. 2014, 1, 43–46. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, M.; Braun, H.G.; Meyer, E. Crystalline structures in ultrathin poly(ethylene oxide)/poly(methyl methacrylate) blend films. Polymer 2003, 44, 5015–5021. [Google Scholar] [CrossRef]
  32. Lad, K.; Maaroof, M.; Raval, K.G.; Pratap, A. Fractal Growth Kinetics during crystallization of amorphous Cu50Zr50. Prog. Cryst. Growth Charact. Mater. 2002, 45, 15–19. [Google Scholar] [CrossRef]
  33. Takayama, H.; Bohnen, K.P.; Fulde, P. Magnetic surface anisotropy of transition metals. Phys. Rev. B 1976, 14, 2287. [Google Scholar] [CrossRef]
  34. Gradmann, U. Magnetic surface anisotropies. J. Magn. Magn. Mater. 1986, 54–57, 733–736. [Google Scholar] [CrossRef]
Figure 1. Geometry of systems obtained based on the diffusion limited aggregation (DLA) algorithm with different ratio of the number of particles on the surface (NS) to the volume (NV).
Figure 1. Geometry of systems obtained based on the diffusion limited aggregation (DLA) algorithm with different ratio of the number of particles on the surface (NS) to the volume (NV).
Materials 14 01443 g001
Figure 2. Reverse magnetization curves for all examined objects with NS/(NS + NV) = 0.7 and S = 1.
Figure 2. Reverse magnetization curves for all examined objects with NS/(NS + NV) = 0.7 and S = 1.
Materials 14 01443 g002
Figure 3. Coercive field HC as a function of the scaling factor for all studied cases with S = 1.
Figure 3. Coercive field HC as a function of the scaling factor for all studied cases with S = 1.
Materials 14 01443 g003
Figure 4. Spin configurations (cross-section) of the system with KS = 5 × 10−4 eV, KV = 5 × 10−4 eV, n = 1, S = 1 and NS/(NS + NV) = 0.7 in the external magnetic field equal to (a) 0 T, (b) −1.6 T, (c) −2 T and (d) −2.2 T.
Figure 4. Spin configurations (cross-section) of the system with KS = 5 × 10−4 eV, KV = 5 × 10−4 eV, n = 1, S = 1 and NS/(NS + NV) = 0.7 in the external magnetic field equal to (a) 0 T, (b) −1.6 T, (c) −2 T and (d) −2.2 T.
Materials 14 01443 g004
Figure 5. Spin configurations (cross-section) of the system with KS = 5 × 10−3 eV, KV = 0 eV, n = 1, S = 1 and NS/(NS + NV) = 0.7 in the external magnetic field equal to (a) 6.2 T, (b) 0 T, (c) −1.6 T and (d) −3.6 T.
Figure 5. Spin configurations (cross-section) of the system with KS = 5 × 10−3 eV, KV = 0 eV, n = 1, S = 1 and NS/(NS + NV) = 0.7 in the external magnetic field equal to (a) 6.2 T, (b) 0 T, (c) −1.6 T and (d) −3.6 T.
Materials 14 01443 g005
Figure 6. Simulated reverse magnetization curves for the systems with KS = 5 × 10−4 eV, KV = 0 eV and S = 0.2 as well as S = 0.5.
Figure 6. Simulated reverse magnetization curves for the systems with KS = 5 × 10−4 eV, KV = 0 eV and S = 0.2 as well as S = 0.5.
Materials 14 01443 g006
Figure 7. Coercive field HC as a function of the scaling factor for the systems presented in Figure 6.
Figure 7. Coercive field HC as a function of the scaling factor for the systems presented in Figure 6.
Materials 14 01443 g007
Figure 8. Dependencies of the S value on the coercive field and the estimated energy product for the systems with KS = 5 × 10−4 eV, KV = 0 eV, NS/(NS + NV) = 0.7 and n = 1 as well as n = 10.
Figure 8. Dependencies of the S value on the coercive field and the estimated energy product for the systems with KS = 5 × 10−4 eV, KV = 0 eV, NS/(NS + NV) = 0.7 and n = 1 as well as n = 10.
Materials 14 01443 g008
Figure 9. Simulated hysteresis loop for the fractal with the highest |BH|max i.e., S = 0.7, n = 1, KS = 5 × 10−4 eV, KV = 0 eV and NS/(NS + NV) = 0.7. The a–d points correspond to the spin configurations in Figure 10.
Figure 9. Simulated hysteresis loop for the fractal with the highest |BH|max i.e., S = 0.7, n = 1, KS = 5 × 10−4 eV, KV = 0 eV and NS/(NS + NV) = 0.7. The a–d points correspond to the spin configurations in Figure 10.
Materials 14 01443 g009
Figure 10. Spin configurations (cross-section) of the system presented in Figure 9 in the external magnetic field equal to (a) 1 T, (b) 0 T, (c) −0.4 T and (d) −0.6 T.
Figure 10. Spin configurations (cross-section) of the system presented in Figure 9 in the external magnetic field equal to (a) 1 T, (b) 0 T, (c) −0.4 T and (d) −0.6 T.
Materials 14 01443 g010
Table 1. Values of the simulation parameters and scaling roles used in this work.
Table 1. Values of the simulation parameters and scaling roles used in this work.
ParameterBase ValuesScaling Rules
Scaling factorfrom 1 to 100-
Distance between nodes0.28 nmr’ = r × n
Dipolar constant2.15 × 10−7 eVnm3D’ = D
Anisotropy constant (surface)from 0 eV to 5 × 10−3 eVK’S = KS × n2
Anisotropy constant (volume)from 0 eV to 5 × 10−4 eVK’V = KV × n3
Spinfrom 0.1 to 1S’ = S × n3
Exchange integral parameter1.5 × 10−3 eVJ’ = J × n−4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ziółkowski, G.; Chrobak, D.; Chełkowska, G.; Zivotsky, O.; Chrobak, A. Enhancement of Hard Magnetic Properties in Fraktal-Like Nano and Mesoscopic Grains. Materials 2021, 14, 1443. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14061443

AMA Style

Ziółkowski G, Chrobak D, Chełkowska G, Zivotsky O, Chrobak A. Enhancement of Hard Magnetic Properties in Fraktal-Like Nano and Mesoscopic Grains. Materials. 2021; 14(6):1443. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14061443

Chicago/Turabian Style

Ziółkowski, Grzegorz, Dariusz Chrobak, Grażyna Chełkowska, Ondrej Zivotsky, and Artur Chrobak. 2021. "Enhancement of Hard Magnetic Properties in Fraktal-Like Nano and Mesoscopic Grains" Materials 14, no. 6: 1443. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14061443

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