Enabling Ab-Initio Molecular Dynamics under Bias: The CP2K+SMEAGOL Interface for Integrating Density Functional Theory and Non-Equilibrium Green Functions

Resumo Técnico: Habilitando Dinâmica Molecular Ab-Initio sob Viés via a Interface CP2K+SMEAGOL

Definição do Problema
Embora a Teoria do Funcional da Densidade (DFT) combinada com Funções de Green de Não-Equilíbrio (NEGF) seja uma abordagem padrão para modelar o transporte sob potenciais de viés externos, as metodologias convencionais são limitadas a sistemas estáticos. As implementações existentes operam tipicamente dentro de um arcabouço canônico, calculando propriedades de transporte em posições atômicas fixas. Isso cria uma lacuna significativa entre a compreensão teórica e os resultados macroscópicos de experimentos, particularmente para sistemas eletroquímicos (EC) onde as transformações são intrinsecamente grand canônicas. Para descrever com precisão o ambiente natural das transformações EC, as simulações requerem dinâmica molecular (MD) sob viés e uma descrição de contorno aberto explícita de elétrons que possam entrar e sair da célula computacional. Além disso, as implementações atuais de DFT-NEGF muitas vezes carecem da capacidade de calcular forças induzidas pela corrente necessárias para otimização de geometria ou dinâmica, ou são restritas a códigos DFT específicos (por exemplo, SIESTA) que carecem de recursos presentes em pacotes de última geração otimizados para simulações de fase condensada.

Metodologia
Os autores apresentam uma nova interface entre o popular pacote de estrutura eletrônica de mistura Gaussiana/onda plana CP2K e o código de NEGF SMEAGOL. Esta interface, denominada CP2K+SMEAGOL, foi projetada para ser uma biblioteca independente capaz de ser vinculada a qualquer software DFT, reutilizando o código original do SMEAGOL o máximo possível.

As principais implementações técnicas incluem:

• Forças Induzidas pela Corrente: A interface implementa o cálculo de forças não conservativas atuando nos íons sob viés. Essas forças são derivadas da derivada temporal dos momentos atômicos, utilizando os gradientes do Hamiltoniano completo em relação às posções atômicas. Crucialmente, a implementação inclui o termo de força decorrente da derivada da matriz de densidade ($\Omega_{ij}$), que contabiliza o componente de "vento eletrônico" (transferência de momento de elétrons para íons). Este termo é frequentemente negligenciado em outras implementações, mas mostra-se essencial para a correção qualitativa.
• Construção da Matriz de Densidade: Para ocupar corretamente os estados ligados dentro da janela de viés, a interface atualiza o cálculo da matriz de densidade. Em vez de usar uma única avaliação de contorno ou uma simples média de contribuições esquerda/direita, ela emprega um esquema de duplo contorno ponderado (seguindo Brandbyge et al.) que soma integrais de ambos os eletrodos com pesos determinados pelas contribuições de não-equilíbrio.
• Configuração da Simulação: O método permite a dinâmica molecular de Born-Oppenheimer sob viés. O sistema é modelado com uma região de "molécula estendida" central conectada a eletrodos semi-infinitos. As auto-energias dos eletrodos são tratadas como estáticas e independentes da densidade de carga da região central, uma suposição válida quando a região central é suficientemente grande para blindar mudanças de carga.

Principais Resultados e Validação
Os autores validaram a interface CP2K+SMEAGOL contra cálculos de referência realizados com SIESTA+SMEAGOL através de diversos sistemas:

1. Forças de Viés Zero (Fio de Au Infinito): A interface reproduz com sucesso as forças atômicas em viés zero, coincidindo com cálculos apenas de CP2K. Ela também demonstra que negligenciar o termo da matriz de densidade de energia ($\Omega$) leva a forças qualitativamente incorretas.
2. Capacitor de Placas Paralelas: A interface reproduz a queda de potencial esperada através de um capacitor, mostrando um potencial constante nos eletrodos e uma queda linear no vácuo. Confirma que o esquema de duplo contorno ponderado é necessário para produzir densidades de carga simétricas e potenciais de Hartree corretos em alto viés, enquanto avaliações de contorno único produzem resultados assimétricos e incorretos.
3. Otimização de Geometria da Junção Au-H$_2$-Au: Sob viés aplicado (até 1,5 V), a interface prevê corretamente o alongamento da ligação H-H e os deslocamentos atômicos, consistentes com resultados anteriores de SIESTA+SMEAGOL. As mudanças são atribuídas à corrente elétrica, em vez do campo elétrico.
4. Dinâmica Molecular de um Fio de Au Solvatado: Os autores realizaram as primeiras simulações de dinâmica molecular DFT-NEGF em um sistema de fase condensada de grande escala (um fio monoatômico de Au solvatado por 166 moléculas de água).
   • Estabilidade: A dinâmica foi estável com deriva de energia mínima a longo prazo, apesar da natureza de sistema aberto.
   • Migração de Elétrons: Em uma trajetória de 200 fs a 1 V de viés, a simulação observou efeitos de migração de elétrons, onde o mínimo da densidade eletrônica local deslocou-se do átomo central do fio para o átomo à direita, consistente com a blindagem da carga positiva do eletrodo direito.
   • Transmissão: Os espectros de transmissão para o fio solvatado mostraram similaridade qualitativa com um fio infinito ideal, com diferenças atribuídas a canais de tunelamento adicionais abertos pelas moléculas de água.

Desempenho
O principal gargalo computacional é a avaliação da densidade eletrônica como uma integral da função de Green, exigindo 96 computações de pontos de energia (64 no eixo real, 32 no eixo complexo) por passo de SCF, substituindo uma única diagonalização de Hamiltoniano.

• Paralelização: O código utiliza MPI para computações no espaço real e pontos k, e OpenMP para loops DO e rotinas LAPACK. O desempenho escala bem até 32 threads OpenMP.
• Custo: Para um sistema de fio de Au solvatado, o cálculo CP2K+SMEAGOL foi aproximadamente 3,2 vezes mais lento por passo de SCF do que um cálculo CP2K padrão, apesar de utilizar 8 vezes mais núcleos (2048 vs. 256). Os autores observam que isso representa um aumento de uma ordem de magnitude no custo, tornando trajetórias longas de MD desafiadoras sem aceleração adicional (por exemplo, via ScaLAPACK ou campos de força de aprendizado de máquina).

Significância e Alegações
O artigo afirma que o CP2K+SMEAGOL é a implementação DFT-NEGF mais rica em recursos e completa disponível para o CP2K até o momento, superando a implementação nativa do CP2K, que é limitada ao transporte no ponto $\Gamma$ sem forças.

A principal significância deste trabalho é a demonstração de que a DFT-NEGF pode ser usada para realizar simulações de dinâmica molecular sob viés para sistemas de fase condensada de grande escala sob condições operacionais realistas. Tanto quanto os autores sabem, estas representam as primeiras tais cálculos realizados. A interface permite o estudo de reações eletrocatalíticas e sistemas eletroquímicos onde a interação entre o potencial aplicado, o fluxo de corrente e a dinâmica atômica é crítica, preenchendo a lacuna entre modelos teóricos microscópicos e condições experimentais macroscópicas.