Nuclear scattering phase shifts with factorized geometric-time RODEO on a quantum processor

Resumo Técnico: Deslocamentos de Fase de Espalhamento Nuclear com RODEO de Tempo Geométrico Fatorizado em um Processador Quântico

Enunciado do Problema
Reconstruir os deslocamentos de fase de espalhamento nuclear a partir de espectros computados quânticamente apresenta um desafio significativo. Embora os resolvedores de autovalores quânticos possam estimar níveis de energia discretos em um sistema confinado (uma "armadilha"), conectar esses estados discretos a observáveis de espaço livre (espalhamento de espaço livre) requer processamento pós-clássico, especificamente a Expansão de Alcance Efetivo Modificada (MERE). O desempenho de um circuito quântico neste contexto não pode ser julgado apenas pelo erro quadrático médio (RMS) de níveis de energia individuais. Erros espectrais podem ser amplificados, suprimidos ou reorganizados durante a subsequente interpolação de energia e extrapolação de confinamento zero. Portanto, um circuito quântico com redução de recursos deve ser avaliado por sua capacidade de preservar a informação necessária para a reconstrução final do espalhamento, em vez de apenas pela fidelidade espectral.

Metodologia
Os autores reconstroem os deslocamentos de fase elástica s-wave de nêutron–próton (np) utilizando dados do processador quântico IBM Aachen. O estudo compara duas implementações do algoritmo RODEO (Repetitive Operation for Discrete Energy Optimization):

1. Direct R10 Dinâmico: Uma implementação de base utilizando dez ciclos com passos de tempo de quantil Gaussiano fixos. Para cada energia de consulta, a evolução controlada é ressintetizada e compilada separadamente.
2. FG–R6 de Tempo Geométrico Fatorizado (Factorized Geometric-Time RODEO): Uma implementação comprimida utilizando seis ciclos com passos de tempo amostrados geometricamente. Este método emprega a separação exata da fase de consulta, onde a evolução do Hamiltoniano é fixa a uma energia de referência ($E_{ref}$) e compilada uma única vez. A fase dependente da consulta ($\Delta E = E_q - E_{ref}$) é vinculada ao qubit ancila após a compilação.

Ambos os métodos utilizam um registrador de três qubits (dois qubits de sistema e uma ancila reciclada) no backend IBM Aachen. O modelo físico é uma interação de poço quadrado s-wave np esquemática. Hamiltonianos efetivos Okubo–Lee–Suzuki (OLS) de quatro dimensões, construídos classicamente, são usados para representar o sistema. O experimento mede quatro níveis de energia positiva em cinco intensidades de armadilha diferentes ($\omega = 3,6$ a $4,0$ MeV), gerando vinte entradas espectrais. Essas entradas são alimentadas em uma análise MERE clássica para extrapolar para o espaço livre e calcular os deslocamentos de fase sobre uma grade de energia de 0,1–30,0 MeV.

Principais Contribuições

• Compressão de Circuito via Fatorização: O artigo introduz o circuito FG–R6, que combina o reuso de ancila com uma estratégia de evolução temporal fatorizada. Ao separar a fase de consulta da evolução do Hamiltoniano, o circuito reduz o número de ciclos de evolução controlada de dez para seis e elimina a necessidade de recompilar o operador de evolução para cada energia de consulta.
• Benchmarking de Hardware: O estudo fornece um benchmark completo derivado de hardware para espalhamento nuclear, indo além de testes de nível único para um conjunto completo de vinte entradas espectrais necessárias para a reconstrução de espaço livre.
• Análise de Troca entre Recurso vs. Precisão: O trabalho demonstra que reduzir a profundidade do circuito e a contagem de portas de dois qubits não degrada necessariamente a precisão do observável de espalhamento final, mesmo que o erro RMS dos níveis de energia intermediários aumente.

Resultados

• Redução de Recursos: Em relação à implementação Direct R10 Dinâmica, a implementação FG–R6 reduz a profundidade média do circuito compilado em 33,1% (de 459 para 307) e a contagem de portas de dois qubits em 37,2% (de 113 para 71).
• Incerteza de Energia: O FG–R6 alcança uma incerteza de energia de finite-shot mediana menor (0,534 keV) comparado ao R10 (0,610 keV), apesar de possuir uma amplitude ajustada mediana maior e um deslocamento central RMS ($\Delta E_c$) maior em relação à energia de referência (1,648 keV vs. 0,793 keV).
• Precisão do Deslocamento de Fase de Espalhamento: Na grade completa de 0,1–30,0 MeV, o desvio máximo do deslocamento de fase reconstruído em relação à referência MERE de energia exata diminui significativamente de 0,842° (R10) para 0,163° (FG–R6). O desvio máximo em relação à solução analítica do poço quadrado também cai de 0,981° para 0,156°.
• Desempenho Dependente de Energia: A melhoria é impulsionada principalmente pela região de extrapolação de baixa energia (abaixo de 10 MeV). No subconjunto de interpolação comum de 10,0–30,0 MeV, a implementação Direct R10 produz, na verdade, resíduos ligeiramente menores que o FG–R6.
• Sensibilidade da Forma do Ajuste: O estudo observa que a sensibilidade da forma do ajuste permanece apreciável. A dispersão nos resultados entre diferentes ordens polinomiais para os ajustes de energia e armadilha é substancial (por exemplo, uma meia-amplitude de 20,12° na grade total para R10 vs. 6,75° para FG–R6), indicando que a escolha do modelo de extrapolação impacta significativamente o resultado final.

Significância e Alegações
O artigo afirma que este benchmark de espaço reduzido conecta com sucesso a compressão de circuito RODEO aos observáveis do contínuo nuclear. A principal significância reside em demonstrar que o desempenho do circuito deve ser avaliado através da reconstrução do observável de espalhamento e seu intervalo de energia específico, em vez de apenas pelos erros de RMS espectrais. A implementação FG–R6 é apresentada como um compromisso de precisão–recurso validado para o espectro e janela específicos testados, e não como uma solução de ciclo mínimo universal.

Os autores declaram explicitamente que estes resultados não demonstram vantagem quântica ou escalabilidade de muitos corpos. Os dados de hardware não determinam separadamente os efeitos da compilação de fase versus amostragem de tempo, nem estabelecem a reprodutibilidade de execuções independentes. Este trabalho serve como um ponto de partida para estender estes métodos para espaços de modelos maiores, interações realistas e dinâmicas de reação de poucos corpos, tais como o espalhamento ressonante de ondas parciais superiores (por exemplo, espalhamento nêutron–alfa) e sistemas envolvendo interações de Coulomb de longo alcance.