Point-group symmetry analysis of many-electron wavefunctions on a quantum computer

Resumo Técnico: Análise de simetria de grupo pontual de funções de onda de muitos elétrons em um computador quântico

Declaração do Problema
Os grupos pontuais, que representam operações de simetria espacial em sistemas moleculares, são ferramentas fundamentais na química para a análise de orbitais moleculares e espectroscopia. Embora algoritmos quânticos tenham sido propostos para explorar essas simetrias — como para reduzir a contagem de qubits ou simplificar ansatzes de circuitos —, implementações práticas de operações de simetria de grupo pontual e análises detalhadas de simetria de funções de onda realistas de muitos elétrons permanecem amplamente pouco exploradas. Especificamente, há uma falta de métodos para implementar concretamente operações de simetria e analisar as propriedades de simetria (por exemplo, pesos de representações irredutíveis) de estados realistas e correlacionados de muitos elétrons em hardware quântico de curto prazo.

Metodologia
Os autores propõem um método híbrido quântico-clássico, sem ancilla, para analisar simetrias de grupo pontual de estados de muitos elétrons, aplicável tanto a grupos abelianos quanto não abelianos. O núcleo do método envolve o cálculo do "peso" ($w_\Gamma$) de representações irredutíveis ($\Gamma$) dentro de uma função de onda dada $|\Psi\rangle$. Este peso é definido como o coeficiente ao quadrado da projeção de $|\Psi\rangle$ sobre o subespaço de $\Gamma$.

O quadro teórico baseia-se na fórmula do operador de projeção:
$$w_\Gamma = \frac{d_\Gamma}{|G|} \sum_C \chi^*_\Gamma(C) \sum_{g \in C} \langle\Psi|\hat{g}|\Psi\rangle$$
onde $d_\Gamma$ é a dimensão da representação, $|G|$ é a ordem do grupo, $\chi_\Gamma(C)$ é o caráter da classe de conjugação $C$ e $\hat{g}$ é a operação de simetria.

Para avaliar $\langle\Psi|\hat{g}|\Psi\rangle$ em um computador quântico sem qubits ancilla, os autores utilizam a diagonalização de matrizes de representação. Para orbitais adaptados à simetria (por exemplo, Hartree-Fock), a operação de simetria $U(g)$ é decomposta em rotações de orbitais. Ao diagonalizar a matriz de representação $\bar{D}_\gamma(g)$ para conjuntos de orbitais degenerados, o unitário $U(g)$ é reescrito como um produto de uma rotação de base $\tilde{U}(g)$ e um operador de fase diagonal. O valor esperado é então obtido preparando o estado $|\tilde{\Psi}(g)\rangle = \tilde{U}(g)|\Psi\rangle$, medindo na base computacional e combinando os resultados com fatores de fase apropriados. Esta abordagem evita os circuitos profundos e os requisitos de ancilla de testes de Hadamard ou técnicas de Combinação Linear de Unitários (LCU).

O estudo emprega duas estratégias principais de validação:

1. Simulações Numéricas Clássicas: Utilizando as bibliotecas ffsim e pyscf, o método é aplicado ao benzeno ($D_{6h}$) e ao ferroceno escalonado ($D_{5d}$). Os autores analisam determinantes de Slater únicos (estados fundamentais e estados de excitação simples de Hartree-Fock) e funções de onda correlacionadas (funções Jastrow de Cluster Unitário e Jastrow de Cluster Unitário Local).
2. Demonstração em Hardware: O método é executado no dispositivo quântico supercondutor ibm_kawasaki da IBM (até 32 qubits). O sistema alvo é o benzeno na simetria $D_{2h}$ (um subgrupo de $D_{6h}$). Funções de onda de teste (estados fundamentais e primeiros estados excitados) são preparadas usando cálculos de Grupo de Renormalização de Matriz de Densidade (DMRG), comprimidas em um ansatz de parede de tijolos via técnicas de rede de tensores e mapeadas para circuitos quânticos. Para abordar o ruído, duas técnicas de mitigação de erro são comparadas: Extrapolção de Ruído Zero (ZNE) via dobragem de portas e uma abordagem de mitigação de erro quasi-probabilística (QESEM).

Principais Resultados

• Simulações Clássicas:
  • Para estados fundamentais de Hartree-Fock de benzeno e ferroceno, o método identifica corretamente o peso da representação irredutível totalmente simétrica ($A_{1g}$) como próximo de 1,0, consistente com configurações de camada fechada.
  • Para estados de excitação simples (SE), o método decompõe com sucesso as funções de onda em múltiplas representações irredutíveis. Por exemplo, um estado SE específico no benzeno foi decomposto em $\approx 0,005 B_{1u} + 0,495 B_{2u} + 0,5 E_{1u}$. A soma dos pesos sobre todas as configurações SE independentes produziu resultados invariantes consistentes com as fórmulas de redução da teoria de grupos.
  • Na análise de funções de onda LUCJ correlacionadas, o método revelou que procedimentos de otimização que minimizam o erro de energia (sem regularização) podiam degradar significativamente o peso de simetria (desviando-se do peso correto de $A_{1g}$), enquanto a regularização ajudava a preservar as propriedades de simetria, embora com diferentes compensações de energia. Isso sugere que a métrica é útil para avaliar a qualidade de funções de onda de teste.
• Demonstração em Hardware:
  • No dispositivo ibm_kawasaki, medições brutas mostraram uma degradação significativa dos pesos de simetria à medida que a contagem de qubits aumentava (por exemplo, o peso $A_g$ para o estado fundamental caiu para ~0,54 para 32 qubits sem mitigação).
  • A aplicação de mitigação de erro foi crucial. Embora a ZNE tenha melhorado os resultados, a abordagem quasi-probabilística QESEM produziu a maior fidelidade, reproduzindo os pesos corretos de representação irredutível com um erro de poucos por cento (por exemplo, $0,968 \pm 0,010$ para o estado fundamental em 32 qubits, comparado a um valor exato de $0,999$).
  • O experimento confirmou que o método proposto sem ancilla, combinado com preparação de estado de rede de tensores e mitigação de erro avançada, pode reproduzir fielmente as propriedades de simetria no hardware atual.

Significado e Alegações
O artigo afirma preencher a lacuna entre a projeção de simetria teórica e a implementação prática em hardware quântico. O método proposto serve como uma ferramenta prática para:

1. Análise de Simetria: Analisar as propriedades de simetria de funções de onda de muitos elétrons em simulações realistas de materiais em computadores quânticos de curto prazo e de falha tolerante inicial.
2. Avaliação de Função de Onda: Avaliar a qualidade de funções de onda de teste (por exemplo, em algoritmos variacionais) quantificando sua adesão a simetrias específicas.
3. Benchmarking: Fornecer um benchmark para dispositivos quânticos reais e técnicas de mitigação de erro, uma vez que o método não requer qubits ancilla e possui sobrecarga mínima para orbitais adaptados à simetria.

Os autores enfatizam que, embora a demonstração em hardware tenha utilizado um subgrupo abeliano ($D_{2h}$) e avaliação baseada em Pauli, a metodologia subjacente é geral e aplicável a grupos não abelianos e funções de base arbitrárias. O trabalho sugere que a análise de simetria é um componente viável e necessário para futuras simulações quânticas de materiais complexos.