Electron-molecule scattering via R-matrix variational algorithms on a quantum computer

Resumo Técnico: Espalhamento Elétron-Molécula via Algoritmos Variacionais R-matrix em Computador Quântico

Declaração do Problema
Colisões elétron-molécula são fundamentais para processos naturais e tecnologias, como o processamento por plasma. Embora o método R-matrix seja uma estratégia computacional padrão para esses problemas, ele enfrenta desafios significativos de escalabilidade quando aplicado a sistemas complexos. O problema da região interna, que envolve resolver o Hamiltoniano de muitos elétrons dentro de um volume confinado para capturar correlação e troca eletrônica, escala de forma deficiente com o número de elétrons e o tamanho do conjunto de base. Abordagens convencionais frequentemente dependem de métodos aproximados ou tornam-se computacionalmente proibitivas para soluções de alta precisão. Além disso, algoritmos variacionais padrão de química quântica geralmente visam o estado fundamental, ao passo que o espalhamento R-matrix requer a determinação simultânea de múltiplos autoestados (incluindo pseudoestados excitados e estados do tipo contínuo) dentro de um setor de simetria específico para definir as condições de contorno do espalhamento.

Metodologia
Os autores propõem um framework computacional quântico que mapeia o problema da região interna R-matrix em um Hamiltoniano de qubits usando o Variational Quantum Eigensolver (VQE) e suas variações. A abordagem envolve vários componentes-chave:

1. Mapeamento do Hamiltoniano e Restrições: O Hamiltoniano de $(N+1)$ elétrons para a região interna é mapeado para qubits via transformação de Jordan–Wigner. Uma restrição física crítica é imposta: no máximo um elétron pode ocupar qualquer orbital do contínuo. Isso é alcançado estruturalmente dentro do ansatz variacional usando uma formulação de operador de projeção de número, garantindo que o estado preparado $|\psi(\theta)\rangle$ satisfaça $P|\psi(\theta)\rangle = |\psi(\theta)\rangle$, onde $P$ projeta sobre configurações de ocupação válidas.
2. Arquitetura do Circuito: O algoritmo utiliza um sistema de três registradores: um registrador alvo ($t$), um registrador de contínuo ($c$) e um registrador seletor ($a$). O circuito prepara estados de prova que abrangem uma base truncada de funções de canal antissimetrizadas e configurações ligadas ($L^2$).
   • Simetria de Spin: O ansatz impõe simetria de spin total ($S, M$) usando rotações fixas de Clebsch–Gordan (CG) (implementadas como rotações Givens de dois qubits) para acoplar o spin do alvo com o spin do elétron incidente.
   • Otimização de Subespaço Sequencial (SSO): Para recuperar o espectro completo de autoestados sem a sobrecarga de diagonalização de subespaço ou termos de penalidade, os autores empregam um protocolo SSO. Isso envolve decompor uma rotação $SO(k)$ em uma cascata de rotações Givens de dois qubits. A otimização procede sequencialmente: o autoestado de menor energia é otimizado primeiro, e seus parâmetros são fixados. Rodas subsequentes otimizam o próximo autoestado dentro do complemento ortogonal, efetivamente diagonalizando em blocos o Hamiltoniano um estado de cada vez. Isso permite o uso de um simples valor esperado do Hamiltoniano $\langle H \rangle$ como função de custo, evitando a necessidade de termos de Hamiltoniano ao quadrado ($H^2$) ou medições de sobreposição.
3. Leitura e Amplitudes de Contorno: Um protocolo de soma coerente é introduzido para estimar valores esperados no registrador do sistema, selecionando pós-factivamente o registrador seletor, reduzindo o número de valores esperados distintos necessários. Crucialmente, os parâmetros ótimos do circuito codificam diretamente os coeficientes de expansão ($a_{ijk}$) da função de onda. Esses coeficientes são usados para calcular as amplitudes de contorno R-matrix ($w_{ik}(a)$), que servem como condições de correspondência para o cálculo de espalhamento da região externa.

Contribuições Principais

• Primeira Formulação Quântica da Região Interna R-matrix: Até onde os autores sabem, esta é a primeira aplicação de algoritmos quânticos ao espalhamento elétron-molécula e a primeira formulação do problema da região interna R-matrix em um computador quântico.
• Circuito de Imposição de Simetria: O trabalho introduz o primeiro uso de um circuito de imposição de simetria de Clebsch–Gordan dentro de um ansatz de espalhamento variacional, utilizando blocos Givens fixos para manter a simetria de spin.
• Otimização Multi-Estado Eficiente: O protocolo SSO proposto difere de métodos existentes como SSVQE e MCVQE ao impor ortogonalidade através da estrutura do circuito em vez de termos de penalidade ou medições de sobreposição. Isso elimina a necessidade de medir elementos de Hamiltoniano fora da diagonal ou realizar diagonalização clássica de subespaço, reduzindo a função de custo a avaliações simples de $\langle H \rangle$.
• Extração Direta de Dados de Espalhamento: O método demonstra que parâmetros variacionais ótimos produzem diretamente as amplitudes de contorno R-matrix, contornando a necessidade de pós-processamento adicional para extrair observáveis de espalhamento.

Resultados
O algoritmo foi demonstrado em um problema modelo de espalhamento de elétrons pela molécula de hidrogênio ($H_2$) usando um simulador clássico sem ruído.

• Configuração do Sistema: A simulação usou uma base mínima (STO-3G) com $N=2$ elétrons alvo, resultando em um espaço de prova de dimensão $k=5$ (quatro ramos de canal aberto e uma configuração ligada) dentro do setor de simetria $S=1/2, M=-1/2, B_{1u}$.
• Desempenho: A abordagem SSO recuperou com sucesso o espectro completo de cinco autoestados com erros de energia dentro de $10^{-7}$ Hartree.
• Eficiência de Recursos: O método SSO exigiu apenas 573 avaliações de função de custo envolvendo 92 termos de Pauli. Em contraste, uma otimização simultânea minimizando a soma das variâncias exigiu 13.345 avaliações de 1.080 termos de Pauli, e uma abordagem de Hamiltoniano dobrado de estado único exigiu 1.397 avaliações.
• Profundidade do Circuito: O circuito não otimizado para este exemplo mínimo utilizou 7 qubits, 217 portas CNOT e uma profundidade de 314.

Significado e Alegações
O artigo alega fornecer uma rota concreta para aplicar hardware quântico ao espalhamento elétron-molécula, abordando especificamente o gargalo da região interna do método R-matrix. Os autores enfatizam que sua abordagem:

1. Evita o armazenamento de grandes matrizes de Hamiltoniano representando estados de prova em qubits.
2. Reduz a sobrecarga de medição evitando observáveis de Hamiltoniano ao quadrado e utilizando otimização sequencial.
3. Fornece uma saída "nativa do problema" onde os parâmetros do circuito mapeiam diretamente para quantidades físicas de espalhamento (amplitudes de contorno).

Os autores notam modestamente que, embora não aleguem uma vantagem quântica assintótica comprovada sobre os melhores métodos iterativos clássicos no momento, o método oferece um paradigma híbrido quântico-clássico viável. Eles destacam que a abordagem é particularmente promissora à medida que a dimensão do espaço de prova aumenta, onde a diagonalização clássica escala como $O(k^3)$, enquanto a abordagem quântica escala o número de parâmetros e a profundidade como $O(k^2)$. O trabalho estabelece as bases para futuras extensões a alvos maiores e representações de contínuo mais sofisticadas, potencialmente integrando processadores quânticos em códigos de espalhamento estabelecidos como o UKRmol+.