Power law analysis for temperature dependence of magnetocrystalline anisotropy constants of Nd$_2$Fe$_{14}$B magnets

Phenomenological analysis for the temperature dependence of the magnetocrystalline anisotropy (MA) in rare earth magnets is presented. We define phenomenological power laws applicable to compound magnets using the Zener theory, apply these laws to the magnetocrystalline anisotropy constants (MACs) of Nd$_2$Fe$_{14}$B magnets. The results indicate that the MACs obey the power law well, and a general understanding for the temperature-dependent MA in rare earth magnets is obtained through the analysis. Furthermore, to examine the validity of the power law, we discuss the temperature dependence of the MACs in Dy$_2$Fe$_{14}$B and Y$_2$Fe$_{14}$B magnets as examples wherein it is difficult to interpret the MA using the power law.


I. INTRODUCTION
The complex behaviors of magnetic anisotropy (MA) in rare-earth permanent magnets have attracted considerable research attention, since such magnets are widely used as high-performance magnets that combine high magnetocrystalline anisotropy with reasonable magnetization and Curie temperature. [1][2][3][4] Typical examples of such magnets include R 2 Fe 14 B magnets (where R denotes a rare-earth element), which form the focus of the present study. The MA of such magnets strongly depends on temperature and is strongly influenced by the choice of R. 1, [5][6][7] From previous studies, it has been found that the localized-moment picture is suitable to understand the magnetism of R 2 Fe 14 B magnets. For example, Yamada et al. showed through their systematic analysis of the magnetization process for the R 2 Fe 14 B series that parameters representing crystalline electric fields (CEFs) and effective exchange fields (EXFs) at R sites do not depend significantly on R. 8 Furthermore, we have previously confirmed that an effective Hamiltonian using the parameters suitably describes the experimental temperature dependence of the magnetocrystalline anisotropy constants (MACs) of Nd 2 Fe 14 B magnets. 9,10 However, there remains a gap between our microscopic picture and the macroscopic aspects of MA since previous theoretical studies on the temperature dependence of MA for R 2 Fe 14 B magnets are based on complete diagonalization of the Hamiltonian 9,10 or numerical methods such as the Monte Carlo method. 4,11,12 An effective means of bridging this gap is to develop a power-law theory 5,6,[13][14][15][16][17][18] to establish the relationship between MACs and magnetization. The history of investigations on MA based on the power-law theory has been reviewed by Callen and Callen. 16 In their review article, they provided a quantum-mechanical foundation to the power-law theory and extended the foundation to more general cases such as the interaction between two moments. In deriving these power laws, the important assumptions are that a local magnetic moment picture is valid and that the magnetic structure is spatially uniform as assumed in Zener's original work 14 (in addition, a low temperature assumption is required from a mean field perspective 16,19 ). However, the total magnetic moment of R 2 Fe 14 B magnets is the sum of the R moments and the Fe moments; thus, in the case that the directions of both moments are mutually non-collinear, the temperature dependence of the MACs will not obey the power law, and we can no longer regard the magnetic structure as a uniform one. In this regard, Dy 2 Fe 14 B magnets are typical examples exhibiting high non-collinearity because of the antiferromagnetic coupling existing between the Dy and Fe moments, 1,18 and further, a wide plateau has been found in the temperature-dependence curve of their first-order MAC in the low-temperature range. 20,21 In contrast, in Nd 2 Fe 14 B magnets, it is possible to assume that the Nd moments are collinear with the Fe moments. 18,22 In this article, we aim to provide a simple and general understanding of the temperature dependence of the MA of rare-earth permanent magnets, considering Nd 2 Fe 14 B magnets as an example, using the power-law theory for MACs as developed by Zener. 14 Here, we note that the original theory targets single-element magnets rather than compounds magnets; therefore, it is required to suitably correct the power law to apply it to Nd 2 Fe 14 B magnets. In addition, we emphasize that the presence of a large higher-order MAC at the zero temperature considerably effects on the temperature dependence of the lower-order MACs than it; the Nd 2 Fe 14 B magnet is in this case. Subsequently, we first present an explicit form of the extended power law picture within the framework of the Zener theory; next, we prove the validity of the extended power law and discuss the limitations of our approach. In addition, we briefly discuss the applicability of this law to Y 2 Fe 14 B and Dy 2 Fe 14 B magnets. denotes the mth-order MAC at temperature T , and Θ and Φ denote the zenithal and azimuthal angles of the total magnetization direction, respectively. 1-3 This expansion is an infinite series in a purely mathematical analysis, 23 but a finite-order form is usually assumed for practical purposes. In this study, we consider the thirdorder form as follows: where the angle-independent term, K 0 (T ), is omitted. In general, { K ( ) n m (T ) } includes not only the R-and the Fe-sublattice contributions but also the cross effect between them. However, if the directions of the R and Fe magnetizations are approximately collinear for arbitrary values of Θ and Φ, then one can regard Θ and Φ as the angles of the R or/and Fe magnetizations, and therefore, we can divide the total MA into the R-and Fe-sublattice contributions as Here, by assuming that the localized R moments are established by an effective exchange field induced by the Fe magnetization and that the origin of the R anisotropy is local at each R ion site, K R( ) n m (T ) can be phenomenologically described by means of the Zener theory. Namely, when the free-energy density of the R sublattice is expressed as in terms of a phenomenological real parameter A q k (T ) and the tesseral harmonics Z q k (Θ, Φ), the Zener theory 14,15 yields the Akulov-Zener power law, (T ) } by another form corresponding to Eq. (1): and thus, we obtain the below equations of interest: in which the first term represents the Akulov-Zener-Callen-Callen power law (or today simply called the Callen-Callen law), 13,14,16 and the second and/or the third terms represent the corrections to the Callen-Callen law due to the presence of the higher-order MAC (s) at the zero temperature.
Here, we should notice a connection between Zener's phenomenological theory and a microscopic theory. The Zener theory has been discussed well by using a mean field theory (MFT), 16,19,[24][25][26] and the present phenomenological expressions (4) can be supported within a linear theory for CEFs. Yamada et al. explicitly represented K R( ) n m (0) in terms of the CEF parameters by using the perturbative theory with respect to the CEF. 8 Recently, on the basis of the expression, an MA has been discussed from first principles 4,27-35 although in general, the quantitative validity of the perturbative treatment for the MA of rare-earth magnets is not trivial. 8 In finite temperatures, Kuz'min represented K R( ) n m (T ) in terms of the CEF and the EXF within the same framework, where the temperature dependence of the MACs is related to Kuz'min's generalized Brillouin functions (GBFs), 24,26,36,37 and furthermore, Magnani et al. investigated the effect of a finite spin-orbit interaction on K R( ) n m (T ). 38 As pointed out by Keffer, 19 the temperature dependence of MA predicted by the MFT tends to the Akulov-Zener power law in a low temperature range, and in fact, it has been demonstrated that the GBFs and the power law are in quantitative agreement in the low temperature range. 18,26 Although the MFT does not support the power law quantitatively in a high temperature range, it is worth enough to employ the power law in the whole temperature range if limiting the discussion to a qualitative level.
To understand the temperature dependence of the MA represented by the present power law (4), let us consider a simple case wherein A schematic of µ R (T ) and the resultant µ R (T ) 3 , µ R (T ) 10 , and 8 7 µ R (T ) 3 − µ R (T ) 10 are shown in Fig.  1. In the case that K R 2 (0)/K R 1 (0) 1, the correction resulting from the presence of K R 2 (0) in Eq. (5a) is small, and we confirm the applicability of the Callen-Callen law 3 . An interesting feature is observed in the case where K R 2 (0) is comparable to K R 1 (0). In this case, the correction to K R 1 (T ) from K R 2 (0) immediately occurs on increasing the temperature from absolute zero, as indicated by the dashed line in Fig. 1 thus, the Callen-Callen law for K R 1 (T ) is no longer valid in the low-temperature range. In the high-temperature range satisfying µ R (T ) 3 µ R (T ) 10 , we can determine the power law as K R 1 (T ) K R 1 (0) + 8 7 K R 2 (0) µ R (T ) 3 again; however, the proportional factor differs from that of the Callen-Callen law. Figure 2 illustrates the possible forms of K R 1 (T ) for a fixed K R 2 (T ) in Eqs. (5), wherein we fixed K R 2 (T ) as K R 2 (0) = 20 MJ/m 3 (dashed line) and determined K R 1 (T ), indicated by solid lines, for cases where K R 1 (0) = K R 2 (0)/2, K R 1 (0) = −K R 2 (0)/2, and K R 1 (0) = −8K R 2 (0)/7. From the figure, we can observe a uniaxial magnetic anisotropy over the whole temperature range in the case of K R 1 (0) > 0 and an in-plane magnetic anisotropy in the case of K R 1 (0) < −8K R 2 (0)/7. However, in the case where K R 1 (0) takes on an intermediate value with −8K R 2 (0)/7 < K R 1 (0) < 0, we can observe that the MA changes from in-plane to uniaxial at a particular temperature; that is, spin reorientation transition (SRT) occurs at the temperature indicated by the solid line corresponding to K R 1 (0) = −K R 2 (0)/2 in Fig. 2. Therefore, we can immediately predict from the values of K R 1 (0) and K R 2 (0) as to whether the SRT occurs, although the actual SRT temperature cannot be determined because it depends on µ R (T ). As the form of µ R (T ) is approximately the same for rare-earth permanent magnets, 20 the above analysis is valid for magnets in which non-collinearity is negligible.
The application of the Zener theory to the Fe anisotropy (i.e., K Fe( ) n m (T ) in Eq. (2)) requires caution, because it is not trivial that a local-moment picture for the Fe sublattice is appropriate. This issue is referred in a discussion on the anisotropy of Y 2 Fe 14 B magnets in the next section.

III. TEMPERATURE DEPENDENCE OF MACS IN R2Fe14B MAGNETS
Let us apply the power law (4) to an experiment to study the temperature dependence of the MACs of Nd 2 Fe 14 B magnets. One way to estimate µ R (T ) of an R 2 Fe 14 B magnet experimentally is to subtract the magnetization of Y 2 Fe 14 B magnets from that of the R 2 Fe 14 B magnet. In this regard, Hirosawa et al. have shown that µ R (T ) estimated by this method is reproduced suitably within the mean field theory. 20 Although one can directly use the experimental µ R (T ), here, we use the theoretical µ R (T ) based on the mean field theory as follows. Assuming an effective exchange interaction between the total angular momentum of the R ions and the Fe spins to plot the µ R (T ) curve, we have the conventional result: 20,26 g LSJ := 1 + J(J + 1) + S(S + 1) − L(L + 1) 2J(J + 1) , where J denotes the quantum number with respect to the total angular momentum of the R ion, H R the effective exchange field, k B the Boltzmann constant, and T C the Curie temperature of the R 2 Fe 14 B magnet. The temperature dependent Fe magnetization is phenomenologically expressed by the Kuz'min formula: 39,40 µ Fe (T ) = 1 − 0.
where, the shape parameters have been determined by fitting with the temperature-dependent magnetization of the Y 2 K. 20,21 The experimental results, µ Fe (T ), and the resultant µ Nd (T ) are shown in Fig. 3. Using the above parameters, we demonstrate comparisons of our extended power law with two experiments. Figure 4 compares the calculated result with the results of the experiment conducted by Durst et al., 42 who determined the first-order and the second-order MACs for Nd 2 Fe 14 B magnets. In our study, the solid lines in the figure were obtained using Eqs. (5) and (6)  mental results up to third-order MACs for a Nd 2 Fe 14 B magnet. 22 With the use of these results, the solid lines in the figure were plotted utilizing Eqs. (4) and (6) where K Nd 1 (0) = −279, K Nd 2 (0) = 813, K Nd 2 (0) = 151, K Nd 3 (0) = −368, and K Nd 3 (0) = −155 (the unit is unknown). Further, in this case, it can be observed that the present power law correlates well with the experimental results. From the results of both experiments, we can conclude that the proposed predictions and experimental results are in good agreement. Here, we emphasize that the common parameters in the two comparisons are used in µ Nd (T ), and particularly that the value of H Nd = 350 K is reasonable because the strength of the EXF strongly influences the plateau of the MACs appearing in the lowtemperature range. 9 In addition, we must consider the non-collinearity between the R and Fe moments. Although the MACs of R =Nd obey the present power law well, this is not the case for magnets with a high noncollinearity. For example, Dy 2 Fe 14 B magnets show a high non-collinearity between Dy and Fe moments. 18 In fact, the wide plateau observed in K 1 (T ) up to T ∼ 0.5T C from T = 0 20,21 cannot be explained using the present power law for H Dy = 145 K, as estimated by Yamada et al. 8 This violation is attributed to the high non-collinearity, not meaning that 145 K is a very small value. Therefore, a special treatment beyond the present one is necessary in such a case. A detailed analysis on Dy 2 Fe 14 B magnets will be presented in a separate paper.
Finally, let us consider an interesting example in the context of the above discussion, wherein we have ignored the contribution to the MA from the Fe sublattice, K Fe( ) n m (T ), which appeared in Eq. (2). This contribution is, often assuming that the Fe sublattice has an intrinsic MA, estimated from the MA measured in M 2 Fe 14 B magnets with no M moment (M =Y, La, Ce, and so on). 1 Here, let us consider the Y 2 Fe 14 B magnet referred to in the above discussion. The history of phenomenological investigations on the temperature-dependence of the MA of Y 2 Fe 14 B magnets was reviewed by Cadogan et al. 43 In a similar manner, we attempted to apply the powerlaw-like expression to the temperature-dependent MACs of Y 2 Fe 14 B magnets. As an analogy to Eq. (4a), we assume the following expression (previously assumed by Carr, 15 the same formulae were used in the study on the MA of hexagonal Co), as the Fe contribution to K 1 (T ) in Y 2 Fe 14 B magnets, where we use the same parameters as the corresponding ones in Fig. 3 for µ Fe (T ); thus, the remaining adjustable parameters are K Fe 1 (0), K Fe 2 (0), and K Fe 3 (0). Fitting the form (10) to the experimental data, 20,21 we obtain K Fe 1 (0) = 0.77 MJ/m 3 , K Fe 2 (0) = 1.21 MJ/m 3 , and K Fe 3 (0) = 0.11 MJ/m 3 . The comparison of the expression (10) with the experimental data is shown in Fig.  6, wherein we note that a good fit is obtained. Here, we should emphasize that assuming the presence of finite K Fe 2 (0) and K Fe 3 (0) leads us to a considerably large K Fe 2 (T ) value (and, a small but not negligible K Fe 3 (T ) value) consistent with the treatment analogy using Eqs. (4b) and (4d). However, the presence of such a higherorder MAC has not thus far been reported, as mentioned in an earlier review. 43 We also plot these ghost MACs in the form of dashed and dotted lines in Fig. 6. This situation is the same as that encountered by Zener in the study on the MA of nickel, 14 which suggests that the assumptions of the Zener theory are invalid for Y 2 Fe 14 B magnets. Therefore, the expression (10) along with Eq. The dashed and dotted lines, respectively, represent the second-order and the third-order MACs that should be present for consistency of the Zener theory.
MnBi. 44,45 IV. SUMMARY In summary, we have shown explicit forms of the power law applicable to rare-earth magnets within the framework of the Zener theory, and we have demonstrated that the power law suitably describes the experimental MACs of Nd 2 Fe 14 B magnets. Thus, we obtained a simple and general understanding of the complex temperature dependence of the MA often observed in rareearth magnets. Furthermore, we considered Dy 2 Fe 14 B and Y 2 Fe 14 B magnets as examples regarding the interpretation of MA when the power law failed. We believe that our findings can further contribute to studies on the magnetocrystalline anisotropy of rare-earth magnets.