Effects of Shape Anisotropy on Hard–Soft Exchange-Coupled Permanent Magnets

Exchange-coupled magnets are promising candidates for a new generation of permanent magnets. Here, we investigated the effect of soft magnetic shell thickness and the aspect ratio of the hard magnetic core on the magnetic properties for isolated core/shell cylinder exchange-coupled magnets, as well as the packing effect of the cylindrical array via a micromagnetic simulation method. It was found that the shape anisotropy contributions to the magnetic properties in the cylindrical core/shell exchange-coupled magnets are closely related to the thickness of the soft magnetic shell. When the soft magnetic shell is thin, the magnetic properties are dominated by the hard–soft exchange coupling effects, and the contributions of shape anisotropy are quite limited. When the soft magnetic shell is relatively thick, utilizing shape anisotropy would be an effective method to improve the magnetic performance of hard–soft exchange-coupled magnets. The present work provides an in-depth fundamental understanding of the underlying magnetization reversal mechanism. This work could be useful for designing high-performance permanent magnets and avoiding pitfalls.


Introduction
Rare-earth permanent magnets bearing strong energy density have become vital components of modern electromagnetic technology [1][2][3], which are increasingly demanded in an impressive range of applications, including energy generation and conversion, as well as information storage and processing [4]. Today's energy product holder is Nd 2 Fe 14 B with Dy addition, which exhibits a high coercivity even at 373 K [5]. However, due to the global trade situation, rare-earth is in short supply [6], in particular, Dy, which is expensive and has low reserves even by rare-earth standards. Over the past decades, great efforts have been devoted to developing permanent magnets with higher performance and lower costs to solve this dilemma [7,8].
Exchange-coupled magnets are promising candidates for a new generation of permanent magnets that have drawn increasing attention [9]. These magnets contain two main magnetic phases, i.e., the high-anisotropy hard magnetic phase, and the high-saturationmagnetization soft magnetic phase, which are coupled via exchange interaction on the nanometer scale. The exchange-coupled magnets offer the unique opportunity to simultaneously utilize the high coercivity of the hard magnetic phase and the high saturation magnetization of the soft magnetic phase, and then to achieve high maximum magnetic energy product (up to slightly higher than 1 MJ m −3 , or twice the value for Nd 2 Fe 14 B) [10,11]. Meanwhile, the exchange-couple principle could also provide flexibility in material selection for the constituent hard and soft phases, which represents a possible solution for the supply risk of the strategic rare-earth elements [12]. Many studies have been conducted to improve the coercivity and magnetic energy product of the exchange-coupled magnets by controlling morphological and geometric parameters [8,[13][14][15]. Nevertheless, with the increase in the soft magnetic phase, the remanence of the system increases at the expense of the coercivity [5], and there is a trade-off between the coercivity and remnant magnetization for an optimal magnetic energy product [4], which remains a major challenge for exploiting high-performance permanent magnets.
In a single domain state at the nanometer scale, coercivity substantially higher than the magneto-crystalline anisotropy field could be achieved if the shape anisotropy can be enhanced because the shape anisotropy can provide an additional resource of coercivity by forming an effective anisotropy field [4,[16][17][18]. For example, the Co nanowires with a high coercivity of~1 T can be synthesized via numerous methods, including template-assisted electrodeposition, electroless deposition under an external magnetic field, and chemical synthesis [18]. Thus, the magnetic performance of the traditional hard-soft exchangecoupled magnets could be further improved by utilizing shape anisotropy. Furthermore, the 1-D nanostructured materials (nanowire, or nanorod) with enhanced shape anisotropy could prevent aggregation due to their unique shape characteristics, which can also be beneficial for more sufficient exchange-coupling compared to the traditional exchangecoupled nanoparticles [13].
Currently, micromagnetic simulations are widely used to study the magnetization phenomena on the intermediate scale between the quantum mechanical scale of individual atoms and the macroscale, such as modeling domain structures [19]. Herein, micromagnetic simulation provides a valuable tool to design exchange-coupled magnets, which is very helpful in terms of reducing the research cost and accelerating the design process [20,21]. In this work, according to the typical fabrication process of cylindrical hard core/soft shell exchange-coupled permanent magnets, we performed micromagnetic simulations. The effect of soft magnetic shell thickness and the aspect ratio of the hard magnetic core on the magnetic properties for the isolated core/shell cylinder exchange-coupled magnets was systematically investigated, and the underlying magnetization reversal mechanisms were discussed. In addition, considering that the entire magnets are assembled using amounts of isolated nanostructure units, the packing effect of the cylindrical array was also investigated. Our work may serve as the foundation for exploiting high-performance permanent magnets in the future.

Materials and Methods
Micromagnetic simulations were carried out using the micromagnetic simulator OOMMF (Object Oriented MicroMagnetic Framework) [22] by dynamically solving the Landau-Lifshitz-Gillbert (LLG) equation [23,24] within the framework of the finite difference method [25]: where α is the damping constant and |γ| is the gyromagnetic ratio, satisfying |γ| = |γ| 1+α 2 . The damping constant of α = 0.5 has been used for the simulation of the quasi-static case [18,19,26]. Considering that the typical bottom-up assembly fabrication of cylindrical core/shell nanostructured hard-soft exchange-coupled magnets involves two stagespreparation of the precursor hard magnetic cylinder (nanowire or nanorod), and subsequent deposition of the soft magnetic shell [13]-an isolated cylindrical core/shell nanostructured micromagnetic model was used, as shown in Figure 1a. The core/shell nanostructure has a hard magnetic core surrounded by a soft magnetic shell. The diameter and length of the hard magnetic cylinder is denoted as D and L, respectively, while the thickness of the soft magnetic shell is denoted as t. To investigate the contributions of exchange-coupling and shape anisotropy to the magnetic properties, we calculated the demagnetization curves of the cylindrical core/shell nanostructure with various shell thickness (t = 1, 2, 3, 4, 5 nm) and aspect ratio (L/D = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, and D = 10 nm). Both the applied field H (from +14 T to −14 T) and the magnetocrystalline anisotropic easy axis of the hard magnetic phase are assumed to be paralleled to the longitudinal direction (z-axis direction) of the cylinder. The simulation mesh size of 1 nm × 1 nm × 1 nm for t = 3, 4, 5 nm, and 0.5 nm × 0.5 nm × 0.5 nm for t = 1, 2 nm were used. To examine the effect of packing [12], hexagonal close-packed cylindrical array structures were used, as illustrated in Figure 1b. For effectively exchange-coupling, it is well known that a soft magnetic layer should be thinner than twice the domain-wall width of the hard magnetic phase. The Sm 2 Co 17 exhibits a larger domain-wall width, as compared to SmCo 5 and Nd 2 Fe 14 B, which is more flexible for exchange coupling [13]. Here, we selected Sm 2 Co 17 and Fe 65 Co 35 as the hard magnetic phase and soft magnetic phase, respectively, and the magnetic parameters including anisotropy K 1 , saturation magnetization M s and exchange constant A of Sm 2 Co 17 [27] and Fe 65 Co 35 [28] are as follows: lated the demagnetization curves of the cylindrical core/shell nanostructure w shell thickness (t = 1, 2, 3, 4, 5 nm) and aspect ratio (L/D = 1, 2, 3, 4, 5, 6, 7, 8, 9, 1 10 nm). Both the applied field H (from +14 T to −14 T) and the magnetocrystal tropic easy axis of the hard magnetic phase are assumed to be paralleled to the nal direction (z-axis direction) of the cylinder. The simulation mesh size of 1 nm 1 nm for t = 3, 4, 5 nm, and 0.5 nm × 0.5 nm × 0.5 nm for t = 1, 2 nm were used. T the effect of packing [12], hexagonal close-packed cylindrical array structures as illustrated in Figure 1b. For effectively exchange-coupling, it is well known magnetic layer should be thinner than twice the domain-wall width of the hard phase. The Sm2Co17 exhibits a larger domain-wall width, as compared to S Nd2Fe14B, which is more flexible for exchange coupling [13]. Here, we selecte and Fe65Co35 as the hard magnetic phase and soft magnetic phase, respective magnetic parameters including anisotropy K1, saturation magnetization Ms and constant A of Sm2Co17 [27] and Fe65Co35 [28]

Effect of Soft Magnetic Shell on Magnetic Properties
To explore the effect of the soft magnetic shell thickness on magnetic prop demagnetization curves of the isolated core/shell cylinder with aspect ratio L/ calculated; the results are shown in Figure 2a. It should be noted that the asp defined as the ratio between the diameter and length of the hard magnetic co of the entire core/shell cylinder. Clearly, the shape of the demagnetization curv depends on the thickness of the soft magnetic shell. As the thickness of the sof shell increases, the magnetization reversal is clearly separated into two par slight decrease in magnetization at the lower applied field, and the jump at applied field), which deteriorates the squareness. This is because the mag

Effect of Soft Magnetic Shell on Magnetic Properties
To explore the effect of the soft magnetic shell thickness on magnetic properties, the demagnetization curves of the isolated core/shell cylinder with aspect ratio L/D = 1 were calculated; the results are shown in Figure 2a. It should be noted that the aspect ratio is defined as the ratio between the diameter and length of the hard magnetic core, instead of the entire core/shell cylinder. Clearly, the shape of the demagnetization curves strongly depends on the thickness of the soft magnetic shell. As the thickness of the soft magnetic shell increases, the magnetization reversal is clearly separated into two parts (i.e., the slight decrease in magnetization at the lower applied field, and the jump at the higher applied field), which deteriorates the squareness. This is because the magnetization reversal is expected to start from the soft magnetic phase prior to the magnetization of hard magnetic phase switching, as shown below. The saturation magnetization (M s ) in such soft-hard exchange-coupled magnets usually follows the general mixing rule, i.e., where M s soft and M s hard are the saturation magnetization of the soft magnetic phase and the hard magnetic phase, respectively; V soft and V hard are the volume fraction of the soft magnetic phase and the hard magnetic phase, respectively. As seen from Figure 2a, the M s increases with t increasing. It should be noted that saturation magnetization is an intrinsic material property, which is decisive for the remanence M r and the maximum magnetic energy product (BH) max [29]. The intrinsic coercivity H ci and the saturation magnetization determined from the demagnetization curves are shown as a function of soft magnetic shell thickness t in Figure 2b. As seen from Figure 2b, the intrinsic coercivity H ci decreases monotonously with t increasing. The exchange-coupled soft magnetic phase enhances the saturation magnetization but reduces the coercivity. This trade-off between saturation magnetization and coercivity is a fundamental challenge in the quest for high-performance exchange-coupled magnets.
hard magnetic phase switching, as shown below. The saturation magnetization (Ms) in such soft-hard exchange-coupled magnets usually follows the general mixing rule, i.e., where Ms soft and Ms hard are the saturation magnetization of the soft magnetic phase and the hard magnetic phase, respectively; Vsoft and Vhard are the volume fraction of the soft magnetic phase and the hard magnetic phase, respectively. As seen from Figure 2a, the Ms increases with t increasing. It should be noted that saturation magnetization is an intrinsic material property, which is decisive for the remanence Mr and the maximum magnetic energy product (BH)max [29]. The intrinsic coercivity Hci and the saturation magnetization determined from the demagnetization curves are shown as a function of soft magnetic shell thickness t in Figure 2b. As seen from Figure  2b, the intrinsic coercivity Hci decreases monotonously with t increasing. The exchangecoupled soft magnetic phase enhances the saturation magnetization but reduces the coercivity. This trade-off between saturation magnetization and coercivity is a fundamental challenge in the quest for high-performance exchange-coupled magnets. To gain an in-depth fundamental understanding of the magnetic reversal, we systematically investigated the effect of soft magnetic shell thickness on the magnetic structure evolution during demagnetization. The magnetic structures of different soft magnetic shell thickness prior to coercivity of the core/shell cylinder with L/D = 1 are shown in Figure 3. As can be seen, the magnetic reversal with different soft magnetic shell thickness exhibits three different characteristics. For the cases of t = 1 nm and t = 2 nm, the magnetic reversal begins at the soft magnetic shell, and the magnetization reversal mode is coherent rotation. When t = 3 nm, magnetization reversal occurs in a quasi-coherent rotation, where the partial magnetic moments of the fringe are tilted. However, when t reaches above 4 nm, the magnetization reversal mode evolves into curling [8]. The magnetization reversal mode is determined by the competition between the exchange energy and the demagnetization energy [30]. Generally, the magnetic particles with large size tend towards a reversal in a curling mode, which results in no demagnetization. As the size decreases, the exchange energy density increases, owing to the increase in the relative angle between neighboring magnetic moments. When the size falls into a coherence radius, they tend towards a reversal in coherent-rotation mode. To gain an in-depth fundamental understanding of the magnetic reversal, we systematically investigated the effect of soft magnetic shell thickness on the magnetic structure evolution during demagnetization. The magnetic structures of different soft magnetic shell thickness prior to coercivity of the core/shell cylinder with L/D = 1 are shown in Figure 3. As can be seen, the magnetic reversal with different soft magnetic shell thickness exhibits three different characteristics. For the cases of t = 1 nm and t = 2 nm, the magnetic reversal begins at the soft magnetic shell, and the magnetization reversal mode is coherent rotation. When t = 3 nm, magnetization reversal occurs in a quasi-coherent rotation, where the partial magnetic moments of the fringe are tilted. However, when t reaches above 4 nm, the magnetization reversal mode evolves into curling [8]. The magnetization reversal mode is determined by the competition between the exchange energy and the demagnetization energy [30]. Generally, the magnetic particles with large size tend towards a reversal in a curling mode, which results in no demagnetization. As the size decreases, the exchange energy density increases, owing to the increase in the relative angle between neighboring magnetic moments. When the size falls into a coherence radius, they tend towards a reversal in coherent-rotation mode.

Effect of Shape Anisotropy on the Magnetic Properties
Intuitively, the shape anisotropy can be enhanced by elongating the cylinder. To ex amine the effect of the shape anisotropy of the hard magnetic phase on magnetic proper ties, we calculated the demagnetization curves of the core/shell cylinders with differen aspect ratios. The representative demagnetization curves for the case of t = 1 nm, t = 3 nm and t = 5 nm are shown in Figure 4a-c, respectively. The coercivity Hci and maximum magnetic energy product (BH)max deduced from Figure 4a-c are shown as a function o the aspect ratio in Figure 4d-f, respectively. When the thickness of the soft magnetic shel is only 1 nm, as the aspect ratio increases, the coercivity first increases rapidly and then tends to saturate when the aspect ratio reaches above 3 (Figure 4d), and the demagnetiza tion curves with L/D > 3 are in general overlapped (Figure 4a), suggesting that further enhancing the shape anisotropy would make little contribution to the coercivity. Mean while, the increasing aspect ratio of the hard magnetic core leads to a decrease in Vsoft and thus decreases the Ms. As a result, the (BH)max decreases monotonously with the aspec ratio increasing (Figure 4d). For t = 1 nm, the highest (BH)max (0.62 MJ m −3 , or 78 MGOe and coercivity (3.16 T) can be achieved with aspect ratio L/D = 1 and L/D = 3, respectively In contrast, when t = 3 nm, as the aspect ratio increases, the coercivity first increases and then levels off, and the (BH)max increases monotonously with the aspect ratio increasing (Figure 4e). This indicates that shape anisotropy can provide additional coercivity to com pensate for the coercivity decreasing caused by soft magnetic shell thickness increasing and thus benefit the optimization of the magnetic energy product. The highest (BH)ma

Effect of Shape Anisotropy on the Magnetic Properties
Intuitively, the shape anisotropy can be enhanced by elongating the cylinder. To examine the effect of the shape anisotropy of the hard magnetic phase on magnetic properties, we calculated the demagnetization curves of the core/shell cylinders with different aspect ratios. The representative demagnetization curves for the case of t = 1 nm, t = 3 nm, and t = 5 nm are shown in Figure 4a-c, respectively. The coercivity H ci and maximum magnetic energy product (BH) max deduced from Figure 4a-c are shown as a function of the aspect ratio in Figure 4d-f, respectively. When the thickness of the soft magnetic shell is only 1 nm, as the aspect ratio increases, the coercivity first increases rapidly and then tends to saturate when the aspect ratio reaches above 3 (Figure 4d), and the demagnetization curves with L/D > 3 are in general overlapped (Figure 4a), suggesting that further enhancing the shape anisotropy would make little contribution to the coercivity. Meanwhile, the increasing aspect ratio of the hard magnetic core leads to a decrease in V soft and thus decreases the M s . As a result, the (BH) max decreases monotonously with the aspect ratio increasing (Figure 4d). For t = 1 nm, the highest (BH) max (0.62 MJ m −3 , or 78 MGOe) and coercivity (3.16 T) can be achieved with aspect ratio L/D = 1 and L/D = 3, respectively. In contrast, when t = 3 nm, as the aspect ratio increases, the coercivity first increases and then levels off, and the (BH) max increases monotonously with the aspect ratio increasing (Figure 4e). This indicates that shape anisotropy can provide additional coercivity to compensate for the coercivity decreasing caused by soft magnetic shell thickness increasing and thus benefit the optimization of the magnetic energy product. The highest (BH) max (0.83 MJ m −3 , or 104 MGOe) and coercivity (1.28 T) can be both achieved at the aspect ratio of L/D = 10. For the case of t = 5 nm, owing to the too-large volume fraction of the soft magnetic phase, the collapse of magnetization occurred in the core/shell cylinders with a low aspect ratio L/D ≤ 2 during demagnetization ( Figure 4c). As seen from Figure 4f, with the aspect ratio increasing, the (BH) max increases monotonously, but the coercivity first increases and then decreases. The partial values of the magnetic energy product and coercivity are as follows:  Figure 4f. Here, the squareness is defined as the ratio of the area below the demagnetization curve to the product of remanence M r and coercivity H ci in the second quadrant, which is also called rectangularity [31]. Clearly, the squareness increases monotonously with the aspect ratio increasing. Based on the Stoner-Wohlfarth model, the enhancement of (BH) max may be due to the improvement in squareness of the demagnetization curves [31].
12, x FOR PEER REVIEW 6 of 12 ratio of the area below the demagnetization curve to the product of remanence Mr and coercivity Hci in the second quadrant, which is also called rectangularity [31]. Clearly, the squareness increases monotonously with the aspect ratio increasing. Based on the Stoner-Wohlfarth model, the enhancement of (BH)max may be due to the improvement in squareness of the demagnetization curves [31]. In order to understand the details of magnetization reversal mechanisms in core/shell cylinders associated with different aspect ratios, the revolution of the magnetic structures prior to Hci was monitored. The representative vector maps of magnetic structures in the core/shell cylinders with t = 3 nm are shown in Figure 5. The z-component magnetization in Figure 5 is normalized with respect to its saturation value [32]. Clearly, when the aspect In order to understand the details of magnetization reversal mechanisms in core/shell cylinders associated with different aspect ratios, the revolution of the magnetic structures Nanomaterials 2022, 12, 1261 7 of 12 prior to H ci was monitored. The representative vector maps of magnetic structures in the core/shell cylinders with t = 3 nm are shown in Figure 5. The z-component magnetization in Figure 5 is normalized with respect to its saturation value [32]. Clearly, when the aspect ratio is below the cross-over (L/D < 3), the magnetocrystalline may be dominant over the shape anisotropy, and the reversal of the magnetic moments along the axis of the cylinder occurs in a quasi-coherent way. However, when the aspect ratio increases to 3, the magnetic reversal starts at the two edges of the cylinder, followed by propagation of two domain walls (head-to-head and tail-to-tail) towards the center [33]. This indicates that at this aspect ratio, the shape anisotropy parallel to the axis of the cylinders dominates, which leads to the increase in coercivity.
anomaterials 2022, 12, x FOR PEER REVIEW this aspect ratio, the shape anisotropy parallel to the axis of the cylinde which leads to the increase in coercivity. Figure 5. Vector maps of magnetization structure close to the coercivity of core/shel different aspect ratios L/D. The thickness of the soft magnetic shell is 3 nm. The co and color map indicating the direction of the magnetization are shown next to the i As mentioned above, the coercivity increases with the aspect ratio, fades away when the aspect ratio reaches a critical value. The static magne the ellipsoid sets the upper limit of the shape anisotropy coefficient, which stand the correlation between the aspect ratio and the shape anisotropy, ev ity of the present core/shell cylinders. If no other anisotropy is considered shape anisotropy energy is the total demagnetization energy. In the config ellipsoid ( c > a = b, where a, b, and c are the ellipsoid semi-axes), the coerc proportioned to the Na-Nc, where Nc and Na are the demagnetization factor ing [28]: As mentioned above, the coercivity increases with the aspect ratio, but this effect fades away when the aspect ratio reaches a critical value. The static magnetic behavior of the ellipsoid sets the upper limit of the shape anisotropy coefficient, which helps understand the correlation between the aspect ratio and the shape anisotropy, even the coercivity of the present core/shell cylinders. If no other anisotropy is considered, the magnetic shape anisotropy energy is the total demagnetization energy. In the configuration of the ellipsoid (c > a = b, where a, b, and c are the ellipsoid semi-axes), the coercivity is simply proportioned to the N a -N c , where N c and N a are the demagnetization factors, and satisfying [28]: where m = c/a. The calculated demagnetization factor N a , N c, and N a -N c of the ellipsoid are shown as a function of its aspect ratio in Figure 6. Clearly, the N a -N c first increases rapidly and then increases relatively slowly with the aspect ratio increasing. When the aspect ratio increases approaching 10, further increasing aspect ratio is useless in terms of shape anisotropy optimization [30,33], not to mention the coercivity. As demonstrated from the above results, the shape anisotropy contri magnetic properties in the exchange-coupled magnets are closely related to of the soft magnetic shell. When the soft magnetic shell thickness is very thin properties are dominated by the hard-soft exchange coupling effects, and tions of shape anisotropy would be quite limited. Nevertheless, it is also hig ing to control and manipulate such core/shell exchange coupling magnets, taneously have a large volume fraction of soft magnetic phase and a thin thickness. In the case of a crude core/shell structure with a much thicker shell, the contribution of shape anisotropy emerges and significantly impro netic properties, including the coercivity and magnetic energy product. Co flexibility or feasibility of scalable fabrication, utilizing shape anisotropy w fective method to improve the magnetic performance of hard-soft exch magnets.

The Interaction Effect of Packed Cylindrical Core/Shell Structures
As discussed above, exchange coupling and shape anisotropy has be cally studied for isolated core/shell cylinders. Nevertheless, the magnetic pro entire magnet are determined by those of the unit core/shell cylinder in iso interaction effects due to packing. To examine the effect of packing, we con hexagonal array of the core/shell cylinders, as illustrated in Figure 1b. Arra inders with geometries of L/D = 10 are considered, and no periodic bounda were applied in order to avoid artificial interference effects [34]. This simulat extension of a simple hexagonal array of the core/shell cylinders, which wi As demonstrated from the above results, the shape anisotropy contributions to the magnetic properties in the exchange-coupled magnets are closely related to the thickness of the soft magnetic shell. When the soft magnetic shell thickness is very thin, the magnetic properties are dominated by the hard-soft exchange coupling effects, and the contributions of shape anisotropy would be quite limited. Nevertheless, it is also highly challenging to control and manipulate such core/shell exchange coupling magnets, which simultaneously have a large volume fraction of soft magnetic phase and a thin soft magnetic thickness. In the case of a crude core/shell structure with a much thicker soft magnetic shell, the contribution of shape anisotropy emerges and significantly improves the magnetic properties, including the coercivity and magnetic energy product. Considering the flexibility or feasibility of scalable fabrication, utilizing shape anisotropy would be an effective method to improve the magnetic performance of hard-soft exchange-coupled magnets.

The Interaction Effect of Packed Cylindrical Core/Shell Structures
As discussed above, exchange coupling and shape anisotropy has been systematically studied for isolated core/shell cylinders. Nevertheless, the magnetic properties of the entire magnet are determined by those of the unit core/shell cylinder in isolation and the interaction effects due to packing. To examine the effect of packing, we consider a simple hexagonal array of the core/shell cylinders, as illustrated in Figure 1b. Arrays of unit cylinders with geometries of L/D = 10 are considered, and no periodic boundary conditions were applied in order to avoid artificial interference effects [34]. This simulation is the first extension of a simple hexagonal array of the core/shell cylinders, which will be extended to larger arrays in the future. The demagnetization curves of the hexagonal cylindrical array with different thicknesses of the soft magnetic shell are shown in Figure 7a. The coercivity and maximum magnetic energy product deduced from Figure 7a are shown in Figure 7b. It should be noted that the volume fraction of the soft magnetic phase remains unchanged during packing. Therefore, the assembly of hard core/soft shell cylinders should not experience a reduction in saturation magnetization compared to the isolated unit. However, as seen from Figure 7a, the coercivity of the array is always lower than that in an isolated core/shell cylinder. As a result, the packing effect deteriorates the maximum magnetic energy product, especially when t is larger than 2 nm, as shown in Figure 7b. A similar packing effect (or proximity effect) has also been reported for some magnetic nanowires [6]. The representative magnetic structure evolution of the hexagonal close-packed array of core/shell cylinders with t = 2 nm during demagnetization is shown in Figure 8. In the hexagonal close-packed array, when t ≥ 2 nm, the reversal mechanism of the core/shell cylinder surrounded by an array differs from the magnetic reversal of an isolated cylinder. As seen from Figure 8, the nucleation starts at the soft magnetic shell of the core/shell cylinders in the fringe of the array. Then, a vortex structure is formed, and the center cylinder in the array acts as the core of the vortex. It should be noted that a coherent-rotation mode rather than the curling mode was calculated in the isolated core/shell cylinder with soft magnetic shell thickness of 2 nm. This would be owing to the fact that the single isolated core/shell cylinder only experiences a homogenous external field, while the core/shell cylinder unit in an array is simultaneously affected by the external field and the stray filed of the surrounded cylinders [35], which allows for building a vortex with reduced stray fields in the space between the cylinders [34]. The interaction of each cylinder with the stray field produced by the array would provide an additional antiferromagnetic coupling between the neighboring cylinders, which leads to a decrease in coercivity [36]. Therefore, in order to improve the coercivity of such hard core/soft shell permanent magnets, the units should be separated in a nonmagnetic matric to suppress the proximity effect [12]. The representative magnetic structure evolution of the hexagonal close-packed array of core/shell cylinders with t = 2 nm during demagnetization is shown in Figure 8. In the hexagonal close-packed array, when t ≥ 2 nm, the reversal mechanism of the core/shell cylinder surrounded by an array differs from the magnetic reversal of an isolated cylinder. As seen from Figure 8, the nucleation starts at the soft magnetic shell of the core/shell cylinders in the fringe of the array. Then, a vortex structure is formed, and the center cylinder in the array acts as the core of the vortex. It should be noted that a coherentrotation mode rather than the curling mode was calculated in the isolated core/shell cylinder with soft magnetic shell thickness of 2 nm. This would be owing to the fact that the single isolated core/shell cylinder only experiences a homogenous external field, while the core/shell cylinder unit in an array is simultaneously affected by the external field and the stray filed of the surrounded cylinders [35], which allows for building a vortex with reduced stray fields in the space between the cylinders [34]. The interaction of each cylinder with the stray field produced by the array would provide an additional antiferromagnetic coupling between the neighboring cylinders, which leads to a decrease in coercivity [36]. Therefore, in order to improve the coercivity of such hard core/soft shell permanent magnets, the units should be separated in a nonmagnetic matric to suppress the proximity effect [12].

Conclusions
In summary, we investigated the geometry, including the thickness of the soft magnetic shell and the aspect ratio of the hard magnetic core, and the packing effect of the cylindrical core/shell hard-soft exchange-coupled magnets by micromagnetic simulations. The underlying magnetization reversal mechanisms were revealed. The following main conclusions can be drawn: 1. As the thickness of soft magnetic shell t increases, the saturate magnetization increases monotonously while the coercivity decreases monotonously, and the magnetization reversal changes from coherent rotation mode (t ≤ 2 nm) to quasi-coherent rotation mode (t = 3 nm), and finally to curling mode (t ≥ 4 nm). 2. When the soft magnetic shell is thin (t = 1 nm), as the aspect ratio increases, the (BH)max decreases monotonously, and there is a trade-off between the Hci and (BH)max. The maximum coercivity of 3.16 T can be achieved at the aspect ratio of 3, and further enhancing the shape anisotropy made little contribution to the coercivity. 3. When the soft magnetic shell is relatively thicker (t ≥ 3 nm), as the aspect ratio increases, the (BH)max increases monotonously. For t = 3 nm, the highest (BH)max (0.83 MJ m −3 , or 104 MGOe) and coercivity (1.28 T) can be both achieved at the aspect ratio L/D = 10. 4. The packing effect could lead to a decrease in coercivity and (BH)max in the hard core/soft shell cylinder array, compared to the isolated one. The magnetization reversal of the cylindrical array occurs in an extraordinary curling mode when t ≥ 2 nm, which is different from the isolated one.

Conclusions
In summary, we investigated the geometry, including the thickness of the soft magnetic shell and the aspect ratio of the hard magnetic core, and the packing effect of the cylindrical core/shell hard-soft exchange-coupled magnets by micromagnetic simulations. The underlying magnetization reversal mechanisms were revealed. The following main conclusions can be drawn:

1.
As the thickness of soft magnetic shell t increases, the saturate magnetization increases monotonously while the coercivity decreases monotonously, and the magnetization reversal changes from coherent rotation mode (t ≤ 2 nm) to quasi-coherent rotation mode (t = 3 nm), and finally to curling mode (t ≥ 4 nm).

2.
When the soft magnetic shell is thin (t = 1 nm), as the aspect ratio increases, the (BH) max decreases monotonously, and there is a trade-off between the H ci and (BH) max . The maximum coercivity of 3.16 T can be achieved at the aspect ratio of 3, and further enhancing the shape anisotropy made little contribution to the coercivity.

3.
When the soft magnetic shell is relatively thicker (t ≥ 3 nm), as the aspect ratio increases, the (BH) max increases monotonously. For t = 3 nm, the highest (BH) max (0.83 MJ m −3 , or 104 MGOe) and coercivity (1.28 T) can be both achieved at the aspect ratio L/D = 10.

4.
The packing effect could lead to a decrease in coercivity and (BH) max in the hard core/soft shell cylinder array, compared to the isolated one. The magnetization rever-sal of the cylindrical array occurs in an extraordinary curling mode when t ≥ 2 nm, which is different from the isolated one.
This work could provide important guidelines for developing promising hard-soft exchange-coupled permanent magnets with high magnetic properties.