Many-body post-processing of density functional calculations using the variational quantum eigensolver for Bader charge analysis

Título: Pós-processamento de muitos corpos de cálculos de densidade funcional usando o Variational Quantum Eigensolver para análise de carga de Bader

Autores: Erik Schultheis, Alexander Rehn e Gabriel Breuil (German Aerospace Center - DLR).

---

1. O Problema

A determinação precisa da distribuição de carga eletrônica em materiais é crucial para prever propriedades industriais e comportamentos químicos. Uma métrica comum para isso é a Carga de Bader, que quantifica a carga atômica baseada na densidade eletrônica no espaço real.

No entanto, métodos padrão de Teoria do Funcional da Densidade (DFT), como o Kohn-Sham (KS-DFT), enfrentam desafios significativos:

• Correlação Eletrônica Forte: Materiais com elétrons fortemente correlacionados (ex.: óxidos de metais de transição) são mal descritos pelo DFT padrão devido a erros de auto-interação e à aproximação do funcional de troca-correlação (XC).
• Limitações do DFT+U: Embora o método DFT+U (DFT mais um termo de Hubbard) melhore a descrição desses sistemas, ele requer o ajuste manual e não transferível do parâmetro $U$, o que é computacionalmente custoso e subjetivo.
• Custo Computacional: Métodos de muitos corpos clássicos (como Coupled-Cluster ou CI) que poderiam resolver o problema com maior precisão tornam-se proibitivamente caros para sistemas periódicos devido à escala exponencial ou de alta potência polinomial.

O objetivo deste trabalho é superar essas limitações utilizando algoritmos quânticos híbridos para refinar a análise de carga de Bader em sistemas correlacionados.

---

2. Metodologia

Os autores desenvolveram um fluxo de trabalho de pós-processamento que combina cálculos DFT clássicos com simulações quânticas de muitos corpos. O framework é implementado no pacote de software de código aberto Dopyqo.

O processo segue os seguintes passos:

1. Cálculo DFT Inicial:

   • Realiza-se um cálculo DFT padrão (usando o código Quantum ESPRESSO) para obter os orbitais de Kohn-Sham (KS) e a estrutura eletrônica básica.
   • Utiliza-se pseudopotenciais conservadores de norma (NC) e cálculos no ponto $\Gamma$ (ou malhas k, dependendo do sistema).

2. Construção do Hamiltoniano de Muitos Corpos:

   • A partir dos orbitais KS, constrói-se um Hamiltoniano de muitos corpos em um espaço ativo selecionado.
   • O espaço ativo é definido pelos orbitais próximos à energia de Fermi (geralmente orbitais $d$ ou $f$ de metais de transição), onde a correlação é forte.
   • Aplica-se a aproximação de núcleo congelado (frozen core), onde os elétrons internos são tratados implicitamente, reduzindo o número de elétrons ativos a serem simulados.
   • Os elementos de matriz do Hamiltoniano (cinético, potencial pseudopotencial e interação Coulombiana) são calculados a partir dos orbitais KS e transformados para a base de segundas quantizações.

3. Solução com VQE (Variational Quantum Eigensolver):

   • O Hamiltoniano do espaço ativo é resolvido usando o algoritmo VQE, um algoritmo híbrido quântico-clássico.
   • Utiliza-se um ansatz UCCSD (Unitary Coupled Cluster com Singles and Doubles) que preserva o número de partículas.
   • O estado fundamental quântico é otimizado classicamente para minimizar a energia.

4. Cálculo da Densidade e Cargas de Bader:

   • A partir do estado fundamental otimizado pelo VQE, calcula-se a densidade eletrônica de muitos corpos.
   • A densidade é mapeada em uma grade no espaço real.
   • Um algoritmo especializado integra a densidade dentro das regiões de Bader (delimitadas por superfícies de fluxo zero) para obter as Cargas de Bader (ou Cargas Excessivas de Bader - BECs).

---

3. Contribuições Principais

• Framework Dopyqo: Desenvolvimento e disponibilização de um pacote de software de código aberto que automatiza a construção de Hamiltonianos de muitos corpos a partir de saídas do Quantum ESPRESSO e sua resolução via VQE.
• Método Híbrido Eficiente: Demonstração de como usar recursos quânticos (simulados) para corrigir falhas do DFT em sistemas correlacionados sem a necessidade de parâmetros empíricos como o $U$ do DFT+U.
• Análise de Carga de Bader Quântica: Primeira aplicação sistemática de simulações de muitos corpos via VQE para calcular cargas de Bader em sistemas periódicos, fornecendo uma nova perspectiva sobre a localização de carga.

---

4. Resultados

Os autores validaram o método em dois conjuntos de materiais:

1. Sistemas de MgH2 (Correlação Fraca):

   • Supercélulas de MgH2 dopadas com metais de transição (Ti, Cr, Fe, Ni, Zn).
   • Resultado: As cargas de Bader calculadas pelo método Dopyqo (VQE) foram muito próximas das obtidas pelo DFT padrão e dos dados de referência.
   • Conclusão: Para sistemas fracamente correlacionados, o DFT padrão já é suficiente, e o método quântico não introduz desvios significativos, validando a precisão do código.

2. Óxidos de Metais de Transição (Correlação Forte):

   • Materiais como $CrO_2$, $RhO_2$, $RuO_2$, $MoO_2$, $TiO_2$ e $TiS_2$.
   • Desempenho do DFT: O DFT padrão (NC) falhou em reproduzir as cargas de referência (obtidas via DFT+U ou experimentos indiretos) para sistemas com forte correlação (ex.: $CrO_2$, $TiS_2$), subestimando as cargas nos átomos metálicos.
   • Desempenho do Dopyqo: O método VQE conseguiu recuperar as cargas de Bader com alta precisão, aproximando-se dos valores de referência do DFT+U.
   • Comparação: Em casos como $CrO_2$, o DFT+U exigiu um ajuste de parâmetro $U$, enquanto o Dopyqo obteve resultados superiores ao DFT padrão sem ajuste de parâmetros, capturando a correlação eletrônica intrinsecamente através da solução do Hamiltoniano de muitos corpos.

---

5. Significado e Impacto

• Superação de Limitações do DFT: O trabalho demonstra que é possível obter descrições de alta fidelidade de materiais fortemente correlacionados usando uma abordagem de pós-processamento quântico, evitando os erros de auto-interação do DFT e a subjetividade do DFT+U.
• Viabilidade Prática: Embora os cálculos de muitos corpos sejam mais custosos que o DFT, o uso de um espaço ativo reduzido torna o problema tratável para hardware quântico atual (ou simuladores de alta fidelidade), oferecendo uma via para materiais complexos.
• Fluxo de Trabalho Automatizado: A ferramenta Dopyqo permite que pesquisadores integrem facilmente simulações quânticas em seus fluxos de trabalho de descoberta de materiais existentes, sem a necessidade de reescrever códigos de DFT.
• Futuro: O método sugere que, à medida que o hardware quântico amadurecer, será possível realizar esses cálculos em computadores quânticos reais para prever propriedades de materiais industriais com precisão sem precedentes, especialmente para catalisadores e baterias onde a transferência de carga é crítica.

Em resumo, o artigo estabelece um marco importante na aplicação prática de algoritmos quânticos para a química computacional de materiais, focando especificamente na melhoria da precisão na análise de cargas atômicas em sistemas desafiadores.