Exact chiral symmetry with quantum signal processing

Resumo Técnico: Simetria Quiral Exata com Processamento de Sinal Quântico

Definição do Problema
A simulação quântica de observáveis não perturbativos e fora do equilíbrio em QCD em rede enfrenta desafios significativos em relação ao duplo de férmions (fermion doubling) e à preservação da simetria quiral em formulações hamiltonianas. Embora a relação de Ginsparg-Wilson (GW) e o operador de overlap forneçam soluções elegantes no espaço-tempo euclidiano, um análogo canônico de Hamiltonian para férmions de overlap tem sido ausente. Abordagens existentes, como os férmions de parede de domínio (domain-wall fermions), mapeiam-se bem para configurações hamiltonianas, mas exigem uma dimensão extra explícita, aumentando os custos de qubits. Por outro lado, os férmions de overlap oferecem uma formulação nas dimensões espaciais físicas, mas envolvem interações de todos-para-todos altamente não locais, especificamente a função sinal do Hamiltoniano de Wilson-Dirac, $\varepsilon(h_W)$, que é difícil de implementar eficientemente em hardware quântico. O problema central abordado é como construir algoritmos quânticos eficientes para férmions de rede com simetria quiral que equilibrem custos de memória (qubits) contra complexidade de portas, mantendo a simetria quiral exata (ou controladamente quebrada).

Metodologia
Os autores propõem um algoritmo de Processamento de Sinal Quântico (QSP) para simular o Hamiltoniano de férmion de overlap. A metodologia procede através de etapas fundamentais:

1. Formulação do Hamiltoniano: O trabalho utiliza o Hamiltoniano de overlap $H_{ov} = \psi^\dagger (\Gamma_0 + \varepsilon(h_W)) \psi$, onde $h_W$ é o Hamiltoniano de partícula única de Wilson-Dirac. A função sinal $\varepsilon(X) = X/\sqrt{X^\dagger X}$ é aproximada por um polinômio $p_M(x)$ de grau $M$.
2. Codificação de Bloco (Block Encoding): Os autores constroem uma codificação de bloco do Hamiltoniano de Wilson de partícula única $\hat{h}_W$ (incluindo campos de calibre) usando um operador "preparar" $P$ e um operador "selecionar" $S$. Isso codifica o Hamiltoniano em um operador unitário $U_{\hat{h}_W}$ agindo em um espaço de Hilbert estendido com qubits auxiliares (ancillas).
3. Processamento de Sinal Quântico (QSP): Para implementar a aproximação da função sinal, os autores empregam QSP (especificamente Transformação de Valor Singular Quântica). Isso permite a aplicação de um polinômio de grau $M$, $p_M(U_{\hat{h}_W})$, ao operador codificado em bloco. O polinômio é escolhido para aproximar a função sinal dentro de um gap espectral $[\kappa, 1]$ com um erro $\epsilon_e$.
4. Evolução Temporal: Uma vez que o Hamiltoniano de overlap seja codificado em bloco, o operador de evolução temporal $e^{-iHt}$ é aproximado usando QSP, exigindo um número de portas que escala com o custo da codificação de bloco e o tempo de evolução.

Principais Contribuições e Resultados

• Construção Algorítmica: O artigo fornece um algoritmo baseado em QSP para o Hamiltoniano de overlap que preserva a relação GW até um erro controlável $\epsilon_e$. O operador quiral modificado $\hat{\gamma}_5$ é mostrado para comutar com o Hamiltoniano aproximado até um erro de ordem $\epsilon_e$.
• Análise de Complexidade:
  • Complexidade de Portas: O custo para codificar em bloco o Hamiltoniano de overlap escala como $O(Q \kappa^{-1} \log(1/\epsilon_e))$, onde $Q$ é o número de sítios de rede (vezes os graus de liberdade internos) e $\kappa$ está relacionado ao gap espectral. A complexidade total de portas para a evolução temporal escala como $O(Q \kappa^{-1} \log(1/\epsilon_e) [Qt + \log(1/\epsilon_t)])$.
  • Custo de Qubits: O algoritmo requer $O(Q + \log Q)$ qubits. Esta é uma redução significativa comparada aos férmions de parede de domínio, que requerem $O(Q L_5)$ qubits, onde $L_5$ é a extensão da dimensão extra.
  • Comparação: Embora os férmions de parede de domínio se beneficiem de localidade geométrica (permitindo um custo quase linear $O(Q t \text{polylog})$ via Suzuki-Trotter ou métodos similares), a formulação de overlap incorre em uma profundidade de portas maior devido à natureza não local da aproximação da função sinal. No entanto, a abordagem de overlap oferece uma escalabilidade assintótica superior em termos de memória.
• Interpretação Física de QSP: Os autores demonstram que o grau do polinômio $M$ necessário para atingir um erro $\epsilon_e$ escala como $M = O(\kappa^{-1} \log(1/\epsilon_e))$. Eles identificam esta escalonagem com o tamanho da dimensão extra $L_5$ nos férmions de parede de domínio, onde o erro escala como $e^{-c L_5}$. Assim, a implementação via QSP efetivamente "constrói" uma dimensão extra através da profundidade do circuito, espelhando a correspondência física entre o operador de overlap e a teoria de fronteira dos férmions de parede de domínio.

Significância e Alegações
O artigo alega que as simulações quânticas de férmions de Dirac com simetria quiral exata são "quase gratuitas" no sentido de que aplicar o Hamiltoniano de overlap custa apenas um fator logarítmico a mais (em termos de tolerância de erro) do que o Hamiltoniano de Wilson-Dirac.

A ideia central é que o QSP fornece uma realização algorítmica-quântica da conhecida correspondência entre o operador de overlap e a dimensão extra dos férmions de parede de domínio. A troca (trade-off) é explícita:

• Férmions de parede de domínio: Maior custo de qubits ($O(Q L_5)$), mas menor profundidade de portas devido à localidade.
• Férmions de overlap (via QSP): Menor custo de qubits ($O(Q)$), mas maior profundidade de portas ($O(Q^2)$ na escalonagem de pior caso para evolução temporal devido à não localidade).

Os autores concluem que a escolha entre as formulações depende das restrições de hardware (quantidade de qubits vs. profundidade de circuito). Eles observam que, embora a QCD em rede clássica frequentemente utilize aproximações racionais (ex: Zolotarev) para a função sinal, que são mais eficientes, estas não possuem análogos diretos no QSP padrão, sugerindo uma direção para trabalhos futuros. O artigo também esclarece que os resultados de escalonamento são robustos contra a codificação específica dos campos de calibre, embora os fatores constantes possam variar.