A Givens-exchange ansatz for molecular variational eigensolvers

Resumo Técnico: Um ansatz de troca de Givens para resolvedores variacionais moleculares de autovalores

Definição do Problema
O cálculo preciso das energias do estado fundamental molecular é essencial para a descoberta de fármacos computacional, particularmente para sistemas que envolvem correlação forte ou grandes espaços ativos, onde as aproximações clássicas de estrutura eletrônica tornam-se dispendiosas ou pouco confiáveis. O Solucionador de Autovalores Quânticos Variacionais (VQE) aborda isso otimizando um circuito quântico parametrizado para minimizar o valor de expectativa de um hamiltoniano de qubits. No entanto, a eficácia prática do VQE é limitada pelo design do ansatz:

• Ansatzes motivados quimicamente (ex: UCCSD) codificam estruturas de excitação física, mas frequentemente resultam em circuitos profundos demais para o hardware de curto prazo.
• Ansatzes eficientes de hardware utilizam camadas rasas, mas podem carecer da expressividade para capturar correlações moleculares ou sofrer de baixa treinabilidade.
• Métodos de Busca de Arquitetura Quântica (QAS) (ex: abordagens baseadas em Aprendizado por Reforço como HyRLQAS) automatizam a construção de circuitos para encontrar arquiteturas compactas e de alto desempenho. Embora eficazes, esses métodos introduzem uma sobrecarga computacional significativa através da construção repetida de circuitos e atualizações de política, e as arquiteturas resultantes podem ser difíceis de interpretar quimicamente.

O artigo aborda a necessidade de um ansatz que equilibre precisão, reprodutibilidade e interpretabilidade sem depender do custo computacional de uma busca de arquitetura irrestrita.

Metodologia
Os autores propõem um ansatz de troca de Givens de topologia fixa que evita totalmente a busca de arquitetura. A metodologia fundamenta-se na interpretação geométrica da variedade de Stiefel complexa, onde o ansatz representa uma trajetória unitária ordenada que preserva a normalização do estado.

• Estrutura do Circuito:
  • Inicialização: O circuito parte de um estado de base computacional $|b_0\rangle$ selecionado para possuir o menor valor de expectativa diagonal do hamiltoniano.
  • Camadas: O ansatz consiste em $L=2$ blocos de macros de troca de Givens de todos os pares intercalados com camadas locais de rotação $R_Y$.
  • Macro de Troca de Givens ($G_{ij}$): Uma operação de dois qubits atuando nos qubits $i$ e $j$ que rotaciona o subespaço de excitação única abrangido por $|01\rangle$ e $|10\rangle$, enquanto deixa $|00\rangle$ e $|11\rangle$ invariantes. Em uma base de ocupação de Jordan–Wigner (JW), isso corresponde a uma operação de excitação única que conserva partículas.
  • Sequência: O estado completo é preparado como $|\psi(\theta)\rangle = R_L G_{L-1} \dots G_0 R_0 |b_0\rangle$, onde as camadas $R$ contêm rotações $R_Y$ locais e as camadas $G$ contêm sequências ordenadas de trocas de Givens para todos os pares de qubits.
• Protocolo de Otimização:
  • Os parâmetros (ângulos de rotação $\alpha$ e ângulos de troca $\beta$) são otimizados usando o otimizador Adam para minimizar o valor de expectativa do hamiltoniano.
  • Inicialização: Os ângulos são inicializados a partir de uma distribuição normal de média zero (escala 0.05).
  • Avaliação: A diagonalização exata é usada apenas após a terminação da otimização para calcular as energias do estado fundamental, erros e fidelidades. Ela não faz parte do loop variacional.
  • Reprodutibilidade: Os resultados são reportados como médias de seis sementes (sementes: 17, 11111, 22222, 33333, 44444, 55555) para demonstrar o desempenho típico em vez de selecionar apenas a melhor semente.

Principais Contribuições

1. Definição de uma Topologia Fixa: Os autores definem um ansatz de troca de Givens matematicamente explícito com uma sequência ordenada e fixa de camadas de rotação e troca, eliminando a necessidade de busca de arquitetura.
2. Interpretação da Variedade de Stiefel: O trabalho enquadra o ansatz como uma trajetória estruturada na variedade de Stiefel, fornecendo uma justificativa geométrica para o uso de macros de troca ordenadas como blocos de construção variacionais com significado químico.
3. Protocolo de Benchmarking Rigoroso:
   • O estudo distingue entre hamiltonianos de coeficientes verificados (LiH-6, H2O-8) e candidatos de especificação pública (BeH2-6) ou casos de estresse (LiH-4 mapeado por paridade).
   • Separa as comparações primárias de Jordan–Wigner dos casos de mapeamento de paridade/redução para evitar a confusão entre diferentes paisagens variacionais.
   • Reporta resultados como médias de seis sementes com desvios padrão, oferecendo uma visão transparente da sensibilidade da otimização.
4. Clareza na Contabilização de Recursos: O artigo separa explicitamente as contagens de macros pré-compilação (trocas de Givens) das contagens de portas nativas de hardware (CNOTs), reconhecendo que as macros de Givens requerem compilação para portas nativas.

Resultados
O método foi avaliado nos hamiltonianos LiH-6, BeH2-6 e H2O-8 usando o mapeamento Jordan–Wigner sem redução.

• Precisão Química: Em todas as seis sementes, o método alcançou precisão química (erro < 1.60 mHa) para LiH-6, H2O-8 e o candidato BeH2-6 em cada execução.
• Magnitudes de Erro:
  • LiH-6: Erro médio de $0.000000124$ Ha ($0.000124$ mHa).
  • H2O-8: Erro médio de $0.000128558$ Ha ($0.128558$ mHa).
  • BeH2-6: Erro médio de $0.000002152$ Ha ($0.002152$ mHa).
• Comparação com Baselines de QAS:
  • No LiH-6, o erro médio das seis sementes ($0.000124$ mHa) foi inferior aos erros de ponto publicados pelo HyRLQAS ($0.590$ mHa), TensorRL-QAS, CRLQAS e outros baselines baseados em aprendizado por reforço.
  • No H2O-8, o erro médio ($0.128558$ mHa) também foi inferior aos erros de ponto publicados dos métodos QAS comparados.
• Testes de Estresse:
  • O método mostrou sensibilidade à inicialização e mudanças de mapeamento. Para o caso de estresse LiH-4 mapeado por paridade, 5 de 6 sementes passaram na precisão química, mas a distribuição incluiu soluções locais de alto erro.
  • Para o LiH-6 em geometria alongada (3.4 Å), o desempenho degradou significativamente, com apenas 1/6 sementes passando para o mapeamento JW e 0/6 para o mapeamento de paridade.
• Ablação: Uma ablação de profundidade mostrou que $L=2$ (dois blocos de troca) era necessário para alcançar a precisão química em todas as sementes; $L=1$ foi insuficiente para BeH2 e H2O.

Significância e Alegações
O artigo afirma que um ansatz de troca de Givens de topologia fixa serve como um template de referência preciso, reprodutível e livre de busca para VQEs moleculares.

• Precisão vs. Custo de Busca: O método atinge precisão de estado da arte (menores erros médios que vários métodos QAS) sem a sobrecarga computacional de aprender uma política de posicionamento de portas.
• Compromisso (Trade-off): Os autores reconhecem que seu ansatz utiliza um "orçamento de macro pré-compilação" maior (mais operações de dois qubits) do que alguns circuitos compactos otimizados por busca. Eles enquadram isso não como uma vitória de recursos de hardware, mas como um compromisso: pagar um custo maior de macro-operação para ganhar um prior transparente e sistematicamente definido que evita a instabilidade e o custo da busca de arquitetura.
• Escopo: O método é mais eficaz para hamiltonianos de equilíbrio próximo com interpretações diretas de tipo ocupação (ex: Jordan–Wigner). Seu desempenho degrada em ligações alongadas ou hamiltonianos efetivos mapeados por paridade, destacando a importância do alinhamento geométrico do ansatz com a codificação específica.
• Papel: A principal contribuição é fornecer um baseline contra o qual métodos de busca e futuros variantes compilados podem ser avaliados, demonstrando que templates moleculares simples e fixos podem superar arquiteturas aprendidas complexas em termos de precisão média e reprodutibilidade.