No Arabic abstract
Computational design of more efficient rare earth/transition metal (RE-TM) permanent magnets requires accurately calculating the magnetocrystalline anisotropy (MCA) at finite temperature, since this property places an upper bound on the coercivity. Here, we present a first-principles methodology to calculate the MCA of RE-TM magnets which fully accounts for the effects of temperature on the underlying electrons. The itinerant electron TM magnetism is described within the disordered local moment picture, and the localized RE-4f magnetism is described within crystal field theory. We use our model, which is free of adjustable parameters, to calculate the MCA of the RCo$_5$ (R=Y, La-Gd) magnet family for temperatures 0--600 K. We correctly find a huge uniaxial anisotropy for SmCo$_5$ (21.3 MJm$^{-3}$ at 300 K) and two finite temperature spin reorientation transitions for NdCo$_5$. The calculations also demonstrate dramatic valency effects in CeCo$_5$ and PrCo$_5$. Our calculations provide quantitative, first-principles insight into several decades of RE-TM experimental studies.
The development of permanent magnets containing less or no rare-earth elements is linked to profound knowledge of the coercivity mechanism. Prerequisites for a promising permanent magnet material are a high spontaneous magnetization and a sufficiently high magnetic anisotropy. In addition to the intrinsic magnetic properties the microstructure of the magnet plays a significant role in establishing coercivity. The influence of the microstructure on coercivity, remanence, and energy density product can be understood by {using} micromagnetic simulations. With advances in computer hardware and numerical methods, hysteresis curves of magnets can be computed quickly so that the simulations can readily provide guidance for the development of permanent magnets. The potential of rare-earth reduced and free permanent magnets is investigated using micromagnetic simulations. The results show excellent hard magnetic properties can be achieved in grain boundary engineered NdFeB, rare-earth magnets with a ThMn12 structure, Co-based nano-wires, and L10-FeNi provided that the magnets microstructure is optimized.
Multiscale simulation is a key research tool for the quest for new permanent magnets. Starting with first principles methods, a sequence of simulation methods can be applied to calculate the maximum possible coercive field and expected energy density product of a magnet made from a novel magnetic material composition. Fe-rich magnetic phases suitable for permanent magnets can be found by adaptive genetic algorithms. The intrinsic properties computed by ab initio simulations are used as input for micromagnetic simulations of the hysteresis properties of permanent magnets with realistic structure. Using machine learning techniques, the magnets structure can be optimized so that the upper limits for coercivity and energy density product for a given phase can be estimated. Structure property relations of synthetic permanent magnets were computed for several candidate hard magnetic phases. The following pairs (coercive field (T), energy density product (kJ/m3)) were obtained for Fe3Sn0.75Sb0.25: (0.49, 290), L10 FeNi: (1, 400), CoFe6Ta: (0.87, 425), and MnAl: (0.53, 80).
Magnetocrystalline anisotropy, the microscopic origin of permanent magnetism, is often explained in terms of ferromagnets. However, the best performing permanent magnets based on rare earths and transition metals (RE-TM) are in fact ferrimagnets, consisting of a number of magnetic sublattices. Here we show how a naive calculation of the magnetocrystalline anisotropy of the classic RE-TM ferrimagnet GdCo$_5$ gives numbers which are too large at 0 K and exhibit the wrong temperature dependence. We solve this problem by introducing a first-principles approach to calculate temperature-dependent magnetization vs. field (FPMVB) curves, mirroring the experiments actually used to determine the anisotropy. We pair our calculations with measurements on a recently-grown single crystal of GdCo$_5$, and find excellent agreement. The FPMVB approach demonstrates a new level of sophistication in the use of first-principles calculations to understand RE-TM magnets.
We present a computational study of the compound Y(Co$_{1-x-y}$Fe$_x$Cu$_y$)$_5$ for 0 $leq x,y leq 0.2$. This compound was chosen as a prototype for investigating the cell boundary phase believed to play a key role in establishing the high coercivity of commercial Sm-Co 2:17 magnets. Using density-functional theory, we have calculated the magnetization and magnetocrystalline anisotropy at zero temperature for a range of compositions, modeling the doped compounds within the coherent potential approximation. We have also performed finite temperature calculations for YCo$_5$, Y(Co$_{0.838}$Cu$_{0.162}$)$_5$ and Y(Co$_{0.838}$Fe$_{0.081}$Cu$_{0.081}$)$_5$ within the disordered local moment picture. Our calculations find that substituting Co with small amounts of either Fe or Cu boosts the magnetocrystalline anisotropy $K$, but the change in $K$ depends strongly on the location of the dopants. Furthermore, the calculations do not show a particularly large difference between the magnetic properties of Cu-rich Y(Co$_{0.838}$Cu$_{0.162}$)$_5$ and equal Fe-Cu Y(Co$_{0.838}$Fe$_{0.081}$Cu$_{0.081}$)$_5$, despite these two compositions showing different coercivity behavior when found in the cell boundary phase of 2:17 magnets. Our study lays the groundwork for studying the rare earth contribution to the anisotropy of Sm(Co$_{1-x-y}$Fe$_x$Cu$_y$)$_5$, and also shows how a small amount of transition metal substitution can boost the anisotropy field of YCo$_5$.
As Eu and Gd are zero-orbital-momentum ($L=0$) rare-earth atoms, their crystalline intermetallic alloys illustrate the connection between electron bands and magnetic anisotropy. Here we find out-of-plane anisotropy in 2D atom-thick EuAu$_2$ by X-ray magnetic circular dichroism. Angle-resolved photoemission and density-functional theory prove that this is due to strong $f-d$ band hybridization and Eu$^{2+}$ valence. In contrast, the in-plane anisotropy of the structurally-equivalent GdAu$_2$ is ruled by spin-orbit-split $d$-bands, notably Weyl nodal lines, occupied in the Gd$^{3+}$ state. Irrespectively of $L$, we predict a similar itinerant electron contribution to the anisotropy of analogous compounds.