Periodic Symmetry-Adapted Encoding: Qubit Reduction in Crystalline Electronic Structure

Resumo Técnico: Codificação Periódica Adaptada à Simetria para Estrutura Eletrônica Cristalina

Definição do Problema
Simular a estrutura eletrônica de materiais cristalinos em computadores quânticos apresenta desafios distintos em comparação com sistemas moleculares. Enquanto simulações moleculares frequentemente utilizam simetrias de grupo de ponto para reduzir a contagem de qubits, sistemas periódicos exigem o tratamento de bases de Bloch em múltiplos pontos k, espaços ativos selecionados a partir de bandas de energia em vez de orbitais isolados, e simetrias de grupo de espaço que incluem translações cristalinas. Propostas existentes para simulação quântica periódica concentraram-se amplamente em bases de ondas planas, modelos de Hubbard ou codificações diretas de pontos k, deixando uma lacuna para codificações periódicas centradas em átomos e quimicamente precisas que explorem plenamente a simetria de grupo de espaço para minimizar os requisitos de qubits. O framework de Codificação Adaptada à Simetria (SAE), anteriormente bem-sucedido para moléculas, não havia sido estendido para sistemas periódicos para identificar e remover sistematicamente qubits redundantes decorrentes de simetrias de translação.

Metodologia
Os autores estendem o framework SAE para a estrutura eletrônica periódica construindo um Hamiltoniano de supercélula de ponto $\Gamma$ a partir de um cálculo de pontos k dobrado. A metodologia procede através das seguintes etapas:

1. Construção do Hamiltoniano: Um cálculo de Hartree–Fock com restrição de ponto k (KRHF) é realizado em uma célula primitiva. Este é dobrado em uma representação de supercélula comensurável onde os orbitais moleculares (MOs) resultantes são de valores reais no ponto $\Gamma$. Um espaço ativo é selecionado a partir desses MOs da supercélula (tipicamente preservando os manifolds quase degenerados HOMO/LUMO), e o Hamiltoniano efetivo do espaço ativo é construído usando integrais de ajuste de densidade (density-fitted integrals).
2. Identificação de Simetria: O framework identifica todos os geradores de simetria $Z_2$ aplicáveis do Hamiltoniano de segunda quantização antes do mapeamento para qubits. Esses geradores dividem-se em três classes:
   • Paridade de Spin: Operadores que contam a paridade de elétrons de spin para cima e para baixo (sempre presentes).
   • Operações de Grupo de Ponto: Simetrias espaciais (rotações, reflexões, inversões) da supercélula que atuam como simetrias $Z_2$ nos MOs ativos.
   • Simetrias de Translação Cristalina: Únicas ao cenário periódico, estas são translações de "meia-supercélula" ao longo de eixos de malha dobrada pares. Uma translação por $(N_i/2)a_i$ mapeia a supercélula para si mesma módulo um vetor de rede, atuando como uma simetria Booleana exata no Hamiltoniano do ponto $\Gamma$ dobrado.
3. Redução de Qubits: Os geradores identificados formam um sistema de equações lineares Booleanas ($A\mathbf{a} = \mathbf{c}$) referentes às ocupações de spin-orbital. Uma transformação de Clifford afim é aplicada para projetar o Hamiltoniano no setor de simetria física alvo. Isso remove um qubit para cada gerador independente, reduzindo o tamanho do registro de $N_{so}$ para $N_{so} - n_g$, onde $n_g$ é o número de geradores independentes.
4. Implementação: O método é implementado no pacote de código aberto QuantumSymmetry, exigindo apenas entradas de química periódica padrão (geometria, conjunto de base, malha de pontos k, índices do espaço ativo) sem a necessidade de especificação manual de geradores de simetria.

Contribuições Principais e Resultados
O artigo testa o método SAE periódico em dez sistemas cristalinos (incluindo diamante, silício, MgO, NaCl, CsCl, h-BN, AlN, $\alpha$-quartzo e MgF2) representando grupos de espaço cúbicos, hexagonais, trigonais e tetragonais.

• Redução de Qubits: Em todo o conjunto, o método remove 4–8 qubits. A redução mais significativa é observada na estrutura B2 CsCl, onde o espaço ativo CAS(6,7) é reduzido de 14 para 6 qubits. Este sistema realiza um grupo de simetria Booleana isomorfo a $Z_2^8$ (duas paridades de spin, três meia-translações e três reflexões de grupo de ponto), excedendo o máximo de $Z_2^5$ típico do SAE molecular.
• Compressão de Circuito: A redução propaga-se para o ansatz variacional. O número de parâmetros variacionais UCCSD é reduzido em 3–8$\times$, e as contagens de CNOT são reduzidas em até 309$\times$ (por exemplo, no CsCl, de 22.240 para 72 CNOTs).
• Precisão: Benchmarks de UCCSD-VQE sem ruído contra a diagonalização exata no setor do espaço ativo mostram que as codificações reduzidas preservam as energias alvo com erro bem abaixo da precisão química (erros $< 10^{-6}$ Ha). As codificações reduzidas convergem em menos avaliações do objetivo VQE (2.8–3.9$\times$ menos) em comparação com as codificações Jordan–Wigner (JW) padrão.
• Validação: A diagonalização exata dos Hamiltonianos reduzidos confirma que os menores autovalores coincidem com as energias FCI do setor não reduzido dentro de $10^{-12}$ Ha, verificando que a projeção afim remove apenas coordenadas redundantes sem alterar o espectro físico.

Significância e Alegações
O artigo alega que o SAE periódico converte a informação existente dentro do problema de estado sólido — especificamente translações cristalinas e involuções de grupo de espaço — diretamente em compressão de qubits e de ansatz.

• Além dos Limites Moleculares: Os autores enfatizam que a inclusão de geradores de translação permite que o SAE periódico exceda os limites de geradores do SAE molecular. Enquanto sistemas moleculares são limitados a no máximo cinco geradores Booleanos independentes (duas paridades de spin e três de grupo de ponto), supercélulas periódicas com eixos de malha dobrada pares podem adicionar até três geradores de meia-translação independentes, permitindo um máximo de oito.
• Pré-processamento Exato: O método é apresentado como um passo de pré-processamento exato que não introduz aproximações físicas ao problema do espaço ativo. Dentro do setor de simetria selecionado, o Hamiltoniano reduzido é isoespectral ao bloco correspondente do Hamiltoniano não reduzido.
• Regra de Design: Uma regra de design específica é identificada: eixos de malha dobrada com tamanhos pares expõem simetrias $Z_2$ adicionais (meia-translações), enquanto eixos ímpares não o fazem.
• Aplicabilidade Geral: Os resultados demonstram que as economias significativas não se restringem a cristais cúbicos de alta simetria; materiais não cúbicos (hexagonais, trigonais, tetragonais) também produzem reduções de 4–5 qubits, desde que o espaço ativo respeite as degenerâncias e caracteres de simetria relevantes.

O trabalho conclui que o SAE periódico efetivamente reduz o custo quântico de simular a estrutura eletrônica cristalina ao eliminar graus de liberdade conhecidos por estarem fora do setor alvo, permitindo assim simulações quânticas mais eficientes em hardware de curto prazo (near-term).