Probing the Ground State of the Antiferromagnetic Heisenberg Model on the Kagome Lattice using Geometrically Informed Variational Quantum Eigensolver

Resumo Técnico: Sondando o Estado Fundamental do Modelo de Heisenberg Antiferromagnético em uma Rede de Kagome usando VQE Geometricamente Informado

Enunciado do Problema
O artigo aborda o desafio de encontrar o estado fundamental do modelo de Heisenberg antiferromagnético em uma rede de kagome, um sistema conhecido por sua frustração geométrica e potenciais fases de líquido de spin quântico. Determinar o estado fundamental de tais Hamiltonianos $k$-locais é um problema QMA-hard, tornando a simulação clássica intratável para sistemas grandes devido ao crescimento exponencial do espaço de Hilbert. Embora os Algoritmos Variacionais de Autovalores Quânticos (VQE) ofereçam uma abordagem híbrida quântico-clássica adequada para dispositivos NISQ (Noisy Intermediate-Scale Quantum), eles enfrentam obstáculos significativos: a necessidade de circuitos profundos, suscetibilidade ao ruído (erros coerentes e incoerentes), platôs de barren (barren plateaus) onde os gradientes desaparecem e a dificuldade de projetar ansatzes que sejam simultaneamente expressivos e eficientes em termos de hardware. Especificamente, os autores visam caracterizar as propriedades do estado fundamental, como correlações spin-spin e fatores de estrutura estática, em hardware quântico real sem depender de reduções pesadas específicas do problema ou de mitigação extensiva de erros durante o loop de otimização.

Metodologia
Os autores propõem uma abordagem VQE geometricamente informada adaptada para a rede de kagome, focando em duas geometrias fundamentais: um triângulo único e uma estrela de kagome (seis triângulos que compartilham cantos).

1. Ansatz Geometricamente Informado:
   A inovação central é um ansatz customizado, raso e eficiente em termos de hardware. Diferente do Ansatz de Hamiltoniano Variacional (HVA), que pode ser muito profundo para dispositivos NISQ, este ansatz é construído usando um padrão de rotações de qubit único parametrizadas ($R_y$) e portas de emaranhamento restritas à conectividade nativa do hardware (evitando portas SWAP).

   • Espaço de Parâmetros Euclidiano: Uma característica crítica deste design é que o tensor métrico de Fubini-Study é provadamente diagonal e constante. Isso garante que o espaço de parâmetros seja naturalmente Euclidiano, livre de singularidades e de sobreparametrização. Consequentemente, o gradiente natural quântico é equivalente ao gradiente padrão, simplificando o cenário de otimização.
   • Compatibilidade com Hardware: A profundidade do circuito é minimizada (uma camada de emaranhamento) e o número de parâmetros é reduzido através da alternância na colocação das portas $R_y$.

2. Estratégia de Otimização (AQNGD/AGD):
   Os autores utilizam um otimizador de Descida de Gradiente Adaptativa (AGD). Como o ansatz produz um métrico de Fubini-Study diagonal e constante, o Gradiente Natural Quântico Adaptativo (AQNGD) reduz-se a AGD. A taxa de aprendizado é adaptada dinamicamente usando a regra de backtracking de Armijo. Isso permite passos iniciais grandes para acelerar a convergência, enquanto reduz automaticamente o tamanho do passo conforme o otimizador se aproxima do mínimo, evitando a estagnação em geometrias estreitas. Isso é comparado à Aproximação Estocástica de Perturbação Simultânea (SPSA).

3. Mitigação de Erro (Pós-processamento):
   A mitigação de erro é aplicada após a otimização do VQE para refinar as estimativas finais de energia e propriedades estruturais, em vez de durante o treinamento iterativo.

   • Mitigação de Erro de Leitura (REM): É utilizada uma abordagem de matriz de resposta particionada para reverter erros de SPAM (State Preparation and Measurement), escalando de forma mais eficiente do que a inversão de matriz completa para medições esparsas.
   • Extrapolação de Ruído Zero (ZNE): Os autores aplicam ZNE usando extrapolação quadrática (e regressão polinomial Bayesiana) para extrapolar os resultados para o limite de ruído zero. Eles observam que, embora o ZNE não obedeça estritamente ao princípio variacional de Rayleigh-Ritz, ele é utilizado aqui para melhorar a precisão.

Resultos Principais
O estudo foi conduzido em hardware quântico real (IBM Oslo, IBM Kyoto, IBM Torino e IBM Algiers) e validado via simulações de vetor de estado.

• Convergência e Otimização:

  • Em um sistema de triângulo de 3 qubits, o VQE otimizado por AGD convergiu suavemente para a energia do estado fundamental ($E \approx -3$) com alta fidelidade (99,89% na simulação) usando apenas 19 iterações. Os parâmetros finais corresponderam de perto aos valores teóricos exatos.
  • Em uma estrela de kagome de 12 qubits, o AGD demonstrou convergência superior ao SPSA, exigindo menos iterações para atingir a energia do estado fundamental ($E \approx -18$). A fidelidade final na simulação foi de 96,37%.
  • A adaptação do tamanho do passo via regra de Armijo mostrou-se eficaz para navegar no cenário de custo, iniciando com passos grandes e refinando conforme a convergência se aproximava.

• Resiliência ao Ruído das Propriedades Estruturais:

  • Uma descoberta significativa é que as correlações spin-spin e o fator de estrutura estática ($S(q)$) são resilientes ao ruído. Mesmo sem mitigação de erro, os mapas experimentais de $S(q)$ retiveram a assinatura qualitativa do estado fundamental (por exemplo, a supressão no momento zero indicando ausência de ordem ferromagnética).
  • O coeficiente de correlação de Pearson entre o $S(q)$ experimental ruidoso e o $S(q)$ teórico exato permaneceu extremamente alto (>99,9%).
  • Experimentos de simulação revelaram uma dependência "em degraus" do fator de estrutura estática em relação aos níveis de ruído: a estrutura de correlação permanece robusta até um limiar crítico de ruído, após o qual degrada rapidamente, assemelhando-se a uma transição de fase. Isso contrasta com a energia, que decai monotonicamente com o ruído.

• Eficácia da Mitigação de Erro:

  • A aplicação de REM e ZNE pós-otimização melhorou significamente as estimativas de energia. Para a estrela de kagome, a energia bruta foi de $\approx -14,9$, enquanto o REM sozinho a melhorou para $\approx -17,3$, e o ZNE (quadrático) a aproximou do valor exato de $-18$.
  • Os autores observaram que combinar REM e ZNE às vezes levava ao excesso (overshooting) da energia do estado fundamental (violando o limite variacional), um comportamento consistente com a natureza não-variacional do ZNE. O ZNE sozinho frequentemente forneceu a melhor precisão entre as configurações testadas.

Significância e Alegações
O artigo afirma que um design de ansatz geometricamente informado, que garante um espaço de parâmetros naturalmente Euclidiano, permite o treinamento eficiente de VQE em hardware NISQ atual sem a necessidade de reduções complexas e específicas do problema.

• Robustez: O trabalho demonstra que propriedades físicas significativas, especificamente correlações spin-spin e o fator de estrutura estática, podem ser extraídas de hardware quântico ruidoso com precisão suficiente para caracterizar o estado fundamental, mesmo sem mitigação de erro durante o loop de otimização.
• Eficiência de Otimização: O uso de AGD com tamanhos de passo adaptativos provou ser mais eficiente que o SPSA para este ansatz específico de eficiência de hardware, alcançando uma convergência mais rápida em termos de iterações.
• Escalabilidade: Ao evitar circuitos profundos e mitigação de erro pesada durante o treinamento, a abordagem oferece um caminho escalável para sondar ímãs quânticos frustrados. Os autores sugerem que este framework pode ser estendido para redes de kagome maiores e outros sistemas quânticos frustrados.
• Modéstia: Os autores reconhecem que, embora seu método tenha caracterizado com sucesso pequenas células fundamentais (triângulo e estrela), o estado fundamental geral da rede de kagome infinita permanece uma questão aberta na literatura. Eles não afirmam ter resolvido o problema da rede infinita, mas demonstraram uma metodologia viável para sondar esses sistemas em dispositivos atuais. Eles também observam que, embora o ZNE melhore as estimativas de energia, ele não preserva o princípio variacional, e que a combinação de REM e ZNE requer interpretação cuidadosa.