On Performance and Limitations of NISQ Hardware for Simulations of Quantum Wave Packet Dynamics

Resumo Técnico: Sobre o Desempenho e Limitações de Hardware NISQ para Simulações da Dinâmica de Pacotes de Onda Quântica

Declaração do Problema
A simulação numérica precisa de processos dinâmicos quânticos, como vibrações moleculares, espalhamento e tunelamento, é governada pela equação de Schrödinger dependente do tempo. Abordagens clássicas enfrentam a "maldição da dimensionalidade", onde o custo computacional cresce exponencialmente com o número de graus de liberdade, tornando a propagação de longo tempo intratável para sistemas além de algumas dimensões. Embora computadores quânticos ofereçam uma plataforma natural para simular sistemas quânticos devido à sua escalabilidade exponencial do espaço de estados, a propagação explícita em tempo real de pacotes de onda usando circuitos quânticos digitais permanece pouco explorada em comparação com problemas de estrutura eletrônica. Além disso, os métodos existentes de simulação quântica frequentemente dependem de funções analíticas de energia potencial ou sofrem com escalamento exponencial na decomposição de operadores (tipicamente $\mathcal{O}(4^n)$ para $n$ qubits), limitando sua aplicabilidade ao hardware de escala intermediária ruidosa de curto prazo (NISQ).

Metodologia
Os autores propõem uma codificação baseada em grade da função de onda em um registro de $n$ qubits, onde $M = 2^n$ pontos de grade representam a coordenada espacial discretizada. A evolução temporal é implementada usando a abordagem do operador dividido, que separa os operadores de energia cinética ($\hat{T}$) e potencial ($\hat{V}$).

1. Energia Cinética: O operador cinético é aplicado usando a Transformada de Fourier Quântica (QFT). Isso mapeia a função de onda para a base de momento onde $\hat{T}$ é diagonal, permitindo uma aplicação eficiente com escalamento de profundidade de circuito polinomial de $\mathcal{O}(n^2)$.
2. Energia Potencial: Diferentemente de trabalhos anteriores restritos a potenciais analíticos, este estudo emprega uma representação geral do operador de energia potencial como uma combinação linear de cadeias de Pauli-$Z$ comutantes. Como o potencial é diagonal na base de posição, sua decomposição requer apenas $2^n$ termos (envolvendo apenas operadores $I$ e $Z$) em vez da decomposição completa de Pauli de $\mathcal{O}(4^n)$. Isso reduz o escalamento do operador para $\mathcal{O}(2^n)$ e permite a codificação direta de paisagens de energia potencial discretizadas arbitrárias.
3. Implementação: O propagador global é aproximado via uma decomposição de Trotter de primeira ordem. Simulações foram conduzidas em emuladores clássicos livres de ruído (Qiskit SDK) e executadas em hardware NISQ, especificamente processadores supercondutores da IBM (Torino, Boston, Miami) e o dispositivo de íons presos IonQ Forte.

Principais Contribuições

• Eficiência Algorítmica: O artigo demonstra um método que reduz o escalamento da decomposição do operador de energia potencial de $\mathcal{O}(4^n)$ para $\mathcal{O}(2^n)$, explorando a natureza diagonal do potencial na base de posição.
• Generalidade: A abordagem permite a simulação de potenciais discretizados arbitrários, indo além de sistemas limitados a formas funcionais analíticas conhecidas.
• Avaliação de Hardware: O estudo fornece uma avaliação abrangente da dinâmica de pacotes de onda em contagens de qubits de 2 a 5 (correspondendo a 4 a 32 pontos de grade) em múltiplas plataformas de hardware atuais.

Resultados
A metodologia foi testada em três problemas unidimensionais representativos: propagação de partícula livre, tunelamento através de uma barreira e vibrações de oscilador harmônico.

• Desempenho em Pequena Escala (2–3 Qubits): Todas as plataformas testadas (IBM e IonQ) reproduziram qualitativamente a dinâmica de referência. Para casos de 2 qubits, dispositivos IBM mais recentes (Boston, Miami) mostraram concordância quase quantitativa com o emulador.
• Escalabilidade e Sensibilidade ao Ruído: À medida que as contagens de qubits aumentaram para 4 e 5, o desempenho divergiu significativamente entre as plataformas.
  • Dispositivos Supercondutores da IBM: Os resultados mostraram degradação rápida com o aumento da contagem de qubits. Em simulações de 5 qubits, a função de onda no hardware da IBM colapsou para uma distribuição uniforme pelo terceiro passo de tempo, indicando uma perda da estrutura do pacote de onda devido a erros de porta acumulados e limitações de profundidade de circuito.
  • Dispositivo de Íons Presos IonQ: O processador IonQ Forte manteve concordância significativamente melhor com a referência clássica em todas as contagens de qubits, incluindo o caso de partícula livre de 5 qubits. Esse desempenho é atribuído a maiores fidelidades de portas de dois qubits e conectividade total a todos os qubits inerentes às arquiteturas de íons presos.
• Especificidades do Problema: Nos cenários de tunelamento e oscilador harmônico, o IonQ consistentemente superou os dispositivos supercondutores, particularmente no regime de 5 qubits, onde os dispositivos IBM exibiram desvios substanciais.

Significado e Afirmações
O artigo conclui que, embora o hardware NISQ atual possa reproduzir qualitativamente a dinâmica quântica para tamanhos de sistema pequenos, a fidelidade das simulações dependentes do tempo é altamente sensível à profundidade do circuito e ao ruído do hardware. Os resultados destacam que as arquiteturas de íons presos oferecem atualmente desempenho superior para essas simulações dinâmicas específicas em comparação com processadores supercondutores. Os autores afirmam que implementações escaláveis da dinâmica de pacotes de onda quântica exigirão designs de circuito que minimizem a profundidade e sejam robustos contra o ruído do hardware, enfatizando que a codificação de potencial proposta de $\mathcal{O}(2^n)$ é um passo necessário para utilizar potenciais arbitrários em dispositivos de curto prazo.