Quantum Simulation of Nuclear Shell Model Using GCM-Based Methods on NISQ Devices

Resumo Técnico: Simulação Quântica do Modelo de Camadas Nucleares Usando Métodos Baseados em GCM em Dispositivos NISQ

Definição do Problema
A simulação de sistemas quânticos correlacionados, particularmente sistemas de muitos corpos nucleares, apresenta desafios significativos para computadores clássicos devido ao crescimento exponencial do espaço de Hilbert com o número de coordenadas coletivas. Métodos tradicionais como o Método das Coordenadas Geradoras (GCM) exigem a construção de grandes grades de estados geradores e a resolução de problemas de autovalores generalizados, que se tornam computacionalmente intratáveis à medida que o tamanho do sistema aumenta. Embora o Algoritmo de Eigensolver Variacional Quântico (VQE) e suas variantes adaptativas (ADAPT-VQE) tenham sido propostos para dispositivos de Escala Intermediária com Ruído Quântico (NISQ), eles sofrem com dificuldades de otimização, incluindo platôs estéreis (barren plateaus), convergência para mínimos locais e a necessidade de circuitos quânticos profundos. Além disso, mapeamentos padrão de férmion-para-qubit, como a transformação de Jordan-Wigner (JW) com codificação One-Hot (OH), frequentemente resultam em altas contagens de qubits e circuitos profundos devido a strings de paridade não locais, limitando sua aplicabilidade em hardware atual.

Metodologia
Os autores propõem um framework híbrido quântico-clássico utilizando dois algoritmos principais: o Método de Coordenadas Geradoras Quânticas (QuGCM) e sua variante adaptativa, o Método Inspirado em Coordenadas Geradoras Adaptativas (ADAPT-GCIM).

1. Estrutura Teórica:

   • QuGCM: Esta abordagem constrói a função de onda total como uma superposição de estados geradores não ortogonais. Esses estados são gerados pela aplicação de operadores de Coupled Cluster Unitário (UCC) adaptados por simetria (especificamente UCCSD) a estados de referência de Hartree-Fock (HF). Diferente do VQE, o QuGCM não otimiza parâmetros variacionais; em vez disso, ele prepara esses estados como circuitos quânticos, mede sobreposas (overlaps) de pares e núcleos de Hamiltoniana, e resolve a equação de autovalor generalizada de Hill-Wheeler ($H\psi = EN\psi$) em um computador clássico.
   • ADAPT-GCIM: Para aumentar a eficiência, este método seleciona iterativamente excitações geradoras baseadas em gradientes de energia. Ele expande dinamicamente o subespaço de Hilbert adicionando o operador que mais contribui para a redução da energia, evitando a necessidade de explorar o conjunto completo de base. Isso reduz o número de circuitos quânticos e estados de base necessários.

2. Mapeamentos de Férmion-para-Qubit:
   O artigo compara duas estratégias de codificação:

   • Codificação Jordan-Wigner (JW) com One-Hot (OH): A abordagem padrão onde cada orbital mapeia para um qubit.
   • Codificação Gray Code (GC): Uma ordenação refletida binária onde estados de base consecutivos diferem por um único bit. Os autores demonstram que a codificação GC permite uma representação mais compacta do espaço de Hilbert e gera operadores de Pauli mais simples e locais para excitações, reduzindo assim a profundidade do circuito e a contagem de portas em comparação ao JW.

3. Implementação:
   Os métodos foram implementados para três sistemas nucleares: o deutério (usando o potencial Reid68), $^{38}$Ar (Hamiltoniana de modelo de camadas) e $^{6}$Li (Hamiltoniana de modelo de camadas). As simulações foram realizadas utilizando tanto simuladores de estado vetorial sem ruído quanto simuladores com ruído (FakeBrisbane) para mimetizar o hardware quântico IBM.

Resultos Principais

• Precisão e Convergência: Tanto o QuGCM quanto o ADAPT-GCIM produziram espectros de energia em concordância com os resultados exatos de diagonalização clássica para todos os sistemas testados. O ADAPT-GCIM demonstrou convergência rápida, frequentemente atingindo a energia do estado fundamental em poucas iterações (por exemplo, três iterações para o deutério).
• Resiliência ao Ruído: Sob condições de ruído simulado, os métodos baseados em GCM exibiram maior estabilidade do que os métodos variacionais (VQE/ADAPT-VQE). Especificamente, os resultados do ADAPT-GCIM e do QuGCM permaneceram mais próximos dos valores exatos com menor incerteza estatística (Desvio Absoluto da Mediana) em comparação com seus contrapartes variacionais.
• Eficiência da Codificação GC: A codificação Gray Code superou significativamente a codificação JW. Para o $^{38}$Ar, enquanto o mapeamento JW requer 6 qubits, a codificação GC reduziu o número de qubits necessários para estados específicos (por exemplo, 3 qubits para o estado $0^+$, 2 para o $1^+$, e 1 para o $2^+$) e reduziu drasticamente o número de termos de Pauli e circuitos quânticos necessários para o QuGCM (aproximadamente 27 vezes menor que o JW). Para o $^{6}$Li, a codificação GC reduziu o requisito de qubits de 12 (em JW completo de m-esquema) para 4, permitindo simulações que seriam inviáveis de outra forma.
• Redução de Recursos: A natureza adaptativa do ADAPT-GCIM reduziu drasticamente o número de estados de base necessários. Por exemplo, no $^{38}$Ar, o ADAPT-GCIM exigiu apenas 3–6 estados de base comparado aos 15 totais do QuGCM padrão, diminuindo significativamente a carga computacional.

Significância e Alegações
O artigo afirma que a integração do QuGCM e ADAPT-GCIM com a codificação Gray Code fornece um caminho prático e escalável para simular sistemas nucleares correlacionados em dispositivos NISQ. A principal significância reside em:

• Evitar Armadilhas de Otimização: Ao resolver um problema de autovalor generalizado em vez de minimizar uma função de custo, os métodos contornam os problemas de platôs estéreis e mínimos locais inerentes aos algoritmos variacionais.
• Compatibilidade de Hardware: A combinação de circuitos não variacionais, sem parâmetros (ou com parâmetros fixos), e a codificação eficiente de recursos do GC resulta em circuitos mais rasos com menos portas, tornando a abordagem mais robusta contra o ruído e melhor adaptada às limitações do hardware quântico atual.
• Ponte entre Teoria e Computação: O trabalho estabelece o QuGCM como um análogo quântico do GCM clássico, oferecendo um caminho para estudar correlações nucleares e fenômenos coletivos que são classicamente intratáveis, sem o ônus computacional de grandes grades de base.

Os autores concluem que, embora esses métodos sejam promissores, o trabalho futuro envolve o benchmarking em hardware quântico real e a extensão do framework para calcular observáveis adicionais, como taxas de transição eletromagnética.