Stoichiometry Dependent Properties of Cerium Hydride: An Active Learning Developed Interatomic Potential Study

Resumen Técnico: Propiedades del Hidruro de Cerio mediante Potenciales Interatómicos de Aprendizaje Automático

1. Planteamiento del Problema

El hidruro de cerio ($CeH_X$) posee propiedades físicas fascinantes, incluyendo una contracción de red y un aumento de densidad a medida que aumenta el contenido de hidrógeno. Sin embargo, el control preciso de la estequiometría es experimentalmente difícil, y los métodos ab initio (como la Teoría del Funcional de la Densidad, DFT) son computacionalmente prohibitivos para estudiar propiedades que requieren escalas de tiempo y longitud grandes, como la fusión y la difusión a bajas temperaturas. Además, la falta de datos experimentales densos sobre cómo varían las propiedades con la estequiometría limita la comprensión de los mecanismos fundamentales.

2. Metodología

Los autores desarrollaron un Potencial Interatómico de Aprendizaje Automático (MLIAP) utilizando una estrategia de Aprendizaje Activo (AL) basada en el método "Query-by-Committee" (Consulta por Comité).

• Modelo: Se utilizó una Red Neuronal de Partículas que Interactúan Jerárquicamente (HIPNN). El modelo consta de 2 capas de interacción y 3 capas de átomos con 24 funciones de sensibilidad, utilizando un corte de 5.5 Å.
• Generación de Datos (Active Learning):
  • Se entrenó un comité de 8 modelos HIPNN con divisiones aleatorias de entrenamiento/prueba.
  • Se ejecutaron simulaciones de Dinámica Molecular (MD) clásica para explorar el espacio de fases (temperaturas de 200 K a 4500 K y densidades de 0.5 a 18 g/cm³).
  • Criterio de Incertidumbre: Si la desviación estándar de la fuerza o energía predicha por el comité superaba umbrales (0.2 eV/Å y 0.5 eV, respectivamente), la estructura se enviaba a cálculos DFT para etiquetarla.
• Cálculos de Referencia (DFT): Se utilizaron cálculos DFT al nivel PBE+U (con parámetro Hubbard $U = 4.0$ eV) mediante el código VASP. Se utilizaron 18,068 estructuras en total para el conjunto de datos final.
• Validación: El modelo final mostró errores medios absolutos (MAE) de 0.18 eV/Å para fuerzas y 0.0055 eV para energías.

3. Contribuciones Clave

• Desarrollo de un MLIAP robusto: Creación de un potencial válido para hidruros de cerio con relaciones H/Ce entre 2.0 y 3.0, capaz de cubrir un amplio rango de temperaturas y densidades.
• Mapeo de Propiedades Estequiométricas: Evaluación sistemática de propiedades mecánicas, térmicas y de transporte en una resolución estequiométrica sin precedentes, superando las limitaciones de escalabilidad de la DFT.
• Descubrimiento de Mecanismos Físicos: Identificación de los mecanismos microscópicos que gobiernan las tendencias macroscópicas, específicamente el papel de los sitios octaédricos y los campos de deformación.

4. Resultados Principales

A. Contracción de la Red y Densidad:

• Se confirma que la densidad aumenta y el parámetro de red se contrae al añadir hidrógeno en sitios octaédricos.
• Esta relación es no lineal: la contracción es rápida al inicio, pero se atenúa (tasa decreciente) a medida que la ocupación de sitios octaédricos se acerca a la saturación (3 átomos de H por Ce).
• Mecanismo: La superposición de campos de deformación (strain fields) a altas ocupaciones causa que los efectos adicionales de contracción se cancelen mutuamente.

B. Propiedades Mecánicas (Constantes Elásticas):

• $C_{11}$: Aumenta fuertemente hasta $CeH_{2.5}$ y luego muestra un máximo local con una ligera caída hacia $CeH_3$.
• $C_{12}$: Muestra un aumento lineal directo con el contenido de hidrógeno, indicando una mayor rigidez ante deformaciones laterales.
• $C_{44}$ (Resistencia al corte): Sigue la tendencia de contracción de la red, aumentando con el hidrógeno pero con una tasa decreciente similar a la densidad.
• Respuesta P-V: A bajas presiones, la respuesta es similar entre estequiometrías, pero a altas presiones (>30 GPa), las concentraciones más altas de hidrógeno muestran una resistencia a la compresión significativamente mayor debido al endurecimiento de la red.

C. Temperatura de Fusión:

• La temperatura de fusión aumenta drásticamente al principio con el contenido de hidrógeno y luego se estabiliza.
• Esto se debe al fortalecimiento del enlace de la red y al aumento de la energía cohesiva.
• Limitación: No se pudo predecir la fusión para $X > 2.7$ debido a la separación de fases (el líquido tiende a separarse en Ce y $H_2$) en las escalas de tiempo de simulación.

D. Difusión de Hidrógeno:
Se observaron dos regímenes de temperatura distintos con comportamientos opuestos:

1. Altas Temperaturas: La difusividad es máxima en $CeH_2$ y disminuye al añadir hidrógeno octaédrico. Esto sigue la difusión clásica de red, donde los sitios octaédricos vacíos actúan como vacantes necesarias para el movimiento.
2. Bajas Temperaturas: Se observa un máximo local de difusividad en estequiometrías intermedias (entre 2.5 y 2.6 H/Ce).
   • Mecanismo Propuesto: La ocupación de sitios octaédricos reduce la barrera de activación para la difusión tetraédrica. Esto podría deberse a un mecanismo concertado (empuje de átomos) o a la distorsión de la red que debilita la unión de los átomos de hidrógeno en sitios tetraédricos, facilitando su movimiento.

5. Significado e Impacto

Este estudio demuestra la viabilidad de utilizar potenciales interatómicos guiados por aprendizaje activo para explorar sistemas complejos de materiales donde la experimentación es difícil y la simulación ab initio es inviable.

• Unificación de Tendencias: La mayoría de las propiedades (elasticidad, fusión, rigidez) están gobernadas por la contracción de la red y el fortalecimiento del enlace inducido por el hidrógeno octaédrico.
• Excepción Crítica: La difusión es la única propiedad que no sigue esta tendencia monótona, revelando una complejidad dinámica donde la ocupación de sitios altera las barreras de energía de activación de manera no trivial.
• Aplicabilidad: El modelo desarrollado permite predecir el comportamiento del hidruro de cerio en condiciones extremas (altas temperaturas y presiones), lo cual es crucial para aplicaciones en almacenamiento de hidrógeno y materiales nucleares.

En conclusión, el trabajo no solo proporciona una herramienta computacional precisa para el hidruro de cerio, sino que también desentraña los mecanismos físicos fundamentales que vinculan la microestructura atómica con las propiedades macroscópicas en función de la estequiometría.