Channel-agnostic finite-temperature phase estimation averaged over variable grids: reconstruction of Green's function for dynamical mean-field theory

Resumo Técnico: Estimativa de fase em temperatura finita agnóstica a canais, média sobre grades variáveis

Declaração do Problema
O artigo aborda o desafio de realizar cálculos de Teoria do Funcional da Densidade Dinâmica (DMFT) para sistemas eletrônicos correlacionados em computadores quânticos, especificamente em temperaturas finitas. Embora a Teoria do Funcional da Densidade (DFT) seja bem-sucedida para sistemas fracamente correlacionados, ela falha para materiais fortemente correlacionados, tornando necessário o enfoque DFT+DMFT. Solucionadores DMFT tradicionais em computadores clássicos frequentemente dependem de Monte Carlo Quântico (QMC) ou Diagonalização Exata (ED/FCI). No entanto, a implementação desses métodos em hardware quântico é difícil porque:

1. Complexidade em Temperatura Finita: Em temperaturas finitas, a Função de Green (GF) envolve uma mistura térmica de muitos autoestados de energia iniciais (canais de excitação), ao contrário do caso de temperatura zero, onde apenas o estado fundamental importa.
2. Canais de Excitação Desconhecidos: Em uma medição padrão de Estimativa de Fase Quântica (QPE) sobre um estado térmico, o autoestado inicial específico $|\Psi_{\lambda_0}\rangle$ e o canal de excitação resultante são desconhecidos para cada disparo de medição. Métodos anteriores baseados em QPE para GFs exigiam conhecimento da energia do estado fundamental ou de canais de excitação específicos, o que não é viável em uma mistura térmica.
3. Vazamento Espectral e Viés de Grade: A QPE padrão sofre de vazamento espectral quando os autovalores não se alinham perfeitamente com a grade discreta dos qubits ancila, introduzindo vieses dependentes das configurações específicas da grade (origem e espaçamento).

Metodologia
Os autores propõem um esquema híbrido quântico-clássico denominado QAVG-DMFT (QPE Média sobre Grades Variáveis). O fluxo de trabalho integra um loop clássico DFT+DMFT com uma sub-rotina quântica para resolver o problema de impureza.

• Parte Quântica (QPE Agnóstica a Canais):

  • Preparação de Estado: O algoritmo começa preparando o estado de Gibbs $\rho_{Gibbs}$ do Modelo de Impureza de Anderson (AIM) isolado no computador quântico.
  • Circuito Modificado ($C_{GF-QPE}$): Em vez da QPE padrão, os autores utilizam um circuito modificado equipado com "circuitos parciais de excitação" ($C_{exc}$). Esses circuitos (diagonais $C_m$ e fora da diagonal $C_{mm'}$) criam superposições de excitações de elétrons e buracos.
  • Agnosticismo de Canal: O circuito aplica portas de Evolução em Tempo Real (RTE) controladas antes e após a excitação. Crucialmente, essa configuração permite a extração de amplitudes espectrais e energias de excitação sem necessidade de conhecer o autoestado inicial específico ou o canal de excitação invocado em cada medição. Os resultados da medição produzem histogramas que aproximam a GF diretamente.
  • Grades Variáveis: Para mitigar o vazamento espectral, o esquema emprega $n_{setting}$ circuitos QPE diferentes, cada um com parâmetros de grade distintos (fator de escala $t_0$ e origem $a_{orig}$).

• Parte Clássica (Pós-processamento QAVG):

  • Estrutura de Otimização: O processador clássico coleta histogramas das várias configurações de QPE. Em seguida, ele reconstrói a GF otimizando um conjunto de "parâmetros de teste" $\Lambda$.
  • Parâmetros de Teste: A GF é modelada usando uma representação fictícia na base de Orbitais Naturais (NO). Os parâmetros incluem:
    • Energias de excitação fictícias ($\varepsilon_{\xi\ell}$).
    • Larguras espectrais fictícias ($\Delta E_{\xi\ell}$) usando uma Densidade de Estados (DOS) quadrática para modelar o alargamento de picos.
    • Amplitudes de transição fictícias ($v^{(\nu)}_{\xi\ell}$), parametrizadas via transformações de Householder hiperesféricas para garantir ortonormalidade.
  • Função de Custo: Uma função de custo $F(\Lambda)$ é definida como a distância média $L_1$ não uniformemente ponderada entre as distribuições de probabilidade modeladas (derivadas dos parâmetros de teste) e os histogramas experimentais em todas as configurações de grade. Os pesos priorizam excitações de baixa energia.
  • Otimização: Os parâmetros de teste são otimizados (por exemplo, via métodos de Monte Carlo) para minimizar a discrepância, efetivamente "removendo ruído" dos dados e reconstruindo a GF.

• Loop DMFT: A GF reconstruída é usada para calcular a autoenergia $\Sigma_{loc}$ via equação de Dyson, que atualiza a função de Green da rede para a próxima iteração do DMFT.

Principais Contribuições

1. QPE em Temperatura Finita Agnóstica a Canais: O artigo introduz um design de circuito que extrai a Função de Green em temperatura finita sem conhecimento prévio dos canais de excitação ou autoestados iniciais, superando uma limitação majoritária dos métodos anteriores de GF baseados em QPE.
2. Metodologia QAVG: A extensão da abordagem QAVG (anteriormente para temperatura zero) para sistemas em temperatura finita. Este método usa média sobre grades variáveis e otimização de parâmetros para reconstruir a GF a partir de dados QPE discretizados e ruidosos, reduzindo vieses associados a escolhas específicas de grade.
3. Parametrização por Orbitais Naturais: O uso de Orbitais Naturais e parametrização hiperesférica para amplitudes de transição permite uma representação ortonormal e fisicamente plausível dos canais de excitação fictícios, reduzindo o número de parâmetros necessários em comparação com o número real de canais de excitação.
4. Demonstração em SrVO$_3$: O esquema é aplicado ao Vanadato de Estrôncio (SrVO$_3$), um metal correlacionado prototípico, demonstrando a reconstrução da Função de Green e da Densidade de Estados (DOS) resolvida em momento.

Resultados
Os autores realizaram simulações numéricas usando dados exatos de FCI como proxy para a saída do computador quântico (para isolar a validade do algoritmo de reconstrução do ruído de hardware).

• Precisão de Reconstrução: Em um teste "de um único disparo" (reconstruindo a GF na 10ª iteração de um loop DMFT), o método QAVG capturou com sucesso as características gerais da DOS-FCI usando apenas um pequeno número de canais de excitação fictícios (6 para elétrons, 2 para buracos) em comparação com os milhares de canais reais que contribuem em temperatura finita.
• DMFT Iterativo: Em uma simulação iterativa QAVG-DMFT, o método reproduziu razoavelmente bem a DOS resolvida em momento do SrVO$_3$. Embora discrepâncias menores tenham sido observadas no regime de baixa energia (especificamente vales próximos a $\pm 0,2$ eV), as características espectrais gerais foram preservadas.
• Importância do Ponderamento: O estudo confirmou que o uso de pesos não uniformes (enfatizando excitações de baixa energia) na função de custo produziu uma reprodução espectral melhor do que uma distância $L_1$ uniforme.

Significado e Alegações
O artigo afirma demonstrar um caminho válido para realizar cálculos DMFT em futuros computadores quânticos tolerantes a falhas.

• Eficiência: Ao evitar a necessidade de diagonalização numérica do subespaço correlacionado no computador quântico e contornar a exigência de conhecer canais de excitação específicos, o esquema torna o processo de amostragem eficiente uma vez que o estado de Gibbs é preparado.
• Robustez: A abordagem QAVG alivia vieses decorrentes da resolução finita das grades QPE, ao fazer média sobre configurações variáveis e otimizar parâmetros de teste.
• Independência: Os autores notam que a abordagem QAVG é independente de metodologias específicas de estrutura eletrônica, o que significa que pode ser aplicada para obter a GF de qualquer sistema correlacionado, não apenas daqueles dentro do arcabouço DFT+DMFT.
• Modéstia: Os autores afirmam explicitamente que os resultados atuais são baseados em simulações numéricas usando distribuições de probabilidade exatas (resultados FCI) em vez de dados de hardware quântico real. Eles reconhecem que erros estatísticos e ruído de hardware ainda não são abordados neste estudo específico, mas planejam relatar aplicações em computadores quânticos reais em trabalhos futuros. O artigo não alega ter resolvido o problema completo do DMFT em temperatura finita em dispositivos ruidosos atuais, mas sim propõe um esquema adequado para a era da computação quântica tolerante a falhas.