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

Esta tesis de licenciatura emplea cálculos de la teoría del funcional de la densidad para investigar las propiedades intrínsecas y magnéticas de compuestos de tipo ThMn$_{12}$ sustituidos con Nd, Zr y Ti, identificando materiales cuaternarios prometedores con bajo contenido de Nd como alternativas potenciales a los imanes de tierras raras escasos.

Lo esencial

La Gran Imagen: Encontrar un Imán Más Barato y Ecológico

Imagina que el mundo moderno es una máquina gigante que funciona con electricidad. Para hacer que esta máquina funcione de manera eficiente (como en los coches eléctricos, las turbinas eólicas y los discos duros), necesitamos imanes potentes. Actualmente, el "estándar de oro" para estos imanes es un material llamado Nd-Fe-B (Neodimio-Hierro-Boro).

Piensa en el Neodimio (Nd) como el invitado VIP de una fiesta. Hace que el imán sea increíblemente fuerte, pero es caro, difícil de encontrar y casi todo el mundo depende de un solo país (China) para abastecerse. Esto crea un cuello de botella en la cadena de suministro, como un único puente estrecho que todos intentan cruzar.

El objetivo de esta investigación es encontrar un sustituto para ese invitado VIP. El autor, Nico Yannik Merkt, sugiere usar Zirconio (Zr). El Zirconio es como el vecino confiable y asequible: es abundante, más barato y más fácil de conseguir. La pregunta es: si cambiamos al VIP por el vecino, ¿sigue funcionando la fiesta (el imán) igual de bien?

El Problema: La "Casa Inestable"

El tipo específico de estructura de imán que se estudia se llama ThMn12 (o la "fase 1:12").

• El Plano: Imagina una casa construida con un plano específico donde tienes 1 átomo de tierra rara (el VIP) y 12 átomos de hierro.
• El Problema: Si intentas construir esta casa solo con el VIP (Neodimio) y los 12 átomos de hierro, la casa es inestable. Es como intentar construir un rascacielos sobre una base de arena; se derrumba.
• La Solución: Para que la casa se mantenga en pie, necesitas un "estabilizador". En este caso, los investigadores usan Titanio (Ti). Piensa en el Titanio como las vigas de acero que añades al marco para evitar que la casa se desmorone.

El Experimento: Una Obra de Construcción Virtual

Dado que construir estos imanes en un laboratorio real es costoso y consume mucho tiempo, el autor utilizó superordenadores para simular la construcción. Esto se llama Teoría del Funcional de la Densidad (DFT).

• La Simulación: En lugar de mezclar químicos en un matraz, el ordenador calcula cómo los átomos "sienten" la presencia de los demás. Pregunta: "Si pongo Zirconio aquí, ¿se mantendrá en pie la casa? ¿Qué tan fuerte será la atracción magnética?"

Lo Que Encontró el Ordenador

El artículo recorre varios escenarios de "qué pasaría si" para ver cómo afecta el intercambio de Neodimio por Zirconio al rendimiento del imán. Aquí están los hallazgos clave:

1. Estabilidad (¿Se mantendrá en pie la casa?)

• Neodimio Puro: Sin ayuda, la casa es inestable.
• Zirconio Puro: Sorprendentemente, una casa hecha enteramente de Zirconio y Hierro sí es estable.
• La Mezcla (50/50): Cuando mezclaron la mitad de Neodimio y la mitad de Zirconio, la casa estaba un poco inestable. Necesitaban añadir más "vigas de acero" (Titanio) para mantenerla estable.
• Conclusión: Puedes reemplazar el Neodimio con Zirconio, pero debes tener cuidado con la receta para mantener la estructura estable.

2. Fuerza (¿Qué tan fuerte es el imán?)

• El Intercambio: El VIP (Neodimio) es naturalmente muy magnético. El vecino (Zirconio) es menos magnético.
• El Resultado: Cuando cambiaron el Neodimio por Zirconio, el imán se debilitó ligeramente. Es como reemplazar una cuerda súper fuerte por una ligeramente más delgada. Sin embargo, el imán sigue siendo muy fuerte: lo suficientemente fuerte para ser útil.
• El "Producto de Energía": Esta es una medida de cuánta energía puede almacenar el imán. Los nuevos imanes basados en Zirconio obtuvieron puntuaciones muy altas, superando a algunos tipos de imanes más antiguos y acercándose al campeón actual (Nd-Fe-B).

3. Resistencia al Calor (La Temperatura de Curie)

• Los imanes pierden su poder si se calientan demasiado. La "Temperatura de Curie" es el punto en el que el imán se rinde y deja de funcionar.
• El Hallazgo: Los nuevos imanes de Zirconio pueden soportar el calor casi tan bien como los de Neodimio. No se derretirán ni perderán su poder en un motor eléctrico caliente.

4. Direccionalidad (La "Calle de Sentido Único")

• Un buen imán permanente necesita ser "duro" para desmagnetizarse. Necesita mantener su dirección firmemente, como una calle de sentido único.
• El Hallazgo: Los imanes de Zirconio son muy buenos manteniendo su dirección. De hecho, en algunos cálculos, los imanes de Zirconio fueron incluso mejores manteniendo su dirección que los de Neodimio.

El Veredicto: ¿Es un Ganador?

El artículo concluye que el Zirconio es un sustituto muy prometedor para el Neodimio.

• Las Ventajas: Es más barato, más abundante y menos crítico para las cadenas de suministro. Los imanes resultantes son estables y tienen excelentes propiedades magnéticas.
• Las Desventajas: Los imanes son ligeramente más débiles que los de Neodimio puro y, actualmente, quedan justo por debajo de ser imanes "perfectamente duros" (son "semi-duros").
• El Futuro: El autor sugiere que con un poco más de ajuste (como añadir Nitrógeno o ajustar la receta), estos imanes de Zirconio podrían convertirse en una alternativa real a los costosos imanes de Neodimio que usamos hoy en día.

En resumen: El autor utilizó un ordenador para demostrar que podemos construir un imán fuerte y estable usando al vecino barato y abundante (Zirconio) en lugar del VIP costoso (Neodimio). No es exactamente tan fuerte como la versión VIP todavía, pero está lo suficientemente cerca como para poder revolucionar la forma en que fabricamos imanes para coches eléctricos y energía verde.

Resumen técnico

A continuación se presenta un resumen técnico detallado de la tesis de licenciatura "Cálculo de las propiedades magnéticas de compuestos cuaternarios de tipo ThMn12 con Zr como sustitución de Nd" de Nico Yannik Merkt.

1. Planteamiento del Problema y Motivación

La demanda global de imanes permanentes de alto rendimiento está impulsada por las tecnologías de energía renovable y los vehículos eléctricos. Actualmente, la industria depende en gran medida de los imanes Nd$_2$Fe$_{14}$B. Sin embargo, el Neodimio (Nd) y otros elementos de Tierras Raras (TR) se consideran recursos críticos debido a los monopolios en la cadena de suministro (principalmente China) y a las reservas en declive.

• El Desafío: Desarrollar imanes "pobres en TR" o libres de TR que mantengan un alto rendimiento magnético sin depender de elementos críticos.
• La Brecha Específica: Aunque la estructura de tipo ThMn$_{12}$ (1:12) (RFe$_{12}$) ofrece un menor contenido de TR que el Nd$_2$Fe$_{14}$B, las fases puras RFe$_{12}$ son termodinámicamente inestables. Requieren estabilización (usualmente mediante Titanio, Ti) y a menudo sufren de temperaturas de Curie o anisotropía más bajas en comparación con el Nd$_2$Fe$_{14}$B.
• El Objetivo: Este trabajo investiga el Zirconio (Zr) como sustituto del Nd en la fase 1:12. El Zr es abundante, económico y no crítico. El estudio busca determinar si los compuestos cuaternarios (Zr,Nd)Fe$_{12-y}$Ti$_y$ y quaternarios (Zr,Nd)Fe$_{10}$CoTi con sustitución de Zr pueden servir como alternativas viables y de alto rendimiento a los imanes Nd-Fe-B.

2. Metodología

La investigación emplea métodos computacionales ab initio basados en la Teoría del Funcional de la Densidad (DFT) para predecir las propiedades magnéticas intrínsecas.

• Software y Marco de Trabajo:
  • VASP (Vienna Ab-initio Simulation Package): Utilizado para la relajación estructural y el cálculo de energías totales, momentos magnéticos y anisotropía magnetocristalina.
  • AkaiKKR (MACHIKANEYAMA): Utilizado para calcular las temperaturas de Curie ($T_C$) mediante el método de la función de Green Korringa-Kohn-Rostoker (KKR).
• Enfoques Teóricos:
  • Funcionales de Intercambio-Correlación: Principalmente la Aproximación del Gradiente Generalizado (GGA-PBE). Para la energía de anisotropía magnetocristalina (MAE), también se empleó la Aproximación de Densidad de Espín Local (LSDA), ya que a menudo produce constantes de anisotropía más precisas para estos sistemas.
  • Tratamiento de electrones f: Los electrones 4f del Nd se trataron utilizando el método DFT+U (enfoque de Dudarev, $U = 6$ eV) para corregir los errores de auto-interacción y la fuerte correlación electrónica, que la DFT estándar no captura con precisión.
  • Pseudopotenciales: Método de Ondas Aumentadas por Proyectores (PAW).
• Cálculos Realizados:
  • Estabilidad de Fase: Cálculo de energías de formación y entalpías de solución para identificar configuraciones estables de sustituciones de Ti y Co en sitios de Fe y de Zr en sitios de TR.
  • Propiedades Magnéticas: Momento magnético total ($m_{tot}$), magnetización de saturación ($\mu_0 M_S$) y producto energético máximo ($|BH|_{max}$).
  • Anisotropía: Cálculo de la Energía de Anisotropía Magnetocristalina (MAE) rotando las direcciones de magnetización ([001] a [010]) para derivar la constante de anisotropía ($K_1$) y el campo de anisotropía ($H_a$).
  • Temperatura de Curie ($T_C$): Calculada utilizando las aproximaciones de Campo Medio (MFA) y Desorden de Momento Local (LMD) para tener en cuenta las fluctuaciones térmicas.
  • Factor de Dureza ($\kappa$): Calculado para determinar si el material califica como un imán duro ($\kappa > 1$).

3. Contribuciones Clave

1. Análisis Sistemático de Sustitución con Zr: Este trabajo proporciona la primera investigación teórica integral sobre la sustitución de Nd por Zr en la estructura ThMn$_{12}$ 1:12, avanzando más allá de estudios anteriores que se centraron en Y o Ce.
2. Mecanismos de Estabilización: Confirma que, aunque el (Zr,Nd)Fe$_{12}$ puro es inestable, la adición de Ti (al sitio 8i) y Co (a los sitios 8f/8j) estabiliza la fase. Establece que la sustitución con Zr permite una reducción en la concentración de Ti en comparación con los análogos basados puramente en Nd.
3. Validación Metodológica: El estudio valida el uso de LSDA para calcular la anisotropía en estos sistemas, señalando que GGA+U a menudo subestima $K_1$ para compuestos que contienen Nd, mientras que LSDA proporciona resultados más cercanos a los datos experimentales.
4. Identificación de Candidatos Prometedores: La identificación de (Zr,Nd)Fe$_{11}$Ti y (Zr,Nd)Fe$_{10}$CoTi como candidatos prometedores que equilibran el costo (Nd reducido) con el rendimiento.

4. Resultados Clave

A. Estabilidad y Estructura

• Estabilización con Ti: El Ti ocupa preferentemente el sitio cristalográfico 8i. Una concentración de ~7.7 % at. de Ti (2 átomos por supercelda de 26 átomos) es suficiente para estabilizar la fase (Zr,Nd)Fe$_{12}$.
• Sustitución con Co: El Co prefiere los sitios 8f y 8j. Aunque el Co mejora ligeramente los momentos magnéticos, tiene un efecto estabilizador menor que el Ti.
• Parámetros de Red: La sustitución con Zr conduce a una reducción del volumen de la celda en comparación con los compuestos puros de Nd debido al menor radio atómico del Zr, siguiendo la tendencia: $Nd_2Fe_{24} > ZrNdFe_{24} > Zr_2Fe_{24}$.

B. Momentos Magnéticos y Magnetización de Saturación ($\mu_0 M_S$)

• Efecto del Zr: Reemplazar Nd con Zr reduce el momento magnético total porque el Zr tiene un momento magnético menor (-0.47 $\mu_B$) en comparación con el Nd (-0.27 $\mu_B$ en GGA, pero -3.3 $\mu_B$ cuando se incluyen electrones f mediante DFT+U).
• Valores Calculados:
  • NdFe$_{11}$Ti: $\mu_0 M_S \approx 1.69$ T (DFT+U), coincidiendo con rangos experimentales (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.
• Conclusión: Aunque la sustitución con Zr disminuye ligeramente $M_S$, los valores permanecen lo suficientemente altos para aplicaciones prácticas.

C. Producto Energético Máximo ($|BH|_{max}$)

• Rendimiento: El compuesto cuaternario (Zr,Nd)Fe$_{11.5}$Ti$_{0.5}$ muestra un $|BH|_{max}$ teórico de ~600 kJ/m³, lo que excede el límite teórico del Nd$_2$Fe$_{14}$B (~512 kJ/m³).
• Fase Estable: El compuesto estable (Zr,Nd)Fe$_{11}$Ti produce ~510–525 kJ/m³, comparable al Nd$_2$Fe$_{14}$B y significativamente mayor que el Sm$_2$Co$_{17}$ (294 kJ/m³).
• Efecto del Co: La sustitución con Co aumenta ligeramente $M_S$ pero disminuye marginalmente $|BH|_{max}$ debido a cambios de volumen.

D. Temperatura de Curie ($T_C$)

• Tendencia: $T_C$ disminuye con el aumento del contenido de Ti pero aumenta con la sustitución de Co.
• Valores:
  • 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).
  • Comparación: Estos valores son más altos que los del Nd$_2$Fe$_{14}$B ($T_C \approx 588$ K), lo que los hace adecuados para aplicaciones a temperaturas más elevadas.
• Nota: La aproximación FM sobreestima $T_C$; los valores LMD se consideran más realistas pero aún indican una temperatura de operación viable.

E. Anisotropía Magnetocristalina ($K_1$) y Factor de Dureza ($\kappa$)

• Constante de Anisotropía ($K_1$): La sustitución 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).
  • Significado: El Zr mejora la anisotropía, una propiedad crucial para la coercitividad.
• Factor de Dureza ($\kappa$):
  • (Zr,Nd)Fe$_{11}$Ti produce $\kappa \approx 0.97$ (LSDA).
  • Esto clasifica al material como semi-duro ($\kappa < 1$), mientras que el Nd$_2$Fe$_{14}$B es duro ($\kappa \approx 1.54$).
  • El valor está muy cerca del umbral para imanes duros, lo que sugiere que modificaciones menores (por ejemplo, nitruración) podrían impulsarlo al régimen de imanes duros.

5. Significado y Perspectivas

• Importancia Estratégica: Este trabajo demuestra que el Zirconio es un sustituto viable, de bajo costo y abundante para el Neodimio en imanes de tipo ThMn$_{12}$ sin sacrificar propiedades magnéticas críticas como la magnetización de saturación o la temperatura de Curie.
• Potencial de Rendimiento: Las propiedades calculadas sugieren que los imanes basados en (Zr,Nd)Fe$_{11}$Ti podrían competir con los Nd-Fe-B en términos de producto energético y estabilidad térmica, al tiempo que ofrecen una cadena de suministro más segura.
• Direcciones Futuras:
  • Verificación Experimental: Las predicciones teóricas requieren síntesis (por ejemplo, mediante Fusión Selectiva por Láser) y validación experimental.
  • Optimización: Dado que el factor de dureza está ligeramente por debajo de 1.0, el trabajo futuro debería explorar la nitruración (formación de fases R-Fe-N) o un ajuste adicional de Co/Ti para impulsar $\kappa > 1.0$.
  • Expansión: La metodología puede aplicarse a otras estructuras pobres en TR (fases 2:14, 2:17) u otros elementos abundantes (por ejemplo, La, Ce) para reducir aún más la dependencia de tierras raras críticas.

En conclusión, la tesis identifica exitosamente a los compuestos (Zr,Nd)Fe$_{12-y}$Ti$_y$ como candidatos prometedores para la próxima generación de imanes permanentes sostenibles y de alto rendimiento.