Bachelorthesis: Calculation of the magnetic properties of quarternary ThMn$_{12}$-type compounds with Zr as a substitution for Nd

Cette thèse de licence utilise des calculs de la théorie de la fonctionnelle de la densité pour étudier les propriétés intrinsèques et magnétiques des composés de type ThMn$_{12}$ substitués au Nd avec Zr et Ti, identifiant des matériaux quaternaires prometteurs à faible teneur en Nd comme des alternatives potentielles aux aimants aux terres rares rares.

L'essentiel

La Grande Image : Trouver un Aimant Moins Cher et Plus Écologique

Imaginez que le monde moderne est une gigantesque machine fonctionnant à l'électricité. Pour faire fonctionner cette machine efficacement (comme dans les voitures électriques, les éoliennes et les disques durs), nous avons besoin d'aimants puissants. Actuellement, la « référence absolue » pour ces aimants est un matériau appelé Nd-Fe-B (Néodyme-Fer-Bore).

Considérez le Néodyme (Nd) comme le VIP d'une fête. Il rend l'aimant incroyablement puissant, mais il est cher, difficile à trouver, et presque tout le monde dépend d'un seul pays (la Chine) pour son approvisionnement. Cela crée un goulot d'étranglement dans la chaîne d'approvisionnement, comme un pont étroit unique que tout le monde tente de traverser.

L'objectif de cette recherche est de trouver un substitut à ce VIP. L'auteur, Nico Yannik Merkt, suggère d'utiliser le Zirconium (Zr). Le Zirconium est comme le voisin fiable et abordable : il est abondant, moins cher et plus facile à obtenir. La question est : si nous remplaçons le VIP par le voisin, la fête (l'aimant) fonctionne-t-elle toujours aussi bien ?

Le Problème : La « Maison Instable »

Le type spécifique de structure d'aimant étudié est appelé ThMn12 (ou la « phase 1:12 »).

• Le Plan : Imaginez une maison construite selon un plan spécifique où vous avez 1 atome de terre rare (le VIP) et 12 atomes de fer.
• Le Problème : Si vous essayez de construire cette maison uniquement avec le VIP (Néodyme) et les 12 atomes de fer, la maison est instable. C'est comme essayer de construire un gratte-ciel sur un fondation de sable ; elle s'effondre.
• La Solution : Pour que la maison tienne debout, vous avez besoin d'un « stabilisateur ». Dans ce cas, les chercheurs utilisent le Titane (Ti). Imaginez le Titane comme les poutres en acier que vous ajoutez à la structure pour empêcher la maison de s'effondrer.

L'Expérience : Un Chantier de Construction Virtuel

Puisque construire ces aimants dans un vrai laboratoire est coûteux et long, l'auteur a utilisé des superordinateurs pour simuler la construction. Cela s'appelle la Théorie de la Fonctionnelle de la Densité (DFT).

• La Simulation : Au lieu de mélanger des produits chimiques dans un bécher, l'ordinateur calcule comment les atomes « ressentent » les uns les autres. Il se demande : « Si je mets du Zirconium ici, la maison restera-t-elle debout ? Quelle sera la force de l'attraction magnétique ? »

Ce Que L'Ordinateur a Découvert

Le papier examine plusieurs scénarios « et si » pour voir comment le remplacement du Néodyme par le Zirconium affecte les performances de l'aimant. Voici les principales découvertes :

1. Stabilité (La maison va-t-elle tenir ?)

• Néodyme Pur : Sans aide, la maison est instable.
• Zirconium Pur : Étonnamment, une maison faite entièrement de Zirconium et de Fer est stable.
• Le Mélange (50/50) : Lorsqu'ils ont mélangé la moitié Néodyme et la moitié Zirconium, la maison était un peu vacillante. Ils avaient besoin d'ajouter plus de « poutres en acier » (Titane) pour la maintenir stable.
• Conclusion : Vous pouvez remplacer le Néodyme par du Zirconium, mais vous devez faire attention à la recette pour maintenir la structure stable.

2. Force (Quelle est la puissance de l'aimant ?)

• Le Compromis : Le VIP (Néodyme) est naturellement très magnétique. Le voisin (Zirconium) l'est moins.
• Le Résultat : Lorsqu'ils ont remplacé le Néodyme par du Zirconium, l'aimant est devenu légèrement plus faible. C'est comme remplacer une corde super résistante par une corde légèrement plus fine. Cependant, l'aimant reste très puissant – suffisamment puissant pour être utile.
• Le « Produit Énergétique » : C'est une mesure de la quantité d'énergie que l'aimant peut stocker. Les nouveaux aimants à base de Zirconium ont obtenu des scores très élevés, battant certains types d'aimants plus anciens et se rapprochant du champion actuel (Nd-Fe-B).

3. Résistance à la Chaleur (La Température de Curie)

• Les aimants perdent leur puissance s'ils deviennent trop chauds. La « Température de Curie » est le point où l'aimant abandonne et cesse de fonctionner.
• La Découverte : Les nouveaux aimants au Zirconium peuvent supporter la chaleur presque aussi bien que ceux au Néodyme. Ils ne fondront pas et ne perdront pas leur puissance dans un moteur électrique chaud.

4. Directionnalité (La « Rue à Sens Unique »)

• Un bon aimant permanent doit être « dur » à désaimanter. Il doit maintenir sa direction fermement, comme une rue à sens unique.
• La Découverte : Les aimants au Zirconium sont très bons pour maintenir leur direction. En fait, dans certains calculs, les aimants au Zirconium étaient même meilleurs pour maintenir leur direction que ceux au Néodyme.

Le Verdict : Est-ce un Gagnant ?

Le papier conclut que le Zirconium est un substitut très prometteur au Néodyme.

• Les Avantages : Il est moins cher, plus abondant et moins critique pour les chaînes d'approvisionnement. Les aimants résultants sont stables et possèdent d'excellentes propriétés magnétiques.
• Les Inconvénients : Les aimants sont légèrement plus faibles que ceux en Néodyme pur, et ils sont actuellement juste en dessous d'être des aimants « parfaitement durs » (ils sont « semi-durs »).
• L'Avenir : L'auteur suggère qu'avec un peu plus de réglages (comme l'ajout d'Azote ou l'ajustement de la recette), ces aimants au Zirconium pourraient devenir une alternative réelle aux aimants au Néodyme coûteux que nous utilisons aujourd'hui.

En résumé : L'auteur a utilisé un ordinateur pour prouver que nous pouvons construire un aimant puissant et stable en utilisant le voisin bon marché et abondant (Zirconium) au lieu du VIP coûteux (Néodyme). Ce n'est pas tout à fait aussi puissant que la version VIP pour l'instant, mais c'est assez proche pour pouvoir révolutionner la façon dont nous fabriquons des aimants pour les voitures électriques et les énergies vertes.

Résumé technique

1. Énoncé du problème et motivation

La demande mondiale d'aimants permanents haute performance est tirée par les technologies d'énergie renouvelable et les véhicules électriques. Actuellement, l'industrie dépend fortement des aimants Nd$_2$Fe$_{14}$B. Cependant, le Néodyme (Nd) et d'autres éléments de terres rares (TR) sont considérés comme des ressources critiques en raison des monopoles de la chaîne d'approvisionnement (principalement la Chine) et de la diminution des réserves.

• Le défi : Développer des aimants « pauvres en TR » ou sans TR qui maintiennent de hautes performances magnétiques sans dépendre d'éléments critiques.
• La lacune spécifique : Bien que la structure de type ThMn$_{12}$ (1:12) (RFe$_{12}$) offre une teneur en TR inférieure à celle du Nd$_2$Fe$_{14}$B, les phases pures RFe$_{12}$ sont thermodynamiquement instables. Elles nécessitent une stabilisation (généralement via le Titane, Ti) et souffrent souvent de températures de Curie ou d'anisotropie plus faibles par rapport au Nd$_2$Fe$_{14}$B.
• L'objectif : Ce travail étudie le Zirconium (Zr) comme substitut au Nd dans la phase 1:12. Le Zr est abondant, peu coûteux et non critique. L'étude vise à déterminer si les composés quaternaires (Zr,Nd)Fe$_{12-y}$Ti$_y$ et quinaires (Zr,Nd)Fe$_{10}$CoTi substitués au Zr peuvent servir d'alternatives viables et haute performance aux aimants Nd-Fe-B.

2. Méthodologie

La recherche utilise des méthodes de calcul ab initio basées sur la Théorie de la Fonctionnelle de la Densité (DFT) pour prédire les propriétés magnétiques intrinsèques.

• Logiciels et cadre :
  • VASP (Vienna Ab-initio Simulation Package) : Utilisé pour la relaxation structurale et le calcul des énergies totales, des moments magnétiques et de l'anisotropie magnétocristalline.
  • AkaiKKR (MACHIKANEYAMA) : Utilisé pour calculer les températures de Curie ($T_C$) via la méthode de la fonction de Green Korringa-Kohn-Rostoker (KKR).
• Approches théoriques :
  • Fonctionnelles d'échange-corrélation : Principalement l'Approximation du Gradient Généralisé (GGA-PBE). Pour l'énergie d'anisotropie magnétocristalline (MAE), l'Approximation de la Densité de Spin Locale (LSDA) a également été employée, car elle donne souvent des constantes d'anisotropie plus précises pour ces systèmes.
  • Traitement des électrons f : Les électrons 4f du Nd ont été traités via la méthode DFT+U (approche Dudarev, $U = 6$ eV) pour corriger les erreurs d'auto-interaction et la forte corrélation électronique, que la DFT standard ne parvient pas à capturer avec précision.
  • Pseudopotentiels : Méthode des ondes augmentées par projecteur (PAW).
• Calculs effectués :
  • Stabilité de phase : Calcul des énergies de formation et des enthalpies de solution pour identifier les configurations stables des substitutions de Ti et Co sur les sites Fe et de Zr sur les sites TR.
  • Propriétés magnétiques : Moment magnétique total ($m_{tot}$), aimantation à saturation ($\mu_0 M_S$) et produit énergétique maximal ($|BH|_{max}$).
  • Anisotropie : Calcul de l'énergie d'anisotropie magnétocristalline (MAE) en faisant tourner les directions d'aimantation ([001] vers [010]) pour déduire la constante d'anisotropie ($K_1$) et le champ d'anisotropie ($H_a$).
  • Température de Curie ($T_C$) : Calculée en utilisant les approximations du champ moyen (MFA) et du désordre du moment local (LMD) pour tenir compte des fluctuations thermiques.
  • Facteur de dureté ($\kappa$) : Calculé pour déterminer si le matériau qualifie comme un aimant dur ($\kappa > 1$).

3. Contributions clés

1. Analyse systématique de la substitution par Zr : Ce travail fournit la première investigation théorique complète du Zr remplaçant le Nd dans la structure ThMn$_{12}$ 1:12, allant au-delà des études précédentes axées sur Y ou Ce.
2. Mécanismes de stabilisation : Il confirme que bien que le (Zr,Nd)Fe$_{12}$ pur soit instable, l'ajout de Ti (sur le site 8i) et de Co (sur les sites 8f/8j) stabilise la phase. Il établit que la substitution par Zr permet de réduire la concentration de Ti par rapport aux analogues purement à base de Nd.
3. Validation méthodologique : L'étude valide l'utilisation de la LSDA pour le calcul de l'anisotropie dans ces systèmes, notant que GGA+U sous-estime souvent $K_1$ pour les composés contenant du Nd, tandis que la LSDA fournit des résultats plus proches des données expérimentales.
4. Identification de candidats prometteurs : L'identification de (Zr,Nd)Fe$_{11}$Ti et (Zr,Nd)Fe$_{10}$CoTi comme candidats prometteurs équilibrant coût (réduction du Nd) et performance.

4. Résultats clés

A. Stabilité et structure

• Stabilisation par Ti : Le Ti occupe préférentiellement le site cristallographique 8i. Une concentration d'environ 7,7 % at. de Ti (2 atomes par supercellule de 26 atomes) suffit à stabiliser la phase (Zr,Nd)Fe$_{12}$.
• Substitution par Co : Le Co préfère les sites 8f et 8j. Bien que le Co améliore légèrement les moments magnétiques, son effet stabilisateur est moindre que celui du Ti.
• Paramètres de maille : La substitution par Zr entraîne une réduction du volume de la maille par rapport aux composés Nd purs en raison du rayon atomique plus petit du Zr, suivant la tendance : $Nd_2Fe_{24} > ZrNdFe_{24} > Zr_2Fe_{24}$.

B. Moments magnétiques et aimantation à saturation ($\mu_0 M_S$)

• Effet du Zr : Le remplacement du Nd par le Zr réduit le moment magnétique total car le Zr possède un moment magnétique plus faible (-0,47 $\mu_B$) comparé au Nd (-0,27 $\mu_B$ en GGA, mais -3,3 $\mu_B$ lorsque les électrons f sont inclus via DFT+U).
• Valeurs calculées :
  • NdFe$_{11}$Ti : $\mu_0 M_S \approx 1,69$ T (DFT+U), correspondant aux plages expérimentales (1,38–1,70 T).
  • (Zr,Nd)Fe$_{11}$Ti : $\mu_0 M_S \approx 1,62$ T.
  • ZrFe$_{11}$Ti : $\mu_0 M_S \approx 1,51$ T.
• Conclusion : Bien que la substitution par Zr abaisse légèrement $M_S$, les valeurs restent suffisamment élevées pour des applications pratiques.

C. Produit énergétique maximal ($|BH|_{max}$)

• Performance : Le composé quaternaire (Zr,Nd)Fe$_{11.5}$Ti$_{0.5}$ montre un $|BH|_{max}$ théorique d'environ 600 kJ/m³, dépassant la limite théorique du Nd$_2$Fe$_{14}$B (~512 kJ/m³).
• Phase stable : Le composé stable (Zr,Nd)Fe$_{11}$Ti produit environ 510–525 kJ/m³, comparable au Nd$_2$Fe$_{14}$B et significativement supérieur au Sm$_2$Co$_{17}$ (294 kJ/m³).
• Effet du Co : La substitution par Co augmente légèrement $M_S$ mais diminue marginalement $|BH|_{max}$ en raison des changements de volume.

D. Température de Curie ($T_C$)

• Tendance : $T_C$ diminue avec l'augmentation de la teneur en Ti mais augmente avec la substitution par Co.
• Valeurs :
  • NdFe$_{11}$Ti : $T_C \approx 652$ K (LMD) / 792 K (FM).
  • (Zr,Nd)Fe$_{11}$Ti : $T_C \approx 652$ K (LMD) / 783 K (FM).
  • Comparaison : Ces valeurs sont plus élevées que celles du Nd$_2$Fe$_{14}$B ($T_C \approx 588$ K), les rendant adaptées à des applications à plus haute température.
• Note : L'approximation FM surestime $T_C$ ; les valeurs LMD sont considérées plus réalistes mais indiquent toujours une température de fonctionnement viable.

E. Anisotropie magnétocristalline ($K_1$) et facteur de dureté ($\kappa$)

• Constante d'anisotropie ($K_1$) : La substitution par Zr augmente significativement $K_1$.
  • NdFe$_{11}$Ti : $K_1 \approx 1,68$ MJ/m³ (LSDA).
  • (Zr,Nd)Fe$_{11}$Ti : $K_1 \approx 2,41$ MJ/m³ (LSDA).
  • ZrFe$_{11}$Ti : $K_1 \approx 2,47$ MJ/m³ (LSDA).
  • Signification : Le Zr améliore l'anisotropie, une propriété cruciale pour la coercivité.
• Facteur de dureté ($\kappa$) :
  • (Zr,Nd)Fe$_{11}$Ti donne $\kappa \approx 0,97$ (LSDA).
  • Cela classe le matériau comme semi-dur ($\kappa < 1$), tandis que le Nd$_2$Fe$_{14}$B est dur ($\kappa \approx 1,54$).
  • La valeur est très proche du seuil des aimants durs, suggérant que de légères modifications (par exemple, la nitruration) pourraient le faire basculer dans le régime des aimants durs.

5. Importance et perspectives

• Importance stratégique : Ce travail démontre que le Zirconium est un substitut viable, peu coûteux et abondant au Néodyme dans les aimants de type ThMn$_{12}$ sans sacrifier les propriétés magnétiques critiques telles que l'aimantation à saturation ou la température de Curie.
• Potentiel de performance : Les propriétés calculées suggèrent que les aimants à base de (Zr,Nd)Fe$_{11}$Ti pourraient rivaliser avec les Nd-Fe-B en termes de produit énergétique et de stabilité thermique, tout en offrant une chaîne d'approvisionnement plus sécurisée.
• Axes futurs :
  • Vérification expérimentale : Les prédictions théoriques nécessitent une synthèse (par exemple, via fusion sélective par laser) et une validation expérimentale.
  • Optimisation : Puisque le facteur de dureté est légèrement inférieur à 1,0, les travaux futurs devraient explorer la nitruration (formation de phases R-Fe-N) ou un ajustement supplémentaire de Co/Ti pour pousser $\kappa > 1,0$.
  • Extension : La méthodologie peut être appliquée à d'autres structures pauvres en TR (phases 2:14, 2:17) ou à d'autres éléments abondants (par exemple, La, Ce) pour réduire davantage la dépendance aux terres rares critiques.

En conclusion, la thèse identifie avec succès les composés (Zr,Nd)Fe$_{12-y}$Ti$_y$ comme des candidats prometteurs pour la prochaine génération d'aimants permanents durables et haute performance.