Efficient Quantum Implementation of Dynamical Mean Field Theory for Correlated Materials

Resumo Técnico: Implementação Quântica Eficiente da Teoria do Campo Médio Dinâmico para Materiais Correlacionados

Enunciação do Problema
A descrição teórica precisa de materiais com elétrons fortemente correlacionados permanece um desafio formidável na física da matéria condensada e na química computacional. Métodos de incorporação, particularmente a Teoria do Campo Médio Dinâmico (DMFT), abordam isso mapeando modelos de rede complexos em problemas de impureza efetivos acoplados a um banho. No entanto, a aplicação prática da DMFT é limitada pelo custo computacional de resolver o modelo de impureza, especificamente o cálculo da função de Green (GF) da impureza. Solucionadores tradicionais enfrentam obstáculos significativos: o Monte Carlo Quântico (QMC) sofre com problemas de sinal, a Diagonalização Exata (ED) é restringida pelo crescimento exponencial do espaço de Hilbert, e os métodos de rede de tensor apenas aliviam parcialmente essas limitações. Embora a computação quântica ofereça um caminho potencial, as abordagens variacionais existentes são frequentemente afetadas por platôs áridos e ruído de hardware, e a evolução temporal tipicamente requer circuitos profundos com milhares de portas de dois qubits, excedendo as capacidades dos atuais dispositivos quânticos de escala intermediária ruidosa (NISQ).

Metodologia
Este trabalho propõe um framework híbrido quântico-clássico para DMFT que aproveita circuitos de baixa profundidade e representações eficientes de estados para superar as limitações atuais de hardware. A metodologia é construída sobre três pilares centrais:

1. Representação de Subespaço via Estados Gaussianos Fermiônicos (FGS):
   Em vez de preparação variacional de estado, o estado fundamental da impureza é aproximado como uma combinação linear de Estados Gaussianos Fermiônicos (SGS). Esta abordagem permite:

   • Diagonalização Clássica de Subespaço: Os coeficientes do estado fundamental são determinados classicamente diagonalizando o Hamiltoniano dentro de um subespaço de FGS.
   • Continuação de Autovetores: Um único subespaço $S$, construído a partir de um conjunto candidato de FGS, pode ser reutilizado em múltiplas iterações da DMFT à medida que os parâmetros do banho variam, desde que as mudanças sejam graduais.
   • Preparação Quântica Exata: Os FGS podem ser preparados exatamente e eficientemente em hardware quântico usando operadores de evolução fermiônica livre.

2. Compressão Parcial de Circuitos:
   Para calcular a GF da impureza, o framework avalia funções de correlação evoluídas no tempo. Os autores empregam técnicas de compressão algébrica de circuitos para reduzir significativamente a profundidade do circuito:

   • Exploração de Matchgates: O banho é não interativo (fermiônico livre), permitindo o uso de matchgates (equivalentes a rotações de Givens).
   • Regras de Compressão: Ao utilizar fusão, comutação e propriedades de turnover dos matchgates, os autores comprimem os circuitos de evolução temporal e preparação de estado.
   • Simplificação Estrutural: A estrutura do teste de Hadamard é simplificada deslocando a preparação de estado para o início do circuito e absorvendo partes da evolução de Trotter na preparação de estado ou descartando portas que caem fora do cone de luz da medição. Isso reduz a escala da contagem de CNOT de $O(N_q^2)$ por passo de Trotter para $O(N_q(N_q + r N_I))$, onde $N_q$ é o número total de qubits, $r$ é o número de passos de Trotter e $N_I$ é o número de orbitais de impureza.

3. Mitigação de Erros e Processamento de Sinal:
   Reconhecendo o ruído inerente ao hardware atual, o framework integra pós-processamento clássico:

   • Mitigação de Hardware: As técnicas incluem twirling de Pauli (compilação aleatorizada), desacoplamento dinâmico (sequências de pulsos XY) e pós-seleção baseada na conservação do número de partículas.
   • Desruído e Extensão via PSD: A função de correlação é tratada como uma função semidefinida positiva (PSD). Dados ruidosos são projetados na matriz PSD mais próxima para remover ruído, e o sinal é estendido no domínio do tempo usando extensão PSD. Isso aprimora a resolução de frequência sem exigir recursos quânticos adicionais, permitindo a extração de funções de Green de Matsubara via transformada de Fourier.

Resultados Chave
O artigo demonstra a viabilidade desta abordagem através de simulações livres de ruído e experimentos em hardware:

• Convergência da Base: Em simulações livres de ruído, a representação SGS mostrou-se capaz de reproduzir fielmente o estado fundamental da impureza e a função de Green. O posto $\chi$ necessário da SGS foi encontrado como uma fração muito pequena da dimensão do espaço de Hilbert selecionado por partículas (por exemplo, $<1\%$ para sistemas com até 7 orbitais de banho), mesmo à medida que o tamanho do sistema aumentava.
• Autoconsistência da DMFT: O framework convergiu com sucesso o loop da DMFT para um modelo de impureza única em uma rede de Bethe. A densidade de estados (DOS) resultante e o peso de quasipartícula ($Z$) corresponderam estreitamente às referências de ED através da transição metal-isolante, validando que o subespaço SGS permanece preciso ao longo do processo iterativo da DMFT.
• Demonstração em Hardware: Os autores executaram o algoritmo no processador quântico Sherbrooke da IBM para um sistema com um orbital de impureza e três orbitais de banho (8 qubits + 1 ancilla).
  • Eles extraíram com sucesso a função de Green da impureza para uma fase isolante ($U=5.337$).
  • Os dados de hardware mitigados de erro, após desruído e extensão via PSD, recuperaram as características de frequência proeminentes da solução exata.
  • Embora artefatos espúrios de alta frequência tenham aparecido devido ao ruído, eles foram identificados como estando fora da largura de banda física e puderam ser descartados.

Significância e Alegações
Os autores afirmam que este trabalho estabelece um caminho prático para o uso de computadores quânticos como solucionadores de impureza na DMFT para materiais correlacionados. A significância reside em:

• Viabilidade de Hardware: Ao combinar SGS com compressão parcial de circuitos, o método reduz as contagens de portas para níveis gerenciáveis pelos atuais dispositivos NISQ, evitando os circuitos profundos tipicamente necessários para evolução temporal.
• Robustez: A integração de desruído e extensão via PSD permite a extração de informações dinâmicas de dados ruidosos, preenchendo a lacuna entre as capacidades atuais de hardware e a precisão requerida para a ciência dos materiais.
• Escalabilidade: A abordagem não é limitada pelo tamanho do sistema da mesma forma que a ED, nem por problemas de sinal como o QMC. Os autores sugerem que este framework poderia levar a uma vantagem quântica precoce, particularmente para simular interações retardadas e estender-se para DMFT de múltiplas impurezas (aglomerados), desde que a parametrização do banho seja bem otimizada.

O artigo conclui que, embora limitações relacionadas ao ruído e à necessidade de pós-processamento clássico permaneçam, a metodologia proposta reduz efetivamente tanto os custos clássicos quanto os quânticos, convidando a uma exploração adicional de frameworks híbridos para química quântica e pesquisa em materiais.