5 On the applicability of the formula $ {\cal K}_\mathrm {eff} = {\cal K}_\infty + 6 {\cal K}_s/D$

The results presented above show that in the most general case, studied here, of a multi-spin particle with Néel surface anisotropy, this formula is not applicable, for the following reasons:
(i)
It assumes that the overall anisotropy of the particle remains uniaxial. However, we have shown that the surface anisotropy induces an additional cubic contribution.
(ii)
It assumes that the surface anisotropy always enhances that in the core. In this and the previous sections we saw that both situations can arise.
(iii)
It is implicitly based on the hypothesis that the core and surface anisotropies are additive contributions. As we have seen above for large $ k_s$ the energy barrier indeed can be represented as a sum of the effective cubic and uniaxial anisotropies. However, the cubic anisotropy term is proportional to $ k_s^2$, which is inconsistent with formula (3.5).
(iv)
It assumes linear dependence of energy barriers on the parameter $ 1/D$, or equivalently $ N_s/{\cal N}$.

Consequently, spherical or octahedral particles cannot be described by formula (3.5), since in this case:

(i)
No term linear in $ k_s$ is obtained.
(ii)
No term scales as the ratio of the surface-to-volume number of spins $ N_s/{\cal N}$.
However, in the case of elongated particles with a not too large surface anisotropy, i.e., $ \vert\zeta\vert < 1$, the energy barriers are independent of the effective cubic anisotropy. In this case, for weakly ellipsoidal particles, for example, we may write

$\displaystyle \Delta E_{EOSP} = k_\mathrm{ua}^\mathrm{eff}\approx k_c N_c/{\cal N} + A\vert k_s\vert/{\cal N}^{1/3}$ (52)

where $ A$ is a parameter that depends on the particle's elongation and surface arrangement, and which is positive for $ k_s < 0$ and negative in the opposite case. Hence, the behavior is as predicted by this formula (3.5). The approximately linear behavior in $ N_s/{\cal N}$ was also observed in the case of large surface anisotropy $ \zeta»1$ (see Fig. 3.1). However, in this case at $ {\cal N} \rightarrow \infty$, the ``uniaxial anisotropy term" $ \cal{K}_{\infty}$ is renormalized by the effective cubic anisotropy $ k_\mathrm{ca}^\mathrm{eff} \thicksim k_s^2$. In Fig. 3.8 we plot the energy barriers of small ellipsoidal particles with sc structure, aspect ratio 2:3, and $ k_s < 0$ from Fig. 3.7.
Figure 3.8: Linear fit of energy barriers, versus the ratio of the surface-to-total number of spins $ N_S/{\cal N}$, of ellipsoidal particles with aspect ratio 2:3 with $ k_s = -0.041$ (circles) and $ k_s = -0.1125$ (squares).
\includegraphics[totalheight=0.35\textheight]{barrier_ellips_1.eps}
For such particles, the formula (3.5) should be modified as $ {\cal K}_\mathrm{eff}={\cal K}{\infty} + \vert{\cal K}_s\vert N_s/{\cal N}$. Accordingly, in Fig. 3.8 we plot the energy barrier against $ N_s/{\cal N}$. These data are well fitted to straight lines, especially when small particle sizes are removed. We note that in the case of relatively small surface anisotropy $ k_s = -0.041$ (though 17 times larger than that in the core), the full core anisotropy $ {\cal K}_{\infty} = K_c/v$ ( $ v$ is the atomic volume) can be extracted. However, for the larger surface anisotropy $ k_s = -0.1125$, $ {\cal K}_{\infty}$ is renormalized by the surface contribution, defined by Eq. (2.15). On the other hand, it is not possible to extract the value of $ k_s$, since the exact proportionality coefficient of Eq. (2.17) (the value of A in Eq. 3.10) is dependent on the particles surface arrangement and elongation. The effective anisotropy constant $ {\cal K}_s$ obtained from this fit is much smaller than the input value, namely, for $ k_s/k_c = 45$ we obtain from the fit $ (k_s/k_c)_\mathrm{eff}=4.3$.
Rocio Yanes