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

Diese Bachelorarbeit nutzt Dichtefunktionaltheorie-Rechnungen, um die intrinsischen und magnetischen Eigenschaften von Nd-substituierten ThMn$_{12}$-Typ-Verbindungen mit Zr und Ti zu untersuchen und vielversprechende Nd-arme quaternäre Materialien als potenzielle Alternativen zu seltenen Erdmagneten zu identifizieren.

Das Wesentliche

Das große Ganze: Einen günstigeren, umweltfreundlicheren Magneten finden

Stellen Sie sich die moderne Welt als eine riesige Maschine vor, die mit Elektrizität läuft. Damit diese Maschine effizient funktioniert (wie bei Elektroautos, Windkraftanlagen und Festplatten), benötigen wir leistungsstarke Magnete. Derzeit ist das „Goldstandard"-Material für diese Magnete eine Verbindung namens Nd-Fe-B (Neodym-Eisen-Bor).

Stellen Sie sich Neodym (Nd) als den VIP-Gast auf einer Party vor. Er macht den Magneten unglaublich stark, ist aber teuer, schwer zu finden, und fast jeder ist darauf angewiesen, dass nur ein einziges Land (China) ihn liefert. Dies erzeugt einen Engpass in der Lieferkette, ähnlich wie eine einzige schmale Brücke, die alle überqueren wollen.

Das Ziel dieser Forschung ist es, einen Ersatz für diesen VIP-Gast zu finden. Der Autor, Nico Yannik Merkt, schlägt die Verwendung von Zirkonium (Zr) vor. Zirkonium ist wie der zuverlässige, erschwingliche Nachbar, der reichlich vorhanden, günstiger und leichter zu beschaffen ist. Die Frage lautet: Wenn wir den VIP gegen den Nachbarn austauschen, funktioniert die Party (der Magnet) dann immer noch genauso gut?

Das Problem: Das „instabile Haus"

Die spezifische Art der untersuchten Magnetstruktur wird ThMn12 (oder die „1:12-Phase") genannt.

• Der Bauplan: Stellen Sie sich ein Haus vor, das nach einem bestimmten Bauplan errichtet wird, bei dem Sie 1 Seltenerd-Atom (den VIP) und 12 Eisen-Atome haben.
• Das Problem: Wenn Sie versuchen, dieses Haus nur mit dem VIP (Neodym) und den 12 Eisen-Atomen zu bauen, ist das Haus instabil. Es ist, als würde man versuchen, einen Wolkenkratzer auf einem Sandfundament zu errichten; es stürzt ein.
• Die Lösung: Damit das Haus steht, benötigen Sie einen „Stabilisator". In diesem Fall verwenden die Forscher Titan (Ti). Denken Sie an Titan als die Stahlträger, die Sie in den Rahmen einfügen, damit das Haus nicht auseinanderfällt.

Das Experiment: Eine virtuelle Baustelle

Da der Bau dieser Magnete in einem echten Labor teuer und zeitaufwendig ist, verwendete der Autor Supercomputer, um den Bau zu simulieren. Dies wird als Dichtefunktionaltheorie (DFT) bezeichnet.

• Die Simulation: Anstatt Chemikalien in einem Becherglas zu mischen, berechnet der Computer, wie die Atome „einander spüren". Er fragt: „Wenn ich hier Zirkonium einsetze, bleibt das Haus stehen? Wie stark wird der magnetische Zug sein?"

Was der Computer herausfand

Das Papier durchläuft mehrere „Was-wäre-wenn"-Szenarien, um zu sehen, wie sich der Austausch von Neodym gegen Zirkonium auf die Leistung des Magneten auswirkt. Hier sind die wichtigsten Erkenntnisse:

1. Stabilität (Wird das Haus stehen?)

• Reines Neodym: Ohne Hilfe ist das Haus instabil.
• Reines Zirkonium: Überraschenderweise ist ein Haus, das vollständig aus Zirkonium und Eisen besteht, stabil.
• Die Mischung (50/50): Als sie die Hälfte Neodym und die Hälfte Zirkonium mischten, war das Haus etwas wackelig. Sie mussten mehr „Stahlträger" (Titan) hinzufügen, um es stabil zu halten.
• Fazit: Sie können Neodym durch Zirkonium ersetzen, müssen aber das Rezept sorgfältig wählen, um die Struktur stabil zu halten.

2. Stärke (Wie stark ist der Magnet?)

• Der Kompromiss: Der VIP (Neodym) ist von Natur aus sehr magnetisch. Der Nachbar (Zirkonium) ist weniger magnetisch.
• Das Ergebnis: Als sie Neodym gegen Zirkonium austauschten, wurde der Magnet etwas schwächer. Es ist, als würde man ein superstarkes Seil durch ein etwas dünneres ersetzen. Dennoch ist der Magnet immer noch sehr stark – stark genug, um nützlich zu sein.
• Der „Energieprodukt": Dies ist ein Maß dafür, wie viel Energie der Magnet speichern kann. Die neuen Zirkonium-basierten Magneten erzielten sehr hohe Werte, schlugen einige ältere Magnettypen und kamen dem aktuellen Champion (Nd-Fe-B) nahe.

3. Hitzebeständigkeit (Die Curie-Temperatur)

• Magnete verlieren ihre Kraft, wenn sie zu heiß werden. Die „Curie-Temperatur" ist der Punkt, an dem der Magnet aufgibt und aufhört zu funktionieren.
• Die Erkenntnis: Die neuen Zirkonium-Magnete können Hitze fast genauso gut verkraften wie die Neodym-Magnete. Sie werden in einem heißen Elektromotor nicht schmelzen oder ihre Kraft verlieren.

4. Richtungsabhängigkeit (Die „Einbahnstraße")

• Ein guter Permanentmagnet muss „hart" sein, um entmagnetisiert zu werden. Er muss seine Richtung festhalten, wie eine Einbahnstraße.
• Die Erkenntnis: Die Zirkonium-Magnete sind sehr gut darin, ihre Richtung zu halten. Tatsächlich waren die Zirkonium-Magnete in einigen Berechnungen sogar besser darin, ihre Richtung zu halten, als die Neodym-Magnete.

Das Urteil: Ist es ein Gewinner?

Das Papier kommt zu dem Schluss, dass Zirkonium ein sehr vielversprechender Ersatz für Neodym ist.

• Die Vorteile: Es ist günstiger, reichlicher vorhanden und weniger kritisch für Lieferketten. Die resultierenden Magnete sind stabil und weisen hervorragende magnetische Eigenschaften auf.
• Die Nachteile: Die Magnete sind etwas schwächer als die reinen Neodym-Magnete, und sie liegen derzeit knapp darunter, „perfekt harte" Magnete zu sein (sie sind „halbhart").
• Die Zukunft: Der Autor schlägt vor, dass diese Zirkonium-Magnete mit ein wenig mehr Feintuning (wie dem Hinzufügen von Stickstoff oder der Anpassung des Rezepts) zu einer realen Alternative zu den teuren Neodym-Magneten werden könnten, die wir heute verwenden.

Kurz gesagt: Der Autor nutzte einen Computer, um nachzuweisen, dass wir einen starken, stabilen Magneten mit dem günstigen, reichlich vorhandenen Nachbarn (Zirkonium) statt des teuren VIPs (Neodym) bauen können. Es ist noch nicht ganz so stark wie die VIP-Version, aber nah genug, um die Art und Weise, wie wir Magnete für Elektroautos und grüne Energie herstellen, revolutionieren zu können.

Technische Zusammenfassung

1. Problemstellung und Motivation

Die globale Nachfrage nach Hochleistungs-Permanentmagneten wird durch erneuerbare Energietechnologien und Elektrofahrzeuge getrieben. Derzeit verlässt sich die Industrie stark auf Nd$_2$Fe$_{14}$B-Magnete. Neodym (Nd) und andere Seltene-Erden-(SE)-Elemente gelten jedoch aufgrund von Monopolen in der Lieferkette (hauptsächlich China) und schwindenden Reserven als kritische Ressourcen.

• Die Herausforderung: Entwicklung von „SE-armen" oder SE-freien Magneten, die eine hohe magnetische Leistung beibehalten, ohne auf kritische Elemente angewiesen zu sein.
• Die spezifische Lücke: Während die ThMn$_{12}$-Typ-(1:12)-Struktur (RFe$_{12}$) einen geringeren SE-Gehalt als Nd$_2$Fe$_{14}$B aufweist, sind reine RFe$_{12}$-Phasen thermodynamisch instabil. Sie erfordern eine Stabilisierung (üblicherweise durch Titan, Ti) und leiden oft unter niedrigeren Curie-Temperaturen oder Anisotropie im Vergleich zu Nd$_2$Fe$_{14}$B.
• Das Ziel: Diese Arbeit untersucht Zirkonium (Zr) als Ersatz für Nd in der 1:12-Phase. Zr ist reichlich vorhanden, kostengünstig und nicht kritisch. Die Studie zielt darauf ab, zu bestimmen, ob Zr-substituierte quaternäre (Zr,Nd)Fe$_{12-y}$Ti$_y$- und quinäre (Zr,Nd)Fe$_{10}$CoTi-Verbindungen als tragfähige, hochleistungsfähige Alternativen zu Nd-Fe-B-Magneten dienen können.

2. Methodik

Die Forschung verwendet ab-initio-Rechenmethoden auf Basis der Dichtefunktionaltheorie (DFT), um intrinsische magnetische Eigenschaften vorherzusagen.

• Software und Rahmenwerk:
  • VASP (Vienna Ab-initio Simulation Package): Verwendet für strukturelle Relaxation und Berechnung von Gesamtenergien, magnetischen Momenten und magnetokristalliner Anisotropie.
  • AkaiKKR (MACHIKANEYAMA): Verwendet zur Berechnung der Curie-Temperaturen ($T_C$) mittels der Korringa-Kohn-Rostoker (KKR)-Green-Funktions-Methode.
• Theoretische Ansätze:
  • Austausch-Korrelations-Funktionale: Hauptsächlich die Verallgemeinerte Gradientennäherung (GGA-PBE). Für die magnetokristalline Anisotropieenergie (MAE) wurde zusätzlich die Lokale Spindichtennäherung (LSDA) eingesetzt, da sie für diese Systeme oft genauere Anisotropiekonstanten liefert.
  • Behandlung von f-Elektronen: Nd 4f-Elektronen wurden mit der DFT+U-Methode (Dudarev-Ansatz, $U = 6$ eV) behandelt, um Selbstwechselwirkungsfehler und starke Elektronenkorrelationen zu korrigieren, die von der Standard-DFT nicht korrekt erfasst werden.
  • Pseudopotentiale: Projektor-erweiterte Wellen-Methode (PAW).
• Durchgeführte Berechnungen:
  • Phasenstabilität: Berechnung von Bildungsenergien und Lösungsenthalpien zur Identifizierung stabiler Konfigurationen von Ti- und Co-Substitutionen an Fe-Stellen sowie Zr an SE-Stellen.
  • Magnetische Eigenschaften: Gesamtmagnetisches Moment ($m_{tot}$), Sättigungsmagnetisierung ($\mu_0 M_S$) und maximale Energieprodukt ($|BH|_{max}$).
  • Anisotropie: Berechnung der magnetokristallinen Anisotropieenergie (MAE) durch Rotation der Magnetisierungsrichtungen ([001] bis [010]), um die Anisotropiekonstante ($K_1$) und das Anisotropiefeld ($H_a$) abzuleiten.
  • Curie-Temperatur ($T_C$): Berechnet unter Verwendung der Molekularfeldnäherung (MFA) und des Modells für lokale Momenten-Unordnung (LMD), um thermische Fluktuationen zu berücksichtigen.
  • Härtefaktor ($\kappa$): Berechnet, um zu bestimmen, ob das Material als Hartmagnet qualifiziert ($\kappa > 1$).

3. Hauptbeiträge

1. Systematische Zr-Substitutionsanalyse: Diese Arbeit bietet die erste umfassende theoretische Untersuchung von Zr als Ersatz für Nd in der 1:12 ThMn$_{12}$-Struktur und geht über frühere Studien hinaus, die sich auf Y oder Ce konzentrierten.
2. Stabilisierungsmechanismen: Es wird bestätigt, dass zwar reines (Zr,Nd)Fe$_{12}$ instabil ist, die Zugabe von Ti (auf der 8i-Stelle) und Co (auf den 8f/8j-Stellen) die Phase jedoch stabilisiert. Es wird festgestellt, dass die Zr-Substitution eine Reduzierung der Ti-Konzentration im Vergleich zu reinen Nd-basierten Analoga ermöglicht.
3. Methodische Validierung: Die Studie validiert die Verwendung von LSDA zur Berechnung der Anisotropie in diesen Systemen und stellt fest, dass GGA+U $K_1$ für Nd-haltige Verbindungen oft unterschätzt, während LSDA Ergebnisse liefert, die näher an experimentellen Daten liegen.
4. Identifizierung vielversprechender Kandidaten: Die Identifizierung von (Zr,Nd)Fe$_{11}$Ti und (Zr,Nd)Fe$_{10}$CoTi als vielversprechende Kandidaten, die Kosten (reduziertes Nd) und Leistung in Einklang bringen.

4. Hauptergebnisse

A. Stabilität und Struktur

• Ti-Stabilisierung: Ti besetzt bevorzugt die 8i-kristallographische Stelle. Eine Konzentration von ~7,7 at.% Ti (2 Atome pro 26-Atome-Supercelle) reicht aus, um die (Zr,Nd)Fe$_{12}$-Phase zu stabilisieren.
• Co-Substitution: Co bevorzugt die 8f- und 8j-Stellen. Während Co die magnetischen Momente leicht verbessert, hat es einen geringeren stabilisierenden Effekt als Ti.
• Gitterparameter: Die Zr-Substitution führt im Vergleich zu reinen Nd-Verbindungen zu einer Verringerung des Zellvolumens aufgrund des kleineren Atomradius von Zr, gemäß dem Trend: $Nd_2Fe_{24} > ZrNdFe_{24} > Zr_2Fe_{24}$.

B. Magnetische Momente und Sättigungsmagnetisierung ($\mu_0 M_S$)

• Einfluss von Zr: Der Ersatz von Nd durch Zr verringert das Gesamtmagnetische Moment, da Zr ein niedrigeres magnetisches Moment (-0,47 $\mu_B$) aufweist als Nd (-0,27 $\mu_B$ in GGA, aber -3,3 $\mu_B$, wenn f-Elektronen über DFT+U einbezogen werden).
• Berechnete Werte:
  • NdFe$_{11}$Ti: $\mu_0 M_S \approx 1,69$ T (DFT+U), übereinstimmend mit experimentellen Bereichen (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.
• Fazit: Obwohl die Zr-Substitution $M_S$ leicht senkt, bleiben die Werte hoch genug für praktische Anwendungen.

C. Maximales Energieprodukt ($|BH|_{max}$)

• Leistung: Die quaternäre Verbindung (Zr,Nd)Fe$_{11.5}$Ti$_{0.5}$ zeigt ein theoretisches $|BH|_{max}$ von ~600 kJ/m³, was die theoretische Grenze von Nd$_2$Fe$_{14}$B (~512 kJ/m³) übersteigt.
• Stabile Phase: Die stabile (Zr,Nd)Fe$_{11}$Ti-Verbindung liefert ~510–525 kJ/m³, vergleichbar mit Nd$_2$Fe$_{14}$B und deutlich höher als Sm$_2$Co$_{17}$ (294 kJ/m³).
• Co-Effekt: Die Co-Substitution erhöht $M_S$ leicht, verringert aber $|BH|_{max}$ aufgrund von Volumenänderungen nur geringfügig.

D. Curie-Temperatur ($T_C$)

• Trend: $T_C$ nimmt mit steigendem Ti-Gehalt ab, steigt jedoch mit Co-Substitution.
• Werte:
  • 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).
  • Vergleich: Diese Werte sind höher als bei Nd$_2$Fe$_{14}$B ($T_C \approx 588$ K), was sie für Anwendungen bei höheren Temperaturen geeignet macht.
• Hinweis: Die FM-Näherung überschätzt $T_C$; die LMD-Werte gelten als realistischer, deuten aber dennoch auf eine tragfähige Betriebstemperatur hin.

E. Magnetokristalline Anisotropie ($K_1$) und Härtefaktor ($\kappa$)

• Anisotropiekonstante ($K_1$): Die Zr-Substitution erhöht $K_1$ signifikant.
  • 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).
  • Bedeutung: Zr verstärkt die Anisotropie, eine entscheidende Eigenschaft für die Koerzitivfeldstärke.
• Härtefaktor ($\kappa$):
  • (Zr,Nd)Fe$_{11}$Ti ergibt $\kappa \approx 0,97$ (LSDA).
  • Dies klassifiziert das Material als halbhart ($\kappa < 1$), wohingegen Nd$_2$Fe$_{14}$B hart ist ($\kappa \approx 1,54$).
  • Der Wert liegt sehr nahe an der Schwelle für Hartmagnete, was darauf hindeutet, dass geringfügige Modifikationen (z. B. Nitrierung) es in den Hartmagnetbereich bringen könnten.

5. Bedeutung und Ausblick

• Strategische Bedeutung: Diese Arbeit zeigt, dass Zirkonium ein tragfähiger, kostengünstiger und reichlich vorhandener Ersatz für Neodym in ThMn$_{12}$-Typ-Magneten ist, ohne kritische magnetische Eigenschaften wie Sättigungsmagnetisierung oder Curie-Temperatur zu beeinträchtigen.
• Leistungspotenzial: Die berechneten Eigenschaften deuten darauf hin, dass (Zr,Nd)Fe$_{11}$Ti-basierte Magnete in Bezug auf Energieprodukt und thermische Stabilität mit Nd-Fe-B konkurrieren könnten, während sie eine sicherere Lieferkette bieten.
• Zukünftige Richtungen:
  • Experimentelle Verifizierung: Die theoretischen Vorhersagen erfordern eine Synthese (z. B. mittels Selective Laser Melting) und experimentelle Validierung.
  • Optimierung: Da der Härtefaktor knapp unter 1,0 liegt, sollte zukünftige Arbeit die Nitrierung (Bildung von R-Fe-N-Phasen) oder eine weitere Feinabstimmung von Co/Ti untersuchen, um $\kappa > 1,0$ zu erreichen.
  • Erweiterung: Die Methodik kann auf andere SE-arme Strukturen (2:14, 2:17 Phasen) oder andere reichlich vorhandene Elemente (z. B. La, Ce) angewendet werden, um die Abhängigkeit von kritischen Seltenen Erden weiter zu verringern.

Zusammenfassend identifiziert die Bachelorarbeit erfolgreich (Zr,Nd)Fe$_{12-y}$Ti$_y$-Verbindungen als vielversprechende Kandidaten für die nächste Generation nachhaltiger, hochleistungsfähiger Permanentmagnete.