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

Resumo Técnico: Capturando Efeitos Quânticos Nucleares na Difusão de Hidrogênio através de MoS₂ via Simulações de Integral de Caminho Aprimoradas por Aprendizado de Máquina

Definição do Problema
O transporte de hidrogênio através de materiais bidimensionais (2D) em camadas, como o dissulfeto de molibdênio (MoS₂), é crítico para aplicações em armazenamento de hidrogênio, células de combustível e separação de isótopos. Embora estudos experimentais tenham demonstrado a difusão interlaminar e o transporte seletivo de isótopos em materiais como MoS₂ e h-BN, descrições teóricas precisas permanecem desafiadoras. Investigações teóricas anteriores basearam-se amplamente em tratamentos clássicos dos núcleos atômicos, negligenciando os Efeitos Quânticos Nucleares (EQNs). No entanto, átomos leves como o hidrogênio exibem comportamentos quânticos pronunciados — incluindo o movimento de ponto zero (ZPE), tunelamento e delocalização — mesmo à temperatura ambiente. Além disso, modelar precisamente esses sistemas requer capturar tanto as interações de van der Waals (vdW) quanto as ligações de hidrogênio, o que necessita de funcionais de troca-correlação avançados (ex: r2SCAN+rVV10) em Teoria do Funcional da Densidade (DFT). O custo computacional de combinar essa DFT de alta fidelidade com Dinâmica de Integral de Caminho (PIMD) e amostragem aprimorada para sistemas grandes e estruturalmente complexos (como bicamadas torcidas) é proibitivo.

Metodologia
Para abordar esses desafios computacionais, os autores desenvolveram um potencial interatômico de aprendizado de máquina (MLIP) baseado na estrutura MACE, treinado em um conjunto de dados compacto gerado usando r2SCAN+rVV10 DFT. Este MLIP foi validado contra cálculos de referência DFT para estabilidade energética e estruturas de bandas de fônons, mostrando boa concordância.

O estudo empregou uma estrutura de simulação híbrida combinando:

1. Metadinâmica Bem Temperada (WTMetaD): Implementada via PLUMED para calcular as barreiras de energia livre ($\Delta F$) para a difusão interlaminar.
2. Dinâmica de Integral de Caminho (PIMD): Implementada via pacote i-PI para incorporar explicitamente os EQNs. Os núcleos foram representados por polímeros de anel com $P=16$ contas para o hidrogênio (e $P=16$ para o deutério, garantindo a convergência) a 300 K.
3. Configurações do Sistema: O estudo investigou politipos de MoS₂ bulk de alta simetria ($H_h^h$, $H_X^h$, $R_M^h$ e fases 3R) e bicamadas de MoS₂ torcidas (TBL-MoS₂) com ângulos de torção de $\theta = 3.89^\circ$, $21.79^\circ$, $38.21^\circ$ e $56.11^\circ$.
4. Análise de Isótopos: Tanto a difusão de prótio (H) quanto a de deutério (D) foram simuladas para quantificar os efeitos de isótopo cinético (KIE).

Principais Contribuições

• Desenvolvimento de uma Estrutura MLIP-PIMD: Os autores estabeleceram um fluxo de trabalho robusto que faz a ponte entre a precisão da DFT e as escalas de comprimento/tempo necessárias para a dinâmica nuclear quântica em sistemas 2D complexos.
• Quantificação de EQNs em MoS₂: O estudo fornece uma avaliação sistemática de como as flutuações quânticas alteram os paisagens de difusão através de diferentes politipos de MoS₂.
• Mecanismo de Seletividade de Isótopos: O trabalho oferece insights atômicos sobre os mecanismos que impulsionam a separação de isótopos, especificamente quantificando a diferença entre as barreiras de energia livre para H e D.
• Paisagens de Difusão Impulsionadas por Moiré: O estudo mapeia o transporte de hidrogênio em bicamadas torcidas, revelando como superredes de moiré criam ambientes de difusão espacialmente variáveis governados pelo empilhamento local (estagiado vs. eclipsado).

Principais Resultados

• Impacto dos EQNs nas Barreiras de Difusão: A inclusão de EQNs reduz substancialmente as barreiras de energia livre para a difusão de hidrogênio em comparação com simulações clássicas. A redução varia de dezenas de meV a quase 60 meV, dependendo da configuração de empilhamento.
  • Para a fase mais estável $H_h^h$, a barreira cai de 122 meV (clássica) para 64 meV (quântica).
  • Para a fase $H_X^h$, a barreira diminui de 45 meV para 18 meV.
  • Esse amolecimento da paisagem de energia livre leva a um aumento significativo no coeficiente de autodifusão ($D$), com Fatores de Correção Quântica (QCF) atingindo até 9.43 para a fase $H_h^h$.
• Efeito de Isótopo Cinético (KIE): Um KIE pronunciado é observado. Embora o deutério exiba barreiras mais altas que o hidrogênio em simulações clássicas, a inclusão de EQNs amplia essa diferença. A barreira de energia livre quântica para o prótio é aproximadamente 35 meV menor que a do deutério (60 meV vs. 95 meV para $H_h^h$), indicando uma forte difusão quântica dependente da massa.
• Dinâmica de Bicamadas Torcidas:
  • Torções de Alto Ângulo ($21.79^\circ, 38.21^\circ$): Estes sistemas carecem de domínios estendidos de alta simetria, resultando em paisagens de difusão dominadas por configurações intermediárias. Consequentemente, as barreiras de difusão são mais altas (~120 meV quânticos) e as difusividades são menores (pelo menos uma ordem de magnitude menores que as fases de alta simetria isoladas).
  • Torções de Baixo Ângulo ($3.89^\circ, 56.11^\circ$): Estes sistemas formam domínios de moiré distintos (estagiados) e nós/solitons (eclipsados). O hidrogênio reside preferencialmente em domínios estagiados, onde as distâncias interlaminares são menores e a coordenação é favorável.
  • Variação Espacial: As barreiras de difusão variam significamente através da superrede de moiré, com diferenças de até 80 meV entre domínios estagiados favoráveis e regiões de soliton/eclipsadas desfavoráveis. O sistema de $56.11^\circ$ exibe barreiras médias mais baixas (150 meV) comparado ao sistema de $3.89^\circ$ (197 meV), atribuído à presença de nós metaestáveis e domínios $H_X^h$ favoráveis.

Significância
O artigo conclui que os EQNs desempenham um papel decisivo no transporte de hidrogênio através de MoS₂ em camadas, mesmo em temperaturas ambientes. Negligenciar esses efeitos leva a uma superestimação das barreiras de difusão e uma subestimação da difusividade. O estudo demonstra que a magnitude do aumento quântico é altamente sensível ao ambiente de empilhamento local. Além disso, a inclusão explícita da dinâmica nuclear quântica é essencial para descrever com precisão os mecanismos de transporte seletivo de isótopos. Finalmente, a pesquisa estabelece que a engenharia de torção em materiais 2D pode ser usada para moldar as paisagens de difusão, criando canais espacialmente seletivos para a migração de hidrogênio. A estrutura MLIP-PIMD apresentada posiciona-se como uma abordagem transferível para investigar o transporte quântico em heteroestruturas de van der Waals complexas relevantes para tecnologias de energia.