Quantum Simulation of Generalized Parton Distributions in the Schwinger Model

Resumo Técnico: Simulação Quântica de Distribuições de Partons Generalizadas no Modelo de Schwinger

Enunciado do Problema
As Distribuições de Partons Generalizadas (GPDs) são fundamentais para compreender a estrutura interna dos hádrons, codificando correlações entre parâmetros de impacto transversais e momento longitudinal. No entanto, o cálculo de GPDs a partir de primeiros princípios na Cromodinâmica Quântica (QCD) é dificultado pelo problema do sinal na QCD de rede euclidiana, o que impede o cálculo direto de funções de correlação de cone de luz em tempo real. Enquanto abordagens clássicas como a Teoria de Momento Efetivo de Grande Momento (LMET) oferecem contornos, elas exigem um aumento extremo de momento e sofrem com correções de potência. A computação quântica oferece uma mudança de paradigma ao permitir a simulação direta da dinâmica hamiltoniana e da evolução em tempo real, potencialmente contornando o problema do sinal. Apesar deste potencial, a simulação de GPDs permanece um desafio formidável, particularmente no que diz respeito à preservação de simetrias discretas (Conjugação de Carga $C$, Paridade $P$, Reversão de Tempo $T$) e à preparação eficiente de estados hádrons com momento não nulo.

Metodologia
Os autores propõem um algoritmo quântico abrangente para simular GPDs no modelo de Schwinger (eletrodinâmica) de (1+1) dimensões (um proxy para a QCD completa). A metodologia é distinguida por várias escolhas técnicas fundamentais:

1. Férmions de Wilson: Ao contrário das simulações quânticas atuais que utilizam predominantemente férmions staggered, este trabalho emprega férmions de Wilson para a discretização de rede. Esta escolha é crítica para a preservação estrita da simetria exata de conjugação de carga ($C$) na rede. Os autores argumentam que a simetria $C$ é vital para o estudo de mésons neutros de carga, pois elimina contaminações não físicas de parte real nas funções de correlação e garante que as GPDs permaneçam funções ímpares da fração de momento $x$.
2. Formulação Hamiltoniana: O modelo de Schwinger é formulado no calibre axial. O Hamiltoniano contínuo é discretizado usando férmions de Wilson, aumentado por um termo de Wilson para contornar o problema do duplo de férmions (fermion doubling) enquanto preserva as simetrias $C$, $P$ e $T$. O campo de gauge temporal é eliminado via lei de Gauss, resultando em um termo de interação não local no Hamiltoniano.
3. Preparação de Estados: O algoritmo utiliza um Algoritmo Evolutivo Quântico Variacional (VQE) com resolução de número quântico. Para preparar estados hádrons com momentos específicos não nulos ($|h(k)\rangle$) e o vácuo ($|\Omega\rangle$), os autores constroem um conjunto de estados de referência ortogonais. Um circuito de injeção de momento, $D(k)$, mapeia um estado de referência para um estado com momento definido. O ansatz do VQE, $U(\theta)$, é projetado para ser preservador de simetria, otimizando uma função de perda multi-estado para extrair simultaneamente o vácuo e os estados excitados hádrons dentro do setor sem carga ($Q=0$).
4. Medição de Função de Correlação: O núcleo do algoritmo avalia a função de correlação de cone de luz invariante por gauge, que inclui linhas de Wilson. Os autores implementam uma versão controlada do operador de translação de momento e links de gauge controlados ($W_0$) usando um qubit auxiliar. Ao medir os valores de expectativa off-diagonais ($\sigma_x$ e $\sigma_y$) do qubit auxiliar após a evolução em tempo real, as funções de correlação espacial $\tilde{H}_{q/h}(N, \xi)$ são extraídas.
5. Análise de Complexidade: Os requisitos de recursos são analisados usando a fórmula de produto de Suzuki-Trotter de segunda ordem para a evolução temporal. A complexidade total escala polinomialmente com o tamanho do sistema $N$ e a precisão inversa $\epsilon$, especificamente como $O(N^p + N^{5/2}/\epsilon^{1/2})$.

Principais Contribuições

• Estrutura Unificada: O artigo estabelece um arcabouço algorítmico quântico unificado capaz de computar tanto Funções de Distribuição de Partons (PDFs) quanto GPDs, em vez de tratá-las como avaliações disjuntas.
• Preservação de Simetria: Demonstra que a implementação de férmions de Wilson em um arcabouço de computação quântica não acarreta maior complexidade computacional em comparação com férmions staggered, mas oferece propriedades de simetria superiores. Especificamente, preserva rigorosamente a simetria $C$ exata, que se mostra essencial para obter resultados fisicamente significativos para mésons neutros.
• Construção Algorítmica: O trabalho detalha a construção de circuitos quânticos específicos para injeção de momento ($D(k)$), implementação de link de gauge e o ansatz VQE preservador de simetria necessário para o modelo de Schwinger.

Resultos
Os autores testaram sua abordagem usando Diagonalização Exata (ED) em uma rede espacial de tamanho $N=18$ (e até $N=22$ para verificações específicas):

• Espectro de Massa: O espectro de massa calculado para o estado ligado $q\bar{q}$ mais leve reproduziu o limite contínuo analítico conhecido ($m_h/g \approx 0.564$) via extrapolação linear, validando a formulação hamiltoniana.
• Funções de Correlação: As funções de correlação de cone de luz espacial $\tilde{H}_{q/h}(N, \xi)$ foram extraídas. Os resultados confirmaram que a parte real da função de correlação desaparece assintoticamente conforme o volume da rede aumenta, consistente com a simetria $C$.
• Comparação com Férmions Staggered: Uma comparação direta mostrou que a razão entre as partes real e imaginária da função de correlação é significativamente suprimida para férmions de Wilson em comparação com férmions staggered. Isso confirma que a formulação de Wilson elimina efetivamente contaminações não físicas decorrentes de efeitos de volume finito e violações de simetria.
• GPDs e PDFs: Transformadas de Fourier das funções de correlação produziram GPDs no espaço de momento. As PDFs resultantes (limite $\xi=0$) apresentaram pico em $x=0.5$, refletindo corretamente a natureza de dois quarks de valência do estado ligado em 1+1 dimensões. As GPDs mostraram-se funções ímpares de $x$, satisfazendo o requisito teórico de que o fator de forma de carga elástica para um estado de carga neutra deve ser nulo.

Significância e Alegações
O artigo afirma que este trabalho fornece um precursor vital para simulações de QCD em (3+1) dimensões. Ao demonstrar que os férmions de Wilson podem ser implementados eficientemente em hardware quântico sem aumentar a complexidade, os autores argumentam que esta abordagem oferece um caminho robusto para simular fenômenos não perturbativos enquanto adere estritamente às simetrias fundamentais. A escala polinomial do algoritmo com o tamanho do sistema e a precisão destaca a vantagem teórica dos computadores quânticos para a simulação direta de correlações de cone de luz. Os autores posicionam este trabalho como um passo fundamental para a generalização destes algoritmos para teorias de gauge não-abelianas (ex: Yang-Mills SU(2) ou SU(3)) e, finalmente, para a QCD completa em futuros hardwares quânticos tolerantes a falhas.