Abstract:Density functional theory (DFT), including its extensions designed to treat strongly correlated localized electron systems such as DFT+$U$ and DFT+dynamical mean field theory, has proven exceedingly useful in studying the magnetic properties of solids. However, materials with rare earths ($R$) have remained a notable exception. The most vital rare-earth magnetic properties, such as magnetocrystalline anisotropy (MA), have been notoriously elusive due to the ubiquitous self-interaction error present in nearly all available DFT flavors. In this work, we show explicitly how the orbital dependence of self-interaction error may contradict Hund's rules and plague MA calculations, and how analyzing DFT metastable states that respect Hund's rules can alleviate the problem. We systematically investigate and discuss several rare-earth-containing families, $R$Co$_5$, $R_2$Fe$_{14}$B, $R$Fe$_{12}$, and $R$Mn$_6$Sn$_6$, to benchmark the MA calculations in DFT+$U$. For all compounds we investigated, we found that our methodology reproduces the magnetic easy-axis, easy-plane, and non-trivial easy-cone anisotropies in full agreement with low-temperature experimental measurements. Besides the fully-numerical ab initio approach, we further illustrate an efficient semi-analytical perturbation method that treats the crystal field as a perturbation in the limit of large spin-orbit coupling. This approach evaluates the rare-earth anisotropy by assessing the dependence of crystal-field energy on spin-quantization axis rotation using $4f$ crystal-field levels obtained from non-spin-orbit calculations. Our analytical method provides a quantitative microscopic understanding of the factors that control MA and can be used for predicting new high-MA materials.

What problem does this paper attempt to address?

### Problems the Paper Attempts to Solve This paper aims to address the self-interaction error (SIE) encountered in density functional theory (DFT) when dealing with rare earth (RE) magnetic materials. Specifically, the magnetocrystalline anisotropy (MA) of rare earth elements is difficult to calculate accurately due to strong self-interaction errors. The paper alleviates this issue by analyzing the metastable states of DFT to comply with Hund's rules and systematically studies several types of compounds containing rare earth elements (such as RCo5, R2Fe14B, RFe12, and RMn6Sn6) to verify the effectiveness of the DFT+U method in calculating magnetocrystalline anisotropy. ### Summary of Main Content 1. **Background Introduction**: - Rare earth elements exhibit significant magnetocrystalline anisotropy (MA) due to their strongly localized 4f orbitals and strong spin-orbit coupling (SOC). - DFT and its extended methods (such as DFT+U and DFT+dynamic mean-field theory) are very useful in studying solid-state magnetism, but they have significant self-interaction errors when dealing with rare earth elements, leading to inaccurate MA calculations. 2. **Problem Analysis**: - Self-interaction errors lead to incorrect occupancy states of the 4f orbitals, particularly violating Hund's rules. - By analyzing the metastable states of DFT, states that comply with Hund's rules can be found, thereby improving the calculation of MA. 3. **Research Methods**: - Systematically studied compounds containing rare earth elements such as RCo5, R2Fe14B, RFe12, and RMn6Sn6, using the DFT+U method to calculate their magnetocrystalline anisotropy. - By enforcing Hund's rules, the correct occupancy states of the 4f orbitals are ensured, thereby improving the accuracy of the calculations. 4. **Results and Discussion**: - For all studied compounds, after enforcing Hund's rules, the DFT+U method can reproduce the experimentally measured easy axis, easy plane, and non-trivial easy cone anisotropies. - A highly efficient semi-analytical perturbation method is proposed, treating the crystal field as a perturbation under strong spin-orbit coupling, to evaluate MA. 5. **Conclusion**: - By enforcing Hund's rules, the DFT+U method can accurately describe the magnetocrystalline anisotropy of rare earth elements without the need for introducing self-interaction correction (SIC) and orbital polarization correction (OPC). - This method provides a quantitative microscopic understanding for predicting new high MA materials. ### Summary This paper systematically studies and improves the DFT+U method to solve the computational difficulties caused by self-interaction errors in rare earth magnetic materials, providing a powerful tool for understanding and predicting new high MA materials.

Toward a first-principles theory of rare-earth ions in crystals

Understanding and optimization of hard magnetic compounds from first principles

Theory of Defect-Induced Crystal Field Perturbations in Rare-Earth Magnets

Theory of defect-induced crystal field perturbations in rare earth magnets

First-principles calculations of magnetic properties for analysis of magnetization processes in rare-earth permanent magnets

Effective point-charge analysis of crystal electric fields -- application to rare-earth pyrochlores and tripod kagome magnets R3Mg2Sb3O14

Atomistic study of rare-earth compounds

Magnetostriction and Magnetism of Rare Earth Intermetallic Compounds: First Principle Study

Theory of magnetism for rare-earth magnets on the Shastry-Sutherland lattice with non-Kramers ions

Atomistic Simulation of Rare Earth-Transition Intermetallic Compounds Derived from Cacu5 Type

A Self-Adaptive First-Principles Approach for Magnetic Excited States

Relevance of $4f$-$3d$ exchange to finite-temperature magnetism of rare-earth permanent magnets: an ab-initio-based spin model approach for NdFe$_{12}$N

Calculated Curie temperatures for rare-earth permanent magnets: ab initio inspection on localized magnetic moments in d-electron ferromagnetism

Systematic extraction of crystal electric-field effects and quantum magnetic model parameters in triangular rare-earth magnets

Magnetocrystalline anisotropy and orbital polarization in ferromagnetic transition metals

Temperature-dependent magnetocrystalline anisotropy of rare earth/transition metal permanent magnets from first principles: The light RCo$_5$ (R=Y, La-Gd) intermetallics

Rare-earth/transition-metal magnets at finite temperature: Self-interaction-corrected relativistic density functional theory in the disordered local moment picture

Electronic Structure Calculations of Rare-Earth-Doped Magnesium Oxide Based on Density Functional Theory

Designing rare-earth free permanent magnets in Heusler alloys via interstitial doping

Performance of Exchange-Correlation Approximations to Density-Functional Theory for Rare-earth Oxides

Pitfalls of exchange–correlation functionals in description of magnetism: Cautionary tale of the FeRh alloy