Capturing Nuclear Quantum Effects in Hydrogen Diffusion through MoS2 via Machine-Learning-Enhanced Path-Integral Simulations

Resumen Técnico: Captura de Efectos Cuánticos Nucleares en la Difusión de Hidrógeno a través de MoS₂ mediante Simulaciones de Integrales de Trayectoria Mejoradas por Aprendizaje Automático

Planteamiento del Problema
El transporte de hidrógeno a través de materiales bidimensionales (2D) estratificados, como el disulfuro de molibdeno (MoS₂), es crítico para aplicaciones en almacenamiento de hidrógeno, celdas de combustible y separación de isótopos. Aunque los estudios experimentales han demostrado la difusión interlaminar y el transporte selectivo de isótopos en materiales como el MoS₂ y el h-BN, las descripciones teóricas precisas siguen siendo un desafío. Las investigaciones teóricas previas se basaron mayoritariamente en tratamientos clásicos de los núcleos atómicos, omitiendo los Efectos Cuánticos Nucleares (NQE). Sin embargo, átomos ligeros como el hidrógeno exhiben comportamientos cuánticos pronunciados —incluyendo el movimiento de punto cero (ZPE), el efecto túnel y la deslocalización— incluso a temperatura ambiente. Además, modelar con precisión estos sistemas requiere capturar tanto las interacciones de van der Waals (vdW) como el enlace de hidrógeno, lo que exige funcionales de intercambio y correlación avanzados (por ejemplo, r2SCAN+rVV10) en la Teoría del Funcional de la Densidad (DFT). El costo computacional de combinar esta DFT de alta fidelidad con la Dinámica de Trayectorias de Punto Integral (PIMD) y el muestreo mejorado para sistemas grandes y estructuralmente complejos (como las bicapas retorcidas) es prohibitivo.

Metodología
Para abordar estos desafíos computacionales, los autores desarrollaron un potencial interatómico basado en aprendizaje automático (MLIP) basado en el marco MACE, entrenado con un conjunto de datos compacto generado mediante r2SCAN+rVV10 DFT. Este MLIP fue validado frente a cálculos de referencia DFT para estabilidad energética y estructuras de bandas de fonones, mostrando una buena concordancia.

El estudio empleó un marco de simulación híbrido que combina:

1. Metadinámica de Temperatura Bien Temperada (WTMetaD): Implementada mediante PLUMED para calcular las barreras de energía libre ($\Delta F$) para la difusión interlaminar.
2. Dinámica de Trayectorias de Punto Integral (PIMD): Implementada mediante el paquete i-PI para incorporar explícitamente los NQE. Los núcleos se representaron mediante polímeros de anillo con $P=16$ cuentas para el hidrógeno (y $P=16$ para el deuterio, asegurando la convergencia) a 300 K.
3. Configuraciones del Sistema: El estudio investigó politipos de MoS₂ masivos de alta simetría ($H_h^h$, $H_X^h$, $R_M^h$ y fases 3R) y bicapas de MoS₂ retorcidas (TBL-MoS₂) con ángulos de torsión de $\theta = 3.89^\circ$, $21.79^\circ$, $38.21^\circ$ y $56.11^\circ$.
4. Análisis de Isótopos: Se simularon tanto la difusión de protio (H) como la de deuterio (D) para cuantificar los efectos de isótopo cinético (KIE).

Contribuciones Clave

• Desarrollo de un Marco MLIP-PIMD: Los autores establecieron un flujo de trabajo robusto que cierra la brecha entre la precisión de DFT y las escalas de longitud/tiempo requeridas para la dinámica nuclear cuántica en sistemas 2D complejos.
• Cuantificación de los NQE en MoS₂: El estudio proporciona una evaluación sistemática de cómo las fluctuaciones cuánticas alteran los paisajes de difusión a través de diferentes politipos de MoS₂.
• Mecanismo de Selectividad de Isótopos: El trabajo ofrece una visión atómica de los mecanismos que impulsan la separación de isótopos, cuantificando específicamente la diferencia entre las barreras de energía libre para H y D.
• Paisajes de Difusión Impulsados por Moiré: El estudio mapea el transporte de hidrógeno en bicapas retorcidas, revelando cómo los superredes de moiré crean entornos de difusión espacialmente variables gobernados por el apilamiento local (escalonado frente a eclipsado) y redes de solitones.

Resultos Clave

• Impacto de los NQE en las Barreras de Difusión: La inclusión de los NQE reduce sustancialmente las barreras de energía libre para la difusión de hidrógeno en comparación con las simulaciones clásicas. La reducción oscila desde decenas de meV hasta casi 60 meV dependiendo de la configuración de apilamiento.
  • Para la fase más estable $H_h^h$, la barrera cae de 122 meV (clásica) a 64 meV (cuántica).
  • Para la fase $H_X^h$, la barrera disminuye de 45 meV a 18 meV.
  • Este ablandamiento del paisaje de energía libre conduce a un aumento significativo en el coeficiente de autodifusión ($D$), con Factores de Corrección Cuántica (QCF) que alcanzan hasta 9.43 para la fase $H_h^h$.
• Efecto de Isótopo Cinético (KIE): Se observa un KIE pronunciado. Mientras que el deuterio exhibe barreras más altas que el hidrógeno en simulaciones clásicas, la inclusión de los NQE amplifica esta diferencia. La barrera de energía libre cuántica para el protio es aproximadamente 35 meV menor que la del deuterio (60 meV frente a 95 meV para $H_h^h$), lo que indica una fuerte difusión cuántica dependiente de la masa.
• Dinámica de Bicapas Retorcidas:
  • Torsiones de Ángulo Alto ($21.79^\circ, 38.21^\circ$): Estos sistemas carecen de dominios extendidos de alta simetría, lo que resulta en paisajes de difusión dominados por configuraciones intermedias. Consecuentemente, las barreras de difusión son más altas (~120 meV cuánticos) y las difusividades son menores (al menos un orden de magnitud menores que las fases de alta simetría aisladas).
  • Torsiones de Ángulo Bajo ($3.89^\circ, 56.11^\circ$): Estos sistemas forman dominios de moiré distintos (escalonados) y nodos/solitones (eclipsados). El hidrógeno reside preferentemente en los dominios escalonados donde las distancias interlaminares son menores y la coordinación es favorable.
  • Variación Espacial: Las barreras de difusión varían significamente a través de la superred de moiré, con diferencias de hasta 80 meV entre los dominios escalonados favorables y las regiones de solitones/eclipsadas desfavorables. El sistema de $56.11^\circ$ exhibe barreras promedio más bajas (150 meV) en comparación con el de $3.89^\circ$ (197 meV), debido a la presencia de nodos metaestables y dominios $H_X^h$ favorables.

Significancia
El artículo concluye que los NQE juegan un papel decisivo en el transporte de hidrógeno a través de MoS₂ estratificado, incluso a temperaturas ambiente. Omitir estos efectos conduce a una sobreestimación de las barreras de difusión y a una subestimación de la difusividad. El estudio demuestra que la magnitud de la mejora cuántica es altamente sensible al entorno de apilamiento local. Además, la inclusión explícita de la dinámica nuclear cuántica es esencial para describir con precisión los mecanismos de transporte selectivo de isótopos. Finalmente, la investigación establece que la ingeniería de torsión en materiales 2D puede utilizarse para moldear los paisajes de difusión, creando canales espacialmente selectivos para la migración de hidrógeno. El marco MLIP-PIMD presentado se posiciona como un enfoque transferible para investigar el transporte cuántico en heteroestructuras de van der Waals complejas relevantes para tecnologías de energía.