Intrinsically weak magnetic anisotropy of cerium in potential hard-magnetic intermetallics

Cerium-based intermetallics are currently attracting much interest as a possible alternative to existing high-performance magnets containing scarce heavy rare-earth elements. However, the intrinsic magnetic properties of Ce in these systems are poorly understood due to the difficulty of a quantitative description of the Kondo effect, a many-body phenomenon where conduction electrons screen out the Ce-4f moment. Here, we show that the Ce-4f shell in Ce–Fe intermetallics is partially Kondo screened. The Kondo scale is dramatically enhanced by nitrogen interstitials suppressing the Ce-4f contribution to the magnetic anisotropy, in striking contrast to the effect of nitrogenation in isostructural intermetallics containing other rare-earth elements. We determine the full temperature dependence of the Ce-4f single-ion anisotropy and show that even unscreened Ce-4f moments contribute little to the room-temperature intrinsic magnetic hardness. Our study thus establishes fundamental constraints on the potential of cerium-based permanent magnet intermetallics.


INTRODUCTION
The last decades have witnessed a rapidly growing demand for energy-efficient technologies, from optimized power generation and conversion devices to new transportation solutions. A key component of motors, generators, sensors, or actuators are highperformance permanent magnets. Identifying and designing improved materials for permanent magnet applications has thus the potential of huge economical and environmental savings 1,2 . Today, the world market for permanent magnets in the more advanced applications is largely dominated by rare-earth-based materials: the best hard-magnetic materials are based on the Nd-Fe-B "2-14-1" system whose extraordinary hard-magnetic properties are further boosted by partially substituting Nd with heavy rare-earth elements like Dy or Tb. In such rare-earth transition metal intermetallics, the electronic structure of the transition metal sublattice provides a large magnetization and high Curie temperature. The intrinsic magnetic hardness κ is due to the rareearth local magnetic moment, which-due to the large spin-orbit coupling acting on the rare-earth electronic states-effectively converts the crystalline anisotropy into single-ion magnetic anisotropy 3 .
Available resources of heavy rare earths are scarce, thus motivating efforts to replace heavy rare-earth intermetallics with so-called rare-earth balance magnets 2 containing more abundant rare-earth elements, in particular Ce 1,[4][5][6][7] . A significant limitation to the magnetic hardness of Ce-transition metal intermetallics is caused by the well-known tendency of Ce-4f electrons to form heavy-electron itinerant bands, as for instance observed in socalled "heavy-fermion" compounds 8 . Indeed, in the heavy-fermion regime the Ce local magnetic moment is screened by conduction electrons, which couple to the single Ce-4f electron to form a nonmagnetic singlet state. Such a heavy-fermion behavior has been observed in various Ce intermetallics, for example, in the CeFe 2 Laves phase 9 , as well as in the RCo 5 (R = rare earth) and R 2 Fe 14 B families of hard-magnetic intermetallics 10,11 compromising, especially, the utility of Ce 2 Fe 14 B for magnetic applications 12 .
In order to possess a significant single-ion magnetic anisotropy, an obvious conditio sine qua non is that the 4f electron of Ce needs to be localized and its magnetic moment not significantly Kondo-screened. A reliable quantitative assessment of localized vs. heavy-fermion behavior of the Ce-4f electrons is thus key to a successful theoretical search for prospective hard-magnetic cerium intermetallics. Such an assessment represents a formidable theoretical challenge. Some qualitative insights can be obtained from density functional theory (DFT) calculations treating Ce-4f as itinerant bands 13 , or estimating the evolution of the Kondo scale 14 from the DFT electronic structure. However, the formation of a heavy-fermion state in Ce-based compounds is a genuine nonperturbative many-body phenomenon, the quantitative description of which requires methods beyond DFT or DFT + U techniques 15,16 . At present, these limitations seriously hamper the possibilities of computational materials design in the field (see e.g., ref. 17 ).
Here, we address this challenge by developing an ab initio approach to the magnetocrystalline anisotropy of Ce-based intermetallics. We describe the electronic structure of magnetic Ce-transition metal intermetallics, using advanced many-body electronic structure methods that take into account the effect of the Kondo screening. Furthermore, we supplement this firstprinciple framework with a variational approach to evaluate the magnitude of the Ce single-ion magnetic anisotropy explicitly taking into account the impact of the Kondo effect. In contrast to previous zero-temperature DFT calculations [17][18][19] , the present study investigates the full temperature dependence of the Ce magnetic anisotropy. Our theoretical predictions are compared to the experimental magnetic anisotropy that we measure in the same wide temperature range.
We focus on the Ce-Fe "1-12" material system 4,7,17,20,21 , which has recently been under scrutiny as a potential hard magnet 16,22,23 . Our calculations show a significant reduction of the Ce magnetic moment due to the Kondo effect. This Kondo screening is found to be very sensitive to substitutions and small-atom insertions. In particular, interstitial nitrogen is predicted to suppress the Ce contribution to the magnetization and magnetic anisotropy. This behavior is in drastic contrast to the one observed in "1-12" intermetallics with localized rare earths, like Nd, where nitrogenation has been found to significantly enhance the uniaxial magnetic anisotropy 24 . Indeed, the unexpected suppression of the Ce magnetization by N interstitials is absent if the Ce moment is treated as purely local and stems from a large enhancement of the Kondo scale by nitrogenation, as schematically shown in Fig. 1b-d. It is confirmed by our measurements of the anisotropy in polycrystalline samples with and without nitrogenation. Such interstitial atoms effectively controlling the state of the Ce-4f shell could be used for fine-tuning and probing the complex manyelectron physics of Ce. Most importantly, our calculations establish the salient fact that even in the absence of Kondo screening, the small size of the Ce-4f spin moment necessarily leads to rather small values of the magnetic anisotropy at elevated temperatures. Our study, therefore, reveals significant limitations for the Ce contribution to the magnetic hardness at temperatures relevant for applications.

RESULTS
Kondo screening of the Ce magnetic moment in "1-12" systems The "1-12" Ce-Fe intermetallics crystallize in the tetragonal lattice structure shown in Fig. 1a. A partial substitution of Fe with another transition metal (usually Ti) is necessary to stabilize this "1-12" phase in the bulk. Hence, we consider the realistic CeFe 11 Ti composition and, for the sake of comparison, also parent thermodynamically unstable CeFe 12 . For the latter, we assume the same crystal structure. The magnetic properties of "1-12" rareearth transition metal intermetallics can be modified by small-atom interstitials, in particular nitrogen. Nitrogenation has been found to significantly enhance the uniaxial magnetic anisotropy in NdFe 12 N and NdFe 11 TiN (ref. 24 ). Nitrogenated Ce-Fe "1-12" compounds have been patented as prospective high-performance magnets 5 , with N interstitials found to significantly increase their Curie temperature and magnetization 5,20,21 . Therefore, we also consider nitrogenated CeFe 11 Ti with N interstitials in the experimentally determined Wyckoff 2b position as shown in Fig. 1a.
The ordered Ce-4f magnetic moment M z in CeFe 12 , CeFe 11 Ti, and CeFe 11 TiN as a function of temperature T is displayed in Fig. 2a as symbols. This Ce moment is induced by the exchange field B ex generated by the ferromagnetically ordered iron sublattice, as schematically shown in Fig. 1b-d. In our ab initio approach, abbreviated below as DMFT(QMC), we take into account the magnetization of Fe within the local spin-density approximation (LSDA), while many-electron effects on the Ce-4f shell are treated within dynamical mean-field theory (DMFT) in Fig. 1 Crystal structure, Ce magnetism, and Kondo effect in Ce-Fe "1-12". a Tetragonal body-centered crystal structure of CeFe 11 TiN (space group I4/mmm). b Quasi-atomic approximation to the Ce magnetism. The Ce-4f spins are aligned antiparallel with respect to the Fe moments due to the exchange field B ex . The magnitude of the Ce magnetic moment in this quasi-atomic picture is determined by B ex and the crystal field splitting on the Ce-4f shell. c Many-electron effects drastically alter this picture: the Kondo screening (green cloud) reduces the Ce magnetic moment concomitantly lowering the Ce contribution to the magnetic anisotropy. d The Kondo effect is strongly sensitive to the local environment: it is dramatically enhanced by nitrogen interstitials thus almost completely suppressing the Ce-4f contribution to the magnetic anisotropy. Fig. 2 Temperature dependence of Ce-4f magnetism in Ce-Fe "1-12". a Calculated temperature dependence of the Ce-4f magnetic moments M z in CeFe 12 (black squares), CeFe 11 Ti (red circles), and CeFe 11 TiN (blue triangles). At low temperatures the Kondo-screened Ce-4f magnetic moments (symbols) are much smaller than their expected quasi-atomic values (solid lines). The Kondo effect is the strongest in CeFe 11 TiN. b Temperature dependence of the anisotropy field H A (Tesla) in CeFe 11 Ti (in red) and CeFe 11 TiN (in blue). In addition to the theoretical and experimental results obtained in this work, values from (d) ref. 20 , (e) ref. 21 , and (f ) ref. 4 are shown for comparison.
conjunction with a numerically exact quantum Monte-Carlo (QMC) approach 25 . With the Fe ferromagnetism included within LSDA, the exchange field B ex acting on the Ce-4f shell is virtually temperature independent (this approximation is justified since we focus on the temperature range up to 300 K, which is well below the Curie temperatures of 485 and 710 K for CeFe 11 Ti and CeFe 11 TiN, respectively). Hence, the M z vs. T curve reflects the temperature dependence of the Ce-4f magnetic susceptibility. One observes a weak T dependence for all three compounds; it is very different from a Curie-like behavior expected for Ce-4f local moments subjected to an exchange field. The obtained values are well below the saturated Russel-Saunders value M z = g J ⋅ J = 2.14μ B for the Ce f 5/2 manifold. In particular, in CeFe 11 TiN the 4f moment is tiny (0.04 μ B ). A small increase of the Ce-4f magnetic moment with respect to pure CeFe 12 is induced by Ti substitution. Overall, the magnitude and temperature dependence of the Ce-4f magnetic moment in the investigated "1-12" intermetallics correspond to a Kondo-screened Fermi-liquid Ce-4f state.
For comparison, the ordered Ce-4f magnetic moment M at z has also been calculated within a quasi-atomic (Hubbard-I) approximation 26,27 to the Ce-4f quantum impurity problem (solid lines in Fig. 2a). In this approach, abbreviated below as DFT + HubI, the Kondo effect is neglected (as schematically shown in Fig. 1b). The Ce magnetic moment evaluated within DFT + HubI exhibits the expected Curie behavior and is significantly larger than the Kondo-screened one at low temperatures. Above T~400 K, the magnitude of the Ce magnetic moments in CeFe 12 and CeFe 11 Ti agree reasonably well between DFT + HubI and DMFT(QMC), indicating this temperature to be the Kondo scale of these compounds. Within the quasi-atomic picture one would expect the largest M at z for CeFe 11 TiN (blue solid line in Fig. 2a), for which the DMFT(QMC) moment remains, however, heavily screened in the whole temperature range up to 600 K.
Magnetic anisotropy of Ce in the quasi-atomic picture In the quasi-atomic approximation (schematically shown in Fig. 1b), the single-ion anisotropy due to the Ce-4f magnetic moments is determined by B ex stemming from the Fe-sublattice and the crystal field (CF) splitting on the Ce-4f shell 28 . While the magnitude of B ex is similar in all three investigated compounds, the CF splitting turns out to be very different (see the crystal field parameters A 0 2 hr 2 i in Table 1). Of particular interest for the Ce-4f magnetic moment at low temperature is the composition of the ground state levels. The CF ground state in CeFe 11 TiN is well separated from excited states and consists almost exclusively of the ±5=2 j i eigenstates, which have the largest jM at z j among the j = 5/2 manifold. CeFe 11 TiN is thus expected to exhibit a strong easy-axis magnetic anisotropy in the quasi-atomic picture. In CeFe 11 Ti, the contribution of the ±5=2 j i eigenstates to the CF ground state is reduced leading to a smaller easy-axis anisotropy. In CeFe 12 , the CF ground state is ±1=2 j icorresponding to an easyplane anisotropy (see Supplementary Fig. 1 for the full CF levels scheme). The peak of the CeFe 12 magnetic moment at~40 K (Fig. 2a) can be explained by the fact that the applied exchange field oriented along the hard axis of CeFe 12 is not sufficient to saturate the ground state moment. Thus, with increasing temperature excited states with a larger moment along the hard axis become thermally populated before the usual Curie decay of the magnetization sets in.
We extract the values of the CF parameters, as well as B ex from our ab initio DFT + HubI calculations and then evaluate the Ce single-ion anisotropy constants. The resulting constants K at 1 (in K per f.u.) are listed in Table 1. At liquid helium temperature, K at 1 is easy-plane in CeFe 12 , and easy-axis in CeFe 11 Ti and CeFe 11 TiN. For T = 300 K, CeFe 11 TiN exhibits a substantially stronger easy-axis anisotropy as compared to CeFe 11 Ti due to a larger value of its low-rank "20" crystal field parameter (see Supplementary Table I for values). Both compounds exhibit a strong reduction of the anisotropy with increasing temperature. The calculated roomtemperature anisotropy constants K at 1 are in very good agreement with the high-T limit analytical expression from ref. 28 , which reads , where g J = 6/7 is the gyromagnetic ratio for the Ce-4f 1 J = 5/2 shell. Using this analytical expression, one may understand the very rapid suppression of the Ce single-ion anisotropy apparent in Table 1. The Zeeman energy 2J|g J − 1|B ex associated to the exchange field B ex~3 00 T gives a characteristic temperature of 2J|g J − 1|B ex /k B~1 40 K, explaining the rapid suppression of the K at 1 values with temperature and the small value at room temperature. For comparison, we note that the same simple estimate for Sm 3+ , with g J = 2/7 gives a five times higher characteristic temperature (700 K), resulting in a 25 times slower reduction of K at 1 vs. T at high temperatures. Hence, as the small magnitude of the spin represents an intrinsic feature of the Ce 3+ ion, i.e., not sensitive to the crystalline environment, the Ce-4f contribution to the magnetic hardness at room temperature and above is also expected to be intrinsically weak. Figure 2b shows our measured anisotropy fields obtained from hard-axis measurements of the textured samples. Our results are consistent with the measurements at room temperature previously reported in refs. 4,20,21 . The anisotropy field is already relatively small at low temperature, and then exhibits a rather slow decay with temperature.
In order to compare theory and experiment, one needs to add the contribution from the Fe-sublattice to the calculated single-ion K at 1 of Ce. As in this work, we do not aim at calculating the Fesublattice magnetic anisotropy, we extract its value and temperature dependence from measurements performed on YFe 11 Ti (ref. 29 ), using the experimental dependence K Fe 1 ðTÞ ¼ 23:3 À 28T=T c (K/f.u.) measured in ref. 29 together with our measured Curie temperatures. At low T, the iron contribution is much smaller than the Ce one and does not affect the picture qualitatively.
Using the relation H A % 2K tot 1 MS , where H A is the anisotropy field, M S the total magnetic moment of the compound saturated in the direction of H A , and K tot 1 ¼ K Fe 1 þ K at 1 , we calculate the corresponding values of H A in the quasi-atomic approximation listed in Table 2 together with the experimental values. The quasi-atomic approximation predicts a strong anisotropy that decreases quickly with temperature, while experimentally one finds only a moderate decrease with T, in particular in CeFe 11 TiN, which bears little resemblance to the predictions of DFT + HubI. Notice that the DFT calculations of ref. 18 treating Ce-4f as delocalized band-like states predict the zero-temperature H A of CeFe 11 Ti to be reduced by nitrogenation by~0.9 T, in disagreement with our measurements and previous experiments 20, 21 finding the difference in H A between CeFe 11 Ti and CeFe 11 TiN x to be in the range of 2.5-3 T (see Table 2). The low-T anisotropy of CeFe 11 Ti seems to be significantly underestimated by DFT 18,19 . Conversely, by treating Ce-4f as an open core Körner et al. 17 obtained H A = 18 T for CeFe 11 Ti, which is, similarly to our DFT + HubI result, strongly overestimated compared to experimental values. Hence, neither DFT nor the fully localized treatment of 4f is able to capture the impact of interstitials on the Ce magnetic anisotropy. As we show below, the theoretical picture is drastically modified when the impact of Kondo screening on the Ce-4f magnetic anisotropy is taken into account within a full many-body framework.
Ce magnetic anisotropy in the presence of Kondo screening A quantitative evaluation of the magnetic anisotropy of a partially Kondo-screened magnetic moment is a formidable problem which-to the best of our knowledge-has not been addressed even for simple model systems so far. In the context of the present work on complex realistic systems, we have developed a variational approach to meet this challenge. In order to capture the interplay of the Kondo interaction (H K ) with the exchange (H B ) and crystal (H CF ) fields in the total impurity Hamiltonian H = H c + H K + H B + H CF for the Ce-4f shell, we assume the ground state of the Ce-4f shell to be given by the following variational manyelectron wave function: where Ψ S is the singlet formed by the Ce-4f total moment coupled to the conduction-electron sea defined by H c ; the second term Ψ J in Eq. (2) represents the 4f ground state in the absence of Kondo screening, i.e., determined by H B + H CF ; α is a variational coefficient. Only the second term, due to its nonzero angular momentum, couples to the exchange field As shown in the Supplementary Note 3, the lowest-order anisotropy constant obtained with the variational ansatz Eq. (2) then reads where K at 1 is the anisotropy constant in the quasi-atomic approximation (DFT + HubI), as detailed above. The value of α 2 is given in accordance with Eq. (3) by M z ðT ¼ 0Þ=M at z ðT ¼ 0Þ and can be extracted by extrapolating the curves in Fig. 2a to zero temperature. The denominator 1 + δ in Eq. (4), where δ ¼ K at 1 =A, is due to the change of α upon the rotation of B ex from the uniaxial B ! ex jjc to in-plane direction and can be evaluated from the dependence of α 2 vs. B ex ||z, calculated by DMFT(QMC) (see Supplementary Note 3 for details).
By multiplying the quasi-atomic K at 1 for Ce with the corresponding reduction factors due to the Kondo screening (Eq. (4)) and adding the contribution of the Fe-sublattice K Fe 1 ðTÞ, we obtain the total DMFT(QMC) anisotropy field listed in Table 2. Strictly speaking, the above variational treatment is applicable only at zero temperature. However, close to room temperature, the cerium contribution is overruled by the iron one, and the final value becomes relatively insensitive to the precise estimate of the Kondo reduction factor. At low temperature, instead, the Kondo reduction is very significant and the Ce contribution to H A is completely suppressed in the nitrogenated compound, but still relatively important in CeFe 11 Ti. The values of the anisotropy field estimated within DMFT(QMC) are in good agreement with our experimental results and previous values reported in refs. 4,20,21 , see Fig. 2b. In particular, the dramatic reduction of the low-T anisotropy field upon nitrogenation predicted by our calculations (by 45%, i.e., 2.6 T) is in quantitative agreement with the reduction observed in experiment. Our prediction for the low-T H A of CeFe 11 Ti (5.8 T) is in the middle of the measured range (5.1-7.4 T), in contrast to the itinerant 4f-in-valence 18,19 and localized 4f-incore 17 DFT limits, which are not able to capture the intermediate partially Kondo-screened state in this system. The itinerant picture with 4f-in-valence 18 performs, however, reasonably well for the strongly screened Ce-4f states in CeFe 11 TiN.
Electronic structure and Ce-4f hybridization As demonstrated above, the magnetic moment and anisotropy of Ce-4f in the CeFe 11 MX intermetallics is heavily influenced by the Kondo effect. Here, we analyze its imprint on the electronic structure of CeFe 11 Ti(N) and identify the origin of the enhanced Kondo screening in the nitrogenated system.
In Fig. 3a, we show the real part of the DMFT self-energy Re Σ(ω) in CeFe 11 Ti and CeFe 11 TiN, respectively, averaged over the six (degeneracy of the J = 5/2 multiplet) Ce-4f orbital contributions with the largest occupation number. The slope of ReΣ(ω) at ω = 0 is a direct proxy for the strength of electronic correlations and determines the quasiparticle weight: The obtained quasiparticle weight is higher in CeFe 11 TiN, as compared to CeFe 11 Ti, indicating a more delocalized character of the Ce-4f states in the nitrogenated compound. This can also be seen from the low occupation of n = 0.86 electrons on the Ce-4f shell in CeFe 11 TiN, deviating stronger from the localized limit with n = 1 electron as compared to n = 0.94 in CeFe 11 Ti. For a further comparison, in CeFe 12 the Ce-4f occupancy is n = 0.93, a bit lower than in CeFe 11 Ti and thus indicating that the necessary Ti substitution slightly favors the localization of the Ce-4f states (for the spectral function of CeFe 12 , we refer the reader to the Supplementary  Fig. 4). The DMFT hybridization function shown in Fig. 3b quantifies the hybridization of the Ce-4f shell with other states of the system. Especially at low energies in the vicinity of the Fermi level, the hybridization function is much larger in CeFe 11 TiN, indicating an increased mixing of the Ce-4f states with neighboring nitrogen states, e.g., N-2p, which is at the origin of the more delocalized behavior of the Ce-4f electrons. The corresponding DMFT(QMC) spectral functions are displayed in Fig. 3c, d. In the characteristic three-peak structure of the Ce-4f spectral function the Kondo peak at the Fermi level is further split into two sub-peaks due to the spin-orbit splitting. The position of the localized features, i.e., the lower and upper Hubbard bands centered at −2 and 5 eV, respectively, is in good agreement with the photoemission (PES) spectra of CeFe 2 (ref. 9 ; we are not aware of any experimental PES for the CeFe 12 materials family). The more delocalized character of the Ce-4f states in CeFe 11 TiN is also reflected in the spectral function, Fig. 3d: a more prominent Kondo  Fig. 3 of this reference.
A. Galler et al.
peak and less pronounced Hubbard bands signal a more delocalized behavior.

DISCUSSION
This work introduces a DMFT-based approach to the Ce magnetization and magnetic anisotropy in ferromagnetic intermetallics. The intermediate valence behavior of Ce in the investigated intermetallics can be captured neither by pure DFT nor by approaches treating the Ce-4f states as purely localized. Indeed, in our work, we show that a fully localized treatment of the Ce-4f states within the quasi-atomic approximation strongly overestimates the Ce contribution to the anisotropy at low temperatures. The problem at hand thus calls for sophisticated many-body approaches, such as DMFT.
We derive a simple formula (Eq. (4)) relating the reduction of the magnetic anisotropy due to the Kondo screening to that of the ordered Ce moment; the latter being more easily accessible within DMFT. Moreover, this formula is of general applicability and can be used, for example, to extract the reduction of the anisotropy from the measured Ce magnetization. Our approach quantitatively captures the impact of Kondo screening on the intrinsic magnetic hardness κ of realistic Ce systems, thus opening the possibility for a theoretical search of Ce-based intermetallics with enhanced κ.
Our calculations reveal a partially Kondo-screened state of the Ce-4f magnetic moment in Ce-Fe intermetallics of the "1-12" family. Nitrogen interstitials strongly enhance the screening of the 4f moment, thus completely suppressing the Ce single-ion contribution to the anisotropy both at liquid helium and room temperature. These theoretical predictions are validated by our experimental measurements finding the anisotropy field in CeFe 11 Ti to be strongly suppressed by nitrogenation, in particular at low temperatures where the Ce contribution is the most significant. The effect of the dramatic Kondo-scale enhancement by nitrogen in the Ce-Fe "1-12" system indicates an unexpected sensitivity of this scale to interstitials in hard-magnetic intermetallics. The Kondo effect at low temperatures might thus be effectively controlled by such interstitials. Therefore, one may systematically search for interstitial elements that, in contrast to nitrogen in Ce-Fe "1-12", suppress the Kondo screening of the 4f moment. The present theoretical framework can be employed in this search.
At the same time, the small magnetic anisotropy of a purely local Ce-4f shell at elevated temperatures found in the present work limits the direct use of cerium-based "1-12" systems as hardmagnetic materials. Cerium can still be successfully used as a partial substitution element for the critical heavy rare-earth elements. In this case, the essential magnetic hardness is provided by heavy rare earths, but their concentration may in some cases be reduced without compromising the magnetic properties as recently shown in ref. 6 on the example of a mixture of La and Ce reducing the material cost of Nd 2 Fe 14 B-type hard-magnetic materials. However, the present study shows that cerium itself should not be expected to provide an essential contribution to the magnetic anisotropy at room temperature and above. Though in the present work we considered only the "1-12" family of Ce-Fe intermetalics, this conclusion based on the intrinsic properties of the Ce-4f 1 shell should be rather independent of the crystalline environment. It is expected thus to be of general validity for magnetic Ce intermetallics, unless one finds a-possibly exoticmechanism allowing the hybridization to enhance the Ce-4f anisotropy in the high-temperature local moment regime (i.e., above the Kondo temperature). Such an enhancement would manifest itself in a larger value of the Ce easy-axis magnetic Fig. 3 Impact of nitrogenation on electronic correlations and spectral function of CeFe 11 Ti. a Real part of the self-energy Σ(ω) in CeFe 11 Ti (in red) and CeFe 11 TiN (in blue). The quasiparticle weight Z and Ce-4f shell occupation n f are also indicated in the corresponding colors. b Imaginary part of the hybridization function Δ(ω). In particular, at the Fermi energy ImΔ(ω) is considerably larger in the nitrogenated compound. (The self-energies (a) and hybridization functions (b) are averaged over the six orbital contributions with the largest occupation number, respectively. For orbital-resolved self-energies, we refer the reader to the Supplementary Fig. 3.) c, d Corresponding spectral functions. The Ce-4f Kondo peak at the Fermi energy is more pronounced in CeFe 11 TiN due to the more delocalized character of its Ce-4f electrons. For completeness, also the Fe-3d, Ti-3d, and N-2p states are shown (the Fe-3d contribution has been rescaled by a factor 1/3). moment obtained within a full DMFT treatment, as compared to the value obtained in the quasi-atomic approximation. In the current "1-12" systems, such an effect is not observed. Given the environmental and economical impact that is at stake if improved permanent magnet materials are found, it might however be worthwhile to systematically search for it.

METHODS Theory
Our ab initio calculations-combined DFT and DMFT 27,30-32 -describe many-body effects in the Ce-4f shell taking into account on equal footing the key competing energy scales of the problem: the CF and exchange field induced by the spin-polarized electronic structure, as well as the spin-orbit coupling. We start from a charge-self-consistent spin-polarized DFT + DMFT calculation 33-37 of our target compounds, where we take into account the magnetization of Fe within the LSDA, with the Fe moment aligned along the c||z axis, while we treat the Ce-4f magnetic moment within the quasi-atomic Hubbard-I approximation 26,27 . In our selfconsistent DFT + DMFT calculations within the Hubbard-I approximation, abbreviated by DFT + HubI, we suppress the self-interaction contribution to the CF splitting, as well as the 4f contribution to the exchange field B ex on the Ce site. In result, the DFT exchange-correlation potential is calculated using a non-spin-polarized and spherical Ce-4f electronic density, while the nonspherical contributions and magnetization density of other states (including non-4f Ce orbitals like Ce-5d) are fully taken into account. This approach provides reliable CFs and B ex in hard-magnetic intermetallics, as previously demonstrated in ref. 16 . We employ the fully localized limit double-counting correction, in conjunction with the nominal n = 1 occupancy for the Ce-4f shell 38 . Since we obtain the Fe exchange field B ex within the DFT + HubI calculations, it is temperature independent. This is a reasonable approximation for the temperature range <300 K that we focus on. The experimental Curie temperatures T c for CeFe 11 Ti and CeFe 11 TiN of 485 and 710 K, respectively, are significantly >300 K. Our tests within the quasi-atomic treatment of the Ce-4f shell show that including the finite-temperature decay B ex leads to a rather small (<10% at 300 K) additional reduction of the Ce moment, as compared to the dependence (Fig. 2a) obtained with the fixed zero-T value of B ex .
Having converged the self-consistent DFT + HubI calculations, we use the resulting charge density to produce the Kohn-Sham one-electron part of the Hamiltonian that includes the correct CF and B ex acting on the Ce-4f shell. We then perform an additional DMFT run using a hybridizationexpansion continuous-time QMC solver in the segment picture 25,39 . We abbreviate this approach as DMFT(QMC). By using in DMFT(QMC) the charge and magnetization density obtained in DFT + HubI, we neglect the feedback of possible Ce-4f delocalization to the Fe magnetization and the exchange field B ex . This treatment is in-line with the usual assumption that the 4f magnetism in the TM-RE hard magnets has a negligible influence on the TM-sublattice magnetism.
In the Ce-4f quantum impurity problem, we employ the basis of f 1 eigenstates of the atomic Hamiltonian, which represents two manifolds, j = 5/2 and j = 7/2, split by the spin-orbit coupling, with additional smaller splittings within the manifolds due to the CF and EF. We treat all 14 Ce-4f orbitals as correlated and employ the density-density approximation for the corresponding Coulomb vertex. We calculate the ordered Ce-4f magnetic moment M z ¼ Tr ðL z þ 2Ŝ z Þρ Â Ã from the converged DMFT density matrix ρ.
The Ce single-ion anisotropy constant within the atomic approximation, K at 1 , is evaluated from the dependence of the anisotropy energy E at a on the azimutal angle θ of the exchange field B ex by fitting it to the form E at a ¼ K at 1 sin 2 θ. The dependence E at a vs. θ is calculated by a numerical diagonalization of the Ce-4f CF Hamiltonian under B ex applied in the corresponding direction. Please refer to the Supplementary Note 2 for further details regarding our ab initio calculations and the calculation of the magnetic anisotropy.

Experiment
Polycrystalline CeFe 11 Ti samples are produced by the suction casting method using commercial purity elements (>99.99%). The as-cast samples are wrapped in Mo foil, sealed in a quartz ampule under argon atmosphere and annealed at 850°C for 12 h followed by quenching in water. Part of the suction-cast plates are ground into powder <20 μm to carry out the nitrogenation. The determination of the crystal structures is done using a room-temperature x-ray powder diffractometer (XRD) and the XRD patterns are fitted with a ThMn 12 -type (I4/mmm, space group no. 139) structure. Thermomagnetic measurements are carried out for the determination of the Curie temperatures. The anisotropy fields are measured on the polycrystalline CeFe 11 Ti and CeFe 11 TiN samples by using the method suggested by Durst and Kronmüller 40 . Applications of this model to ThMn 12 systems have been reported recently for similar material systems 41 . Further details of the sample production and characterization can be found in the Supplementary Note 1.