Symmetry-adapted qubit encoding with complete active space and Bravyi--Kitaev mapping for quantum chemistry on a quantum computer

Resumo Técnico: Codificação de Qubits Adaptada à Simetria com Espaço Ativo Completo e Mapeamento de Bravyi–Kitaev

Definição do Problema
Simulações de química quântica em processadores quânticos de curto prazo (NISQ) e tolerantes a falhas enfrentam restrições significativas de recursos, impulsionadas principalmente pelo número de qubits necessários para codificar estruturas eletrônicas moleculares e pela complexidade (peso de Pauli e profundidade do circuito) dos Hamiltonianos resultantes. Embora a aproximação de Espaço Ativo Completo (CAS) reduza o tamanho do problema ao congelar orbitais de núcleo e descartar orbitais virtuais, as implementações padrão frequentemente falham em explorar totalmente as simetrias para reduzir ainda mais a contagem de qubits. Além disso, diferentes mapeamentos férmion-qubit, como Jordan–Wigner (JW) e Bravyi–Kitaev (BK), oferecem compromissos entre localidade e requisitos de qubits que nem sempre são integrados de forma otimizada com técnicas de redução de simetria.

Metodologia
Os autores propõem um framework unificado chamado Codificação Adaptada à Simetria com Espaço Ativo Completo (SAE-CAS). Esta abordagem trata a aproximação CAS não meramente como uma truncagem do espaço de Hilbert, mas como a imposição de simetrias $Z$ aproximadas para orbitais de núcleo congelados e orbitais virtuais.

1. Framework Teórico:

   • Simetrias Aproximadas: O método estende as codificações adaptadas à simetria (SAEs) exatas para incluir simetrias aproximadas correspondentes a ocupações orbitais fixas (núcleo congelado e virtual). Estas são modeladas como simetrias $Z$ onde qubits específicos são restringidos a autovalores de $-1$ (congelado) ou $+1$ (virtual).
   • Transformações de Clifford Afins: A redução é alcançada via mapas de Clifford afins. Os autores definem uma matriz de caráter orbital binário $S$ representando as simetrias. Ao construir um mapa afim invertível $T$ sobre $\mathbb{F}_2$, eles transformam a base computacional de modo que as restrições de simetria correspondam a valores fixos (tipicamente 0) em qubits específicos.
   • Remoção de Qubits: Uma vez que o Hamiltoniano é conjugado pelo operador Clifford, termos contendo $X$ ou $Y$ nos qubits restringidos são descartados, e os operadores $Z$ restringidos são substituídos por escalares. Isso remove permanentemente os qubits redundantes, reduzindo o sistema de $n$ orbitais de spin para $n-k$ qubits, onde $k$ é o número de simetrias independentes.
   • Prova de Equivalência: Os autores provam que o Hamiltoniano de qubit de espaço ativo resultante é exatamente equivalente ao Hamiltoniano CAS canônico derivado via projeção de núcleo congelado e de orbital virtual na teoria de estrutura eletrônica molecular. O erro de aproximação surge unicamente da escolha do espaço ativo, não da codificação.

2. Integração com Bravyi–Kitaev (SAE-CAS-BK):

   • O framework é estendido para o mapeamento Bravyi–Kitaev (BK). Como a transformação BK é, em si, uma mudança de base de Clifford afim atuando nas bitstrings computacionais, ela pode ser composta com o mapa Clifford SAE-CAS.
   • A codificação SAE-CAS-BK resultante é unitariamente equivalente ao SAE-CAS. Ela preserva a contagem de qubits, o número de termos de Pauli e os parâmetros variacionais, mas altera a estrutura de localidade (peso de Pauli por termo) e a profundidade de circuito/contagens de CNOT resultante.

3. Benchmarking Numérico:

   • O método foi testado em nove moléculas pequenas (incluindo H$_2$O, C$_2$H$_4$, O$_2$, N$_2$ e C$_4$H$_6$) usando uma base STO-3G mínima.
   • Dois ansátzes variacionais foram empregados: Unitary Coupled Cluster with Single and Double excitations (UCCSD) e um ansátze de hardware-efficient shifted-circular-alternating (HE-SCA).
   • Comparações foram feitas contra JW padrão, JW-CAS (CAS com JW, mas sem redução de simetria) e JW-CAS com filtragem de simetria (SF) ao nível do circuito.

Principais Resultados

• Redução de Recursos: O SAE-CAS consistentemente reduz o número de qubits, a contagem de operadores de Pauli, a profundidade do circuito, a contagem de CNOT e os parâmetros variacionais em comparação ao JW-CAS. Por exemplo, no CAS (4,4) de C$_4$H$_4$, o SAE-CAS reduziu a contagem de qubits de 8 (JW-CAS) para 4 e a contagem de Pauli de 97 para 46.
• Convergência e Precisão:
  • UCCSD: O SAE-CAS atingiu energias convergidas comparáveis ao JW-CAS e JW-CAS (SF), mas exigiu significativamente menos iterações de otimização (ex: 21 vs. 109 para N$_2$). Notavelmente, para C$_4$H$_4$ e C$_4$H$_6$, os ansátzes UCCSD de JW-CAS e JW-CAS (SF) falharam em rodar devido à exaustão de memória causada pelo grande número de excitações, enquanto o SAE-CAS obteceu sucesso ao filtrar a lista de excitações na etapa de mapeamento.
  • HE-SCA: O SAE-CAS demonstrou propriedades de convergência superiores. Para todas as moléculas testadas, o SAE-CAS convergiu para as energias de referência CAS dentro dos orçamentos de camadas testados. Em contraste, o JW-CAS falhou em convergir para O$_2$ e CO dentro de 200 camadas. Os autores atribuem isso ao efeito de "platô estéril" (barren plateau) e vazamento de setor no JW-CAS, onde o manifold variacional abrange setores de paridade de spin ou número de partículas incorretos, enquanto o SAE-CAS restringe a busca ao setor de simetria correto nativamente.
  • SAE-CAS vs. SAE-CAS-BK: As duas codificações produziram energias convergidas e contagens de iteração idênticas, confirmando sua equivalência unitária. Diferenças foram observadas apenas na profundidade do circuito e contagens de CNOT, que variaram em menos de 6% para a maioria das moléculas (exceto CO, que viu um aumento de +20% para BK). Os autores observam que a vantagem de localidade assintótica $O(\log n)$ do BK ainda não foi plenamente realizada nestes pequenos tamanhos de espaço ativo.

Significância e Alegações
O artigo afirma que o SAE-CAS oferece um caminho sistemático e composável para simulações moleculares eficientes em recursos. Ao tratar a aproximação CAS como um conjunto de simetrias aproximadas, o método unifica a redução de simetria exata (grupo de ponto e paridade de spin) com a truncagem de espaço ativo. A principal significância reside em:

1. Imposição Nativa de Simetria: Impor simetrias no nível da codificação (em vez de post-hoc no ansátze) remove direções de busca redundantes e impede que o otimizador derive para setores de número de partículas ou spin incorretos.
2. Escalabilidade: O método permite simulações de espaços ativos maiores que seriam intratáveis devido a restrições de memória na construção do ansátze (como visto nos casos C$_4$H$_4$ e C$_4$H$_6$).
3. Flexibilidade: O framework é compatível com qualquer mapeamento férmion-qubit expressável como uma mudança de base de Clifford afim, demonstrado aqui com ambos JW e BK.

Os autores fornecem uma implementação open-source no pacote Python QuantumSymmetry e concluem que o SAE-CAS é uma rota viável para reduzir a contagem de qubits e a complexidade de operadores tanto em processadores quânticos de curto prazo quanto em tolerantes a falhas.