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

Questa tesi di laurea impiega calcoli di teoria del funzionale densità per indagare le proprietà intrinseche e magnetiche dei composti di tipo ThMn$_{12}$ sostituiti con Nd e contenenti Zr e Ti, identificando promettenti materiali quaternari a basso contenuto di Nd come potenziali alternative ai magneti di terre rare scarsi.

L'essenziale

Il Quadro Generale: Trovare un Magnete Più Economico ed Ecologico

Immaginate il mondo moderno come una gigantesca macchina che funziona a elettricità. Per far funzionare questa macchina in modo efficiente (come nelle auto elettriche, nelle turbine eoliche e negli hard disk), abbiamo bisogno di potenti magneti. Attualmente, lo "standard aureo" per questi magneti è un materiale chiamato Nd-Fe-B (Neodimio-Ferro-Boro).

Pensate al Neodimio (Nd) come all'ospite VIP di una festa. Rende il magnete incredibilmente potente, ma è costoso, difficile da reperire e quasi tutti dipendono da un solo paese (la Cina) per la sua fornitura. Questo crea un collo di bottiglia nella catena di approvvigionamento, come un unico ponte stretto che tutti cercano di attraversare.

L'obiettivo di questa ricerca è trovare un sostituto per quell'ospite VIP. L'autore, Nico Yannik Merkt, suggerisce di utilizzare lo Zirconio (Zr). Lo Zirconio è come il vicino affidabile ed economico: è abbondante, più economico e più facile da ottenere. La domanda è: se sostituiamo il VIP con il vicino, la festa (il magnete) funziona ancora altrettanto bene?

Il Problema: La "Casa Instabile"

Il tipo specifico di struttura magnetica in studio è chiamato ThMn12 (o "fase 1:12").

• Il Progetto: Immaginate una casa costruita con un progetto specifico in cui avete 1 atomo di terra rara (il VIP) e 12 atomi di ferro.
• Il Problema: Se provate a costruire questa casa solo con il VIP (Neodimio) e i 12 atomi di ferro, la casa è instabile. È come cercare di costruire un grattacielo su una fondazione di sabbia; crolla.
• La Soluzione: Per far stare in piedi la casa, serve un "stabilizzatore". In questo caso, i ricercatori utilizzano il Titanio (Ti). Pensate al Titanio come alle travi d'acciaio che aggiungete alla struttura per impedire che la casa crolli.

L'Esperimento: Un Cantiere Virtuale

Poiché costruire questi magneti in un vero laboratorio è costoso e richiede tempo, l'autore ha utilizzato supercomputer per simulare la costruzione. Questo è chiamato Teoria del Funzionale della Densità (DFT).

• La Simulazione: Invece di mescolare sostanze chimiche in un becher, il computer calcola come gli atomi "sentono" l'uno l'altro. Si chiede: "Se metto lo Zirconio qui, la casa rimarrà in piedi? Quanto sarà forte l'attrazione magnetica?"

Cosa Ha Trovato il Computer

Il documento esamina diversi scenari "cosa succederebbe se" per vedere come lo scambio tra Neodimio e Zirconio influisca sulle prestazioni del magnete. Ecco le scoperte chiave:

1. Stabilità (La casa rimarrà in piedi?)

• Neodimio Puro: Senza aiuto, la casa è instabile.
• Zirconio Puro: Sorprendentemente, una casa fatta interamente di Zirconio e Ferro è stabile.
• La Miscela (50/50): Quando hanno mescolato metà Neodimio e metà Zirconio, la casa era un po' traballante. Hanno dovuto aggiungere più "travi d'acciaio" (Titanio) per mantenerla stabile.
• Conclusione: Potete sostituire il Neodimio con lo Zirconio, ma dovete fare attenzione alla ricetta per mantenere stabile la struttura.

2. Forza (Quanto è forte il magnete?)

• Il Trade-off: Il VIP (Neodimio) è naturalmente molto magnetico. Il vicino (Zirconio) lo è meno.
• Il Risultato: Quando hanno sostituito il Neodimio con lo Zirconio, il magnete è diventato leggermente più debole. È come sostituire una corda super-resistente con una leggermente più sottile. Tuttavia, il magnete è ancora molto forte — abbastanza forte da essere utile.
• Il "Prodotto Energetico": Questa è una misura di quanta energia il magnete può immagazzinare. I nuovi magneti a base di Zirconio hanno ottenuto punteggi molto alti, battendo alcuni tipi di magneti più vecchi e avvicinandosi al campione attuale (Nd-Fe-B).

3. Resistenza al Calore (La Temperatura di Curie)

• I magneti perdono la loro potenza se si surriscaldano troppo. La "Temperatura di Curie" è il punto in cui il magnete si arrende e smette di funzionare.
• La Scoperta: I nuovi magneti di Zirconio possono gestire il calore quasi tanto bene quanto quelli di Neodimio. Non si scioglieranno né perderanno la loro potenza in un motore elettrico caldo.

4. Direzionalità (La "Strada a Senso Unico")

• Un buon magnete permanente deve essere "duro" da smagnetizzare. Deve mantenere la sua direzione fermamente, come una strada a senso unico.
• La Scoperta: I magneti di Zirconio sono molto bravi a mantenere la loro direzione. In effetti, in alcuni calcoli, i magneti di Zirconio erano addirittura migliori nel mantenere la loro direzione rispetto a quelli di Neodimio.

Il Verdetto: È un Vincitore?

Il documento conclude che lo Zirconio è un sostituto molto promettente per il Neodimio.

• I Pro: È più economico, più abbondante e meno critico per le catene di approvvigionamento. I magneti risultanti sono stabili e possiedono eccellenti proprietà magnetiche.
• I Contro: I magneti sono leggermente più deboli rispetto a quelli di Neodimio puro e attualmente mancano di poco per essere magneti "perfettamente duri" (sono "semi-duri").
• Il Futuro: L'autore suggerisce che con un po' più di aggiustamenti (come aggiungere Azoto o modificare la ricetta), questi magneti di Zirconio potrebbero diventare un'alternativa reale ai costosi magneti di Neodimio che usiamo oggi.

In sintesi: L'autore ha utilizzato un computer per dimostrare che possiamo costruire un magnete forte e stabile usando il vicino economico e abbondante (Zirconio) invece del VIP costoso (Neodimio). Non è esattamente forte quanto la versione VIP, ma è abbastanza vicino da poter rivoluzionare il modo in cui produciamo magneti per auto elettriche ed energia verde.

Sommario tecnico

1. Problema e Motivazione

La domanda globale di magneti permanenti ad alte prestazioni è trainata dalle tecnologie di energia rinnovabile e dai veicoli elettrici. Attualmente, l'industria fa molto affidamento sui magneti Nd$_2$Fe$_{14}$B. Tuttavia, il Neodimio (Nd) e altri elementi delle Terre Rare (RE) sono considerati risorse critiche a causa dei monopoli della catena di approvvigionamento (principalmente la Cina) e delle riserve in declino.

• La Sfida: Sviluppare magneti "poveri di RE" o privi di RE che mantengano alte prestazioni magnetiche senza fare affidamento su elementi critici.
• Il Gap Specifico: Sebbene la struttura di tipo ThMn$_{12}$ (1:12) (RFe$_{12}$) offra un contenuto di RE inferiore rispetto a Nd$_2$Fe$_{14}$B, le fasi pure RFe$_{12}$ sono termodinamicamente instabili. Richiedono stabilizzazione (solitamente tramite Titanio, Ti) e spesso soffrono di temperature di Curie o anisotropie inferiori rispetto a Nd$_2$Fe$_{14}$B.
• L'Obiettivo: Questo lavoro investiga lo Zirconio (Zr) come sostituto del Nd nella fase 1:12. Lo Zr è abbondante, economico e non critico. Lo studio mira a determinare se i composti quaternari (Zr,Nd)Fe$_{12-y}$Ti$_y$ e quinary (Zr,Nd)Fe$_{10}$CoTi sostituiti con Zr possano servire come alternative valide e ad alte prestazioni ai magneti Nd-Fe-B.

2. Metodologia

La ricerca impiega metodi computazionali ab initio basati sulla Teoria del Funzionale della Densità (DFT) per prevedere le proprietà magnetiche intrinseche.

• Software e Framework:
  • VASP (Vienna Ab-initio Simulation Package): Utilizzato per il rilassamento strutturale e il calcolo delle energie totali, dei momenti magnetici e dell'anisotropia magnetocristallina.
  • AkaiKKR (MACHIKANEYAMA): Utilizzato per calcolare le temperature di Curie ($T_C$) tramite il metodo della funzione di Green Korringa-Kohn-Rostoker (KKR).
• Approcci Teorici:
  • Funzionali di Scambio-Correlazione: Principalmente l'Approssimazione del Gradiente Generalizzato (GGA-PBE). Per l'energia di anisotropia magnetocristallina (MAE), è stata impiegata anche l'Approssimazione della Densità di Spin Locale (LSDA), poiché spesso fornisce costanti di anisotropia più accurate per questi sistemi.
  • Trattamento degli elettroni f: Gli elettroni 4f del Nd sono stati trattati utilizzando il metodo DFT+U (approccio di Dudarev, $U = 6$ eV) per correggere gli errori di auto-interazione e la forte correlazione elettronica, che la DFT standard non riesce a catturare accuratamente.
  • Pseudopotenziali: Metodo delle Onde Augmentate dal Proiettore (PAW).
• Calcoli Eseguiti:
  • Stabilità di Fase: Calcolo delle energie di formazione e delle entalpie di soluzione per identificare configurazioni stabili di sostituzioni di Ti e Co sui siti Fe e di Zr sui siti RE.
  • Proprietà Magnetiche: Momento magnetico totale ($m_{tot}$), magnetizzazione di saturazione ($\mu_0 M_S$) e prodotto energetico massimo ($|BH|_{max}$).
  • Anisotropia: Calcolo dell'Energia di Anisotropia Magnetocristallina (MAE) ruotando le direzioni di magnetizzazione ([001] a [010]) per derivare la costante di anisotropia ($K_1$) e il campo di anisotropia ($H_a$).
  • Temperatura di Curie ($T_C$): Calcolata utilizzando l'Approssimazione del Campo Medio (MFA) e i modelli di Disordine del Momento Locale (LMD) per tenere conto delle fluttuazioni termiche.
  • Fattore di Durezza ($\kappa$): Calcolato per determinare se il materiale si qualifica come magnete duro ($\kappa > 1$).

3. Contributi Chiave

1. Analisi Sistematica della Sostituzione con Zr: Questo lavoro fornisce la prima indagine teorica completa sulla sostituzione del Nd con Zr nella struttura ThMn$_{12}$ 1:12, andando oltre i precedenti studi focalizzati su Y o Ce.
2. Meccanismi di Stabilizzazione: Conferma che, sebbene il (Zr,Nd)Fe$_{12}$ puro sia instabile, l'aggiunta di Ti (al sito 8i) e Co (ai siti 8f/8j) stabilizza la fase. Stabilisce che la sostituzione con Zr consente una riduzione della concentrazione di Ti rispetto agli analoghi puri a base di Nd.
3. Validazione Metodologica: Lo studio valida l'uso della LSDA per il calcolo dell'anisotropia in questi sistemi, notando che GGA+U spesso sottostima $K_1$ per i composti contenenti Nd, mentre la LSDA fornisce risultati più vicini ai dati sperimentali.
4. Identificazione di Candidati Promettenti: L'identificazione di (Zr,Nd)Fe$_{11}$Ti e (Zr,Nd)Fe$_{10}$CoTi come candidati promettenti che bilanciano costo (Nd ridotto) e prestazioni.

4. Risultati Chiave

A. Stabilità e Struttura

• Stabilizzazione con Ti: Il Ti occupa preferenzialmente il sito cristallografico 8i. Una concentrazione di ~7,7 at.% di Ti (2 atomi per supercella di 26 atomi) è sufficiente per stabilizzare la fase (Zr,Nd)Fe$_{12}$.
• Sostituzione con Co: Il Co preferisce i siti 8f e 8j. Sebbene il Co migliori leggermente i momenti magnetici, ha un effetto stabilizzante inferiore rispetto al Ti.
• Parametri Reticolari: La sostituzione con Zr porta a una riduzione del volume della cella rispetto ai composti puri di Nd a causa del raggio atomico più piccolo dello Zr, seguendo la tendenza: $Nd_2Fe_{24} > ZrNdFe_{24} > Zr_2Fe_{24}$.

B. Momenti Magnetici e Magnetizzazione di Saturazione ($\mu_0 M_S$)

• Effetto dello Zr: Sostituire il Nd con lo Zr riduce il momento magnetico totale perché lo Zr ha un momento magnetico inferiore (-0,47 $\mu_B$) rispetto al Nd (-0,27 $\mu_B$ in GGA, ma -3,3 $\mu_B$ quando gli elettroni f sono inclusi tramite DFT+U).
• Valori Calcolati:
  • NdFe$_{11}$Ti: $\mu_0 M_S \approx 1,69$ T (DFT+U), corrispondente ai range sperimentali (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.
• Conclusione: Sebbene la sostituzione con Zr abbassi leggermente $M_S$, i valori rimangono sufficientemente elevati per applicazioni pratiche.

C. Prodotto Energetico Massimo ($|BH|_{max}$)

• Prestazioni: Il composto quaternario (Zr,Nd)Fe$_{11.5}$Ti$_{0.5}$ mostra un $|BH|_{max}$ teorico di ~600 kJ/m³, che supera il limite teorico di Nd$_2$Fe$_{14}$B (~512 kJ/m³).
• Fase Stabile: Il composto stabile (Zr,Nd)Fe$_{11}$Ti produce ~510–525 kJ/m³, paragonabile a Nd$_2$Fe$_{14}$B e significativamente superiore a Sm$_2$Co$_{17}$ (294 kJ/m³).
• Effetto del Co: La sostituzione con Co aumenta leggermente $M_S$ ma diminuisce marginalmente $|BH|_{max}$ a causa dei cambiamenti di volume.

D. Temperatura di Curie ($T_C$)

• Tendenza: $T_C$ diminuisce con l'aumentare del contenuto di Ti ma aumenta con la sostituzione con Co.
• Valori:
  • 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).
  • Confronto: Questi valori sono più alti di Nd$_2$Fe$_{14}$B ($T_C \approx 588$ K), rendendoli adatti per applicazioni a temperature più elevate.
• Nota: L'approssimazione FM sovrastima $T_C$; i valori LMD sono considerati più realistici ma indicano comunque una temperatura di esercizio vitale.

E. Anisotropia Magnetocristallina ($K_1$) e Fattore di Durezza ($\kappa$)

• Costante di Anisotropia ($K_1$): La sostituzione con Zr aumenta significativamente $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).
  • Significato: Lo Zr potenzia l'anisotropia, una proprietà cruciale per la coercitività.
• Fattore di Durezza ($\kappa$):
  • (Zr,Nd)Fe$_{11}$Ti produce $\kappa \approx 0,97$ (LSDA).
  • Questo classifica il materiale come semi-duro ($\kappa < 1$), mentre Nd$_2$Fe$_{14}$B è duro ($\kappa \approx 1,54$).
  • Il valore è molto vicino alla soglia per i magneti duri, suggerendo che modifiche minori (ad esempio, nitrurazione) potrebbero spingerlo nel regime dei magneti duri.

5. Significato e Prospettive

• Importanza Strategica: Questo lavoro dimostra che lo Zirconio è un sostituto vitale, a basso costo e abbondante per il Neodimio nei magneti di tipo ThMn$_{12}$ senza sacrificare proprietà magnetiche critiche come la magnetizzazione di saturazione o la temperatura di Curie.
• Potenziale di Prestazioni: Le proprietà calcolate suggeriscono che i magneti a base di (Zr,Nd)Fe$_{11}$Ti potrebbero competere con Nd-Fe-B in termini di prodotto energetico e stabilità termica, offrendo al contempo una catena di approvvigionamento più sicura.
• Direzioni Future:
  • Verifica Sperimentale: Le previsioni teoriche richiedono sintesi (ad esempio, tramite Fusione Selettiva Laser) e validazione sperimentale.
  • Ottimizzazione: Poiché il fattore di durezza è leggermente inferiore a 1,0, il lavoro futuro dovrebbe esplorare la nitrurazione (formazione di fasi R-Fe-N) o ulteriore sintonizzazione Co/Ti per spingere $\kappa > 1,0$.
  • Espansione: La metodologia può essere applicata ad altre strutture povere di RE (fasi 2:14, 2:17) o ad altri elementi abbondanti (ad es. La, Ce) per ridurre ulteriormente la dipendenza dalle terre rare critiche.

In conclusione, la tesi identifica con successo i composti (Zr,Nd)Fe$_{12-y}$Ti$_y$ come candidati promettenti per la prossima generazione di magneti permanenti sostenibili e ad alte prestazioni.